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Abstract 



During the last years quantum graphs have become a paradigm of quantum 
chaos with appUcations from spectral statistics to chaotic scattering and 
wave function statistics. In the first part of this review we give a detailed 
introduction to the spectral theory of quantum graphs and discuss exact 
trace formulae for the spectrum and the quantum-to-classical correspon- 
dence. The second part of this review is devoted to the spectral statistics 
of quantum graphs as an application to quantum chaos. Especially, we 
summarise recent developments on the spectral statistics of generic large 
quantum graphs based on two approaches: the periodic-orbit approach and 
the supersymmetry approach. The latter provides a condition and a proof 
for universal spectral statistics as predicted by random-matrix theory. 
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Chapter 1 

Introduction 



The general mathematical concept of a graph (network) as a set of ele- 
ments which are connected by some relation has found applications in many 
branches of science, engineering, and also social science. A street network 
of a traffic engineer, the network of neurons studied by a neuroscientist, the 
structure of databases in computer science can all be described by graphs. 

Recently the Laplacian on a metric graph has gained a lot of attention in 
physics and mathematics in terms of the diffusion equation or Schrodinger 
equation. They have now become known as quantum graphs but differ- 
ent aspects are studied under various names such as quantum networks or 
quantum wires. They have a long history in mathematics and physics. 

In physics, the first application has probably been in the context of free 
electron models for organic molecules about seventy years ago by Pauling 
[198J, an approach which has been further developed in subsequent years 
[172t ESI [Ml [Ml [73 IIS21 [M] • Quantum graphs have also been applied 
successfully to superconductivity in granular and artificial materials fl2\ . 
acoustic and electromagnetic waveguide networks [1V2\ I181| , the Anderson 
transition in a disordered wire [IHl I222j , quantum Hall systems [TH I15H 1150] , 
fracton excitations in fractal structures \22\ I190j , and mesoscopic quantum 
systems |142^ 11661 1242^ [2i3| 1241 j . Quantum graphs have also been simulated 
experimentally |141j . 

The construction of self-adjoint operators, or wave equations with ap- 
propriate boundary conditions on graphs has first been addressed by Rue- 
denberg and Scherr [209j (see also [204j ). They considered graphs as an 
idealisation of networks of wires or wave guides of finite cross-section in the 
limit where the diameter of the wire is much smaller than any other length 
scale. Similar approaches to graphs as networks of thin wires with a finite 
diameter, or fat quantum graphs as they are now called, have been a topic 
in mathematical physics recently [2111 [2l21 [2081 [Ml [Ml [M]- 

Another interesting aproach which has been discussed mainly by Exner 
and his coworkers is based on leaky graphs [Ml [Ml [Ml [Ml [Ml [Ml 
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11071 11081 1109j . Here, a finite attractive potential in the Schrodinger equation 
is centered on a metric graph. A leaky graph is a generalisation of the 
Schrodinger equation with (5-function potentials. A quantum graph can be 
realised in a limit of infinitely strong attracting potentials of this type. 

Probably the first mathematical approaches to the Laplacian on a metric 
graph were by Roth who derived a trace formula for the spectrum of the 
Laplacian [2071 EQG] and by von Below [30l|3ll|33]. The spectral theory of 
quantum graphs was mainly developed on the basis of the von Neumann 
theory of self-adjoint extensions for formal differential operators [241 [HI [92| 
[1921 [971 [951 IMl [Ml [13 [Ml [Sa [Ml IMl [15^ other 
recent topics in mathematics and mathematical physics include the spectral 
theory of infinite periodic graphs [Ml [Ml IM HOI [93 HZl [IH [IMl [Ml [M] , 
tree graphs ^61 [3 [3 [3 [23ll M\ [189] and Sierpinski graphs [233 1123 [223] . 
scattering and bound states in open graphs [27 1 [Ml[96l[98l [T2Tl[TM[239l[2l0] . 
diffusion and localisation [32l [SI [83 IMl IMl [85], random walks [22011156]. 
approximations of quantum systems by quantum graphs [111] , some inverse 

problems [Ml [Ml [13 113 [Ml [Ml [IMl [Ml [IS3 [IZ3 El [Ml [Ml [IMl [Ml 

I122j . extremal spectral properties [115] . and models of dissipative graphs 
[Ml [Ml IMl [M]- 

Recent review papers on various mathematical approaches to quantum 
graphs can be found in [1681 11691 1170] . Special issues of research journals 
[1711 [ST] and recent conference proceedings [17] have been devoted to quan- 
tum graphs. 

The relevance of quantum graphs to the study of quantum chaos was 
brought to light by the work of Kottos and Smilansky [15911160] . They anal- 
ysed the spectral statistics for simple graphs, and showed that their spectral 
statistics follow very closely the predictions of random-matrix theory. They 
proposed an alternative derivation of the trace formula and pointed out 
its similarity to the famous Gutzwiller trace formula [1321 1133] for chaotic 
Hamiltonian systems. While quantum graphs do not have a deterministic 
classical limit they still share a lot of important properties with classically 
chaotic Hamiltonian systems, e.g. periodic-orbit theory in the semiclassi- 
cal regime is completely analogous to periodic-orbit theory on a quantum 
graph. While semiclassical approaches are approximations the analogous 
approaches for quantum graphs are exact. More importantly, quantum 
graphs are not as resistant to analytical approaches. Following these pi- 
oneering beginnings of Kottos and Smilansky quantum graphs have become 
a new paradigm of quantum chaos and have been applied to various prob- 
lems, including also disorder and diffusion |251[[75] . All aspects of quantum 
chaos have been covered: spectral statistics in finite [23 [33 [13 [13 [Ml [HI 

11^ [^ [^ [TH^ [TM inm lim [^ imi 1^ [^ 1^ [T^ [T?n [T^ 

and infinite periodic structures [85], localisation and wavefunction statistics 
[33 [MSI [211 [IMl [na [2111 [13 [13 [251], chaotic scattering [IMl [113], trans- 
port through chaotic devices [2171 174]. resonances and decay in open chaotic 
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systems p03l [TM [T65] . 

The main purpose of tiiis paper is to review recent approaches to spectral 
statistics in quantum graphs and discuss the relation to quantum chaos 
in general. To make the paper self-contained special effort was made to 
introduce the necessary concepts in a clear and consistent way. 

Recently important progress has been achieved by Miiller et al I186j 
in the periodic-orbit theory for spectral statistics in chaotic Hamiltonian 
systems building on the pioneering work of Sieber and Richter }2U5l l225|l226j . 
A similar approach for quantum graphs has been developed by Berkolaiko, 
Schanz and Whitney |42l [HI I45| . None of these works captured the low 
energy sector of spectral statistics. An alternative approach which so far 
was successfully applied to quantum graphs \127\ 1128] was able to close the 
above mentioned gap, and it will be discussed in detail in the second part 
of this review. 

Several other applications to quantum chaos will not be addressed in 
this review. Functional integration approaches to statistical properties of 
wave functions and localisation are thoroughly reviewed by Comtet, Desbois 
and Texier [76] (see also [7j). Chaotic scattering and transport in quantum 
graphs is discussed in [Ml [Ml [Ml I2IZI Ml- 

This review is arranged in the following way. The first chapters provide a 
detailed introduction to quantum graphs and their spectral theory. Though 
the choice of material is biased by the application to quantum chaos the first 
chapters are rather general. In 11.11 we define graphs and their topological 
description. In Chapter[2]we quantise graphs in a straight forward way. The 
reformulation in Chapter [3] in terms of a unitary quantum evolution map 
will allow us to generalise the quantisation procedure. This point of view 
will be the foundation for the remaining chapters. Chapter [H gives account 
of the classical dynamics that corresponds to a quantum graph and Chapter 
[5| is devoted to spectral theory with the introduction of the trace formula 
and its discussion. 

The remaining chapters are devoted to spectral statistics in quantum 
graphs. After a general introduction of spectral correlation functions (and 
other statistical measures of the spectrum) in Chapter [6] we will give some 
background from quantum chaos and random-matrix theory on universal- 
ity in the spectral fluctuations of complex quantum systems in Chapter [71 
In the last two chapters we present two analytical approaches to universal 
spectral statistics in large quantum graphs. The periodic-orbit approach is 
presented in Chapter [8] which also includes sections on quantum graphs with 
spin and quantum graphs which are coupled to superconductors (Andreev 
graphs). Finally, we summarise the supersymmetry approach with a proof 
of universality in large quantum graphs in Chapter [3 In Appendix [X] we 
add some background on the symmetry classification of quantum systems 
and in Appendix[B]we summarise some relevant results from random-matrix 
theory. 
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1.1 Graphs and their topology 

A graph Q{V,B) consists of V vertices connected by B bonds (or edges). 
A graph with six vertices and ten bonds (V = 6, B = 10) is shown in figure 





Figure 1.1: A graph with V = Q vertices and i? = 10 bonds. 

The graphs are not necessarily embedded in the plane, and the fact 
that in the figure bonds cross each other at points which are not vertices 
is completely immaterial. A physical realisation of a graph is a network of 
coaxial cables (bonds) connected by junctions (vertices). The topology of 
the graph, that is, the way the vertices and bonds are connected is given 
in terms of the V connectivity matrix Cij (sometimes referred to as 
the adjacency matrix) which is defined as 



in if z 7^ j where i and j are connected by m bonds, 
2m if z = j and there are m loops at vertex i 
if i and j are not connected. 

(1.1) 

This definition allows for vertices to be connected by several bonds, and also 
for a vertex to be connected to itself by one or several loops (in which case 
the corresponding diagonal element of the connectivity matrix equals twice 
the number of loops). 

The valency Vi of a vertex i is the number of vertices j connected to 
i, each weighted by the number of parallel bonds (loops). Thus, m parallel 
bonds (m parallel loops) contribute m (2m) to the valency. In terms of the 
connectivity matrix, 

V 

Vi = 

The neighbourhood Fj of the vertex i consists of the vertices j con- 
nected to i. 

The boundary of a subgraph Q C G, F(^), consists of the vertices which 
are not in G but which are in the union of the neighbourhoods of the vertices 
of 
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The number of bonds is expressed by 




(1.3) 



Unless otherwise specified, we shall always consider connected graphs, 
for which the vertices cannot be divided into two non-empty subsets such 
that there is no bond connecting the two subsets. That is, for a connected 
graph the connectivity matrix cannot be brought into a block-diagonal form 
by permuting the vertices. 

There are a few classes of graphs which often appear in the literature. 
They are characterised by their connectivity (see figure 11.21 for some exam- 
ples): 

• Simple graphs are graphs which have no loops and no parallel bonds 
connecting their vertices (no multiply connected vertices). In this 
case, for all i and j, Cij £ {0, 1}, and in particular all the diagonal 
elements vanish Q^, = 0. For simple graphs, the cardinality of T, is 
the valency Vi for each vertex. When we define quantum graphs in 
chapter [2] we will show that every connected quantum graph can be 
turned into a graph of simple topology by adding some vertices without 
changing the spectrum or the wave functions. This will allow us to 
significantly simplify the notation in the remaining chapters where a 
simple topology will be assumed without loss of generality. 

• V -regular graphs are simple graphs whose vertices have the same 
valency v. The simplest v regular graphs are the rings for which 
V = 2 and V = B. A non-trivial ring has at least two vertices. A 
V -regular graph is complete when v = V — 1. 

• Simply connected graphs do not contain any nontrivial ring as a 
subgraph. 

• Tree graphs are simple, connected and simply connected. 

• Star graphs are trees that consist of a main (central) vertex with 
valency v, connected to v peripheral vertices of valency one. 

In many applications it is convenient to refer to bonds directly, and we 
shall use lowercase letters to denote the bonds of the graph. If the graph is 
simple, we can use the end points of the bond as its label : b = = (j, i). 
If a graph is not simple (i,j) will denote the set of all bonds that connect 
the vertices i and j. The round brackets will always be used to denote such 
a set of undirected bonds unless otherwise stated. 
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(d) (e) 



Figure 1.2: Some examples of graphs: (a) star graph (B = W, V = 11), 
(b) ring graph {B = 10, V = 10), (c) v-regular graph with v = A {B = 20, 
F = 10), (d) complete graph {B = 45, F = 10), (e) tree graph {B = 19, 
V = 20). 

All the bonds which emerge from a vertex i form a star, 

5»= U (1-4) 

If the graph is simple, the bonds in a star 5^*) are {b = : j G r(i)}. 

Directed bonds (also referred to as arcs in the literature) are bonds 
on which a direction is specified. For a simple graph we denote them by 
the ordered pair of vertex indices enclosed in square brackets [j, i] , and the 
direction is from the right to the left index. Again, for non-simple graphs 
[j, i] is the set of all directed bonds starting at ?' and ending at j. Lowercase 
Greek letters will be used to distinguish them from the undirected bonds. 
The reverse direction will be denoted by a hat so that e.g., if a G [j, i] then 

Alternatively, we will denote a directed bond a G \j, i] as a pair a = {b, lo) 
of a bond b G (i, i) and a direction index lo = ±1 where uj = +1 if j > i 
and = — 1 if j < i. If the bond is a loop (i = j) the direction ui has to be 
assigned to avoid ambiguities. We will use the notation uja [ba] to refer to 
the direction index (bond) of the directed bond a. If a = (6, to) is a directed 
bond the reverse direction is a = (6, —lj). All directed bonds that start at 
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a vertex i form an outgoing star 



S'. 



+ 



U b''^]- 



(1.5) 



ier(i) 



The incoming star S_ is defined analogously as the set of directed bonds 
that end at i. 

The set of directed bonds [k, I] follows the set of directed bonds [r, s] 
if r = Z. With a G [k,r] and /? G [r, s] we then write a G ^riP) which 
means that the directed bond a follows f3 at the vertex r. A trajectory t = 
{ai, . . . ,an) from vertex i to j is a sequence of directed bonds ai = {huuji) 
such that a;/+i follows ai where ai starts at vertex i and an ends at vertex j. 
The topological length of the trajectory is the number n G N of directed 
bonds in the path. A closed trajectory starts and ends at the same vertex 
i = j and a periodic orbit p = ai, . . . , a„ of period n is the equivalence 
class of closed trajectories that are equal up to cyclic permutation. The code 
of the periodic orbit is the equivalence class ai, . . . , of visited directed 
bonds. A primitive periodic orbit has a code which cannot be written 
as a repetition of a shorter code. Each trajectory t defines a subgraph Qt 
which consists of all bonds and vertices visited by the path. The number 
of different points on periodic periodic orbits with period n is exactly equal 
trC". As a consequence, if n is prime the number of periodic orbits #(n) 
of period n is exactly 



where 7j are the eigenvalues of the connectivity matrix. If n is not prime 
(ll.6p is a good approximation. In the limit of large n, #(n) is dominated by 
the maximum eigenvalue which shows that this number grows exponentially 
with n. In analogy to the theory of dynamical systems we define the topo- 
logical entropy At as the logarithm of the largest eigenvalue. One can 
easily show that for fully connected graphs At = log^ and for star graphs 




(1.6) 



At = 2 




Chapter 2 

The Schrodinger operator on 
graphs 



In chapter [T] we defined and discussed the graphs from a topological point of 
view, where concepts like connectivity and neighbourhood played the main 
roles. At this point we would like to endow the graphs with a metric which 
will enable us to define the Schrodinger operator on the graph. 

We assign the natural metric to the bonds. The position a; of a point 
on the graph is determined by specifying on which bond b it is, and its 
distance Xb from the vertex with the smaller index such that Xb increases in 
direction u; = +1 and decreases in direction iv = —1. If the bond 6 is a loop 
this defines the starting point for x^- The length of a bond is denoted by 
Lb and, < x;, < Lb- The length of a path t = {ai, . . . , On) is the sum over 
all bond lengths Lt = ^^1=1 (where ai = {bi,uJi)) along the path. For a 
bond b = we will also use the notation x^ and for the values of Xb 
at the vertices. That is, if i < j one has = and xj^ = Lb- For later use, 
we define incommensurable (or rationally independent) bond lengths 
on the graph. For these the equation 



where nib £ ^ only has the trivial solution nib = for all b. 

The Schrodinger operator on Q consists of the one dimensional operators 
associated with each bond: 



Here, Wb{xb) is a potential function assumed to be non negative and smooth 
on the interval [0,Lfe]. Ab are real, positive constants. If the graph contains 
a nontrivial ring as a subgraph (that is, it is not simply connected), and Ab 




(2.1) 



b 




(2.2) 
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2. The Schrodinger operator on graphs 



do not vanish on this subgraph, time reversal invariance is broken, as can 
be seen from the eigenfunctions of (12. 2p with Wb = on the hne. 



(ci and C2 are arbitrary constants). The complex conjugate of iIj{x) above 
is not a solution of the same equation - a hallmark of a system which vio- 
lates time reversal invariance. We shall refer to the constants Ab as mag- 
netic fluxes because they play the same role as the vector potentials in the 
Schrodinger operator in higher dimensions. The presence of A implies also 
that the wave functions are intrinsically not symmetric under the reflection 
X ^ —X. 

In the physics literature, and in particular, in the quantum chaos con- 
nection the bond potentials Wb{xb) are commonly set to zero. We included 
them in the general framework for the sake of generality, and also because 
they appear in some of the mathematical literature on the topic. Graphs 
with non vanishing potentials are sometimes referred to as dressed graphs 



Next, we have to identify the space of wave functions and the boundary 
conditions which renders the operator self-adjoint. In physical terms, this 
implies that the evolution induced by the operator conserves probability and 
the vertices cannot be either sinks or sources. In other words, the boundary 
conditions at the vertices should be such that the total probability current 
vanishes when summed over all the bonds which emerge from any of the 
vertices. This is similar to the well known Kirchoff rule in the theory of 
electric networks, and we shall now derive it for quantum graphs. 



2.1 Vertex boundary conditions — self-adjoint ex- 
tension 

We consider the set of functions T> which have the following properties: 
^{x) G T) are continuous and complex valued functions ol x & Q, with 
^{x) = ipbi^b) for X £ b, and < Xf, < Lb- Continuity is implied here also at 
the vertices. This means that at each vertex i the limit lim^^^^i^ ipb{xb) = (t^i 

does not depend on 6 e 5*^*). The functions il'bixb) are complex valued, 
bounded with piecewise continuous and square integrable first derivatives. 





(2.3) 



[HolEaiMlEIS]. 
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The set P is the domain of the positive definite quadratic form 



QAm = [ dx 
Jg 



V + ^ "ffx 




h Abi^b{Xb) 

I axb 



Wb{Xh)\ll^b{Xb)\^ + 



(2.4) 



Here, cj) denotes the V dimensional array of 4>j. A is a positive diagonal 
matrix with elements \j > 0. The role of these parameters will become 
clear in the sequel and their physical significance will be discussed below. 

The self-adjoint extension of the Schrodinger operator, is determined 
by the Rayleigh-Ritz extremum principle. For this purpose we compute 
the variation of the quadratic form (5a[^] with respect to ^ under the 
condition that = Jg<ix\'^{x)\'^ = 1. This constrain is introduced by 

considering the variation of the modified quadratic form Q[^] = Qa[^] — 
/c'^ll^'ll, where k"^ is a positive Lagrange multiplier. The variations of the 
modified quadratic form with respect to both ^' and its complex conjugate 
have to vanish identically. Writing the variation with respect to ^* explicitly, 
and performing partial integration where necessary, we get 



5QA,m = V / ' dxbSrbixb) 

6=1 



1 d 

i 6xb 



+ Ab] + Wb{xb) - k"^ ipb{xb) 



+ 



Hbi^b) + ) ^b{xb) 



+ Y6(t)*\j<i)j 



(2.5) 



A similar expression can be obtained for the variation of Q with respect to 
^{x). Requiring that both variations vanish for every (5\I'(x) and 5'^*{x), we 
find that the domain Vh of the Schrodinger operator consists of functions 
in V, with twice differentiable il)b{xb)-, which satisfy the boundary conditions 



y l<i<V : 



V uj(^ + \Ai 

^ \OXb 



Ipb 



(2.6) 



The derivatives in (12. 6p above are computed at the common vertex i. That 



is at Xfe = = if w = +1 and Xb = 
loop, both directions appear in the sum) 



xl = Lb uj = —I (if 6 is a 
These conditions are obtained 
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by reorganising the second line in (j2.5p according to the vertices and their 
associated (outgoing) stars. 

The eigenfunction are solutions of the bond Schrodinger equations 



which satisfy the boundary conditions (12. 6|) . The spectrum {/c^}^^ is dis- 
crete, non-negative and unbound. It consists of the values of the Lagrange 
multiplier fc^ for which a non-trivial solution is found. The sequence of 
eigenvalues is conveniently arranged by increasing order, so that fc„ < if 
n < m. 

The eigenfunctions (ordered by increasing eigenvalues) have the following 
property. Let T>n denote the subspace of functions in D which are orthogonal 
to the first n — 1 eigenfunctions of H. Then, for any non-zero ^ £ Vn 



Equality holds if and only if $ is the n'th eigenfunction of H. 

At this point it will be good to familiarise oneself with the general for- 
malism above by considering a few examples. 

The boundary conditions (12. 6p take a very simple form for vertices with 
valency v = 2, when the magnetic fluxes on the two bonds are the same 
and the potentials on the two bonds take the same value at the vertex. We 
can now think of the two bonds as two adjacent intervals on the line with 
a common boundary which we choose as the point x = 0. Writing (|2.6p 
explicitly we get 



This is the well known boundary condition one obtains from a "(5-potential" 
of strength A. In physics textbooks this boundary condition is derived by 
integrating the Schrodinger equation along an interval of size 2e centred at 
x = 0, and the self-adjoint extension of the operator is obtained automati- 
cally. This is not the case when we are dealing with a more complex vertex, 
and the road we chose to derive the boundary conditions must be taken. 

Even though the parameters Aj can take arbitrary values, the limiting 
values Aj = {Neumann boundary conditions) or 1/Aj = {Dirichlet 
boundary conditions) are of special interest. 

The importance of Neumann boundary conditions (also known as Kirch- 
hoff boundary conditions in the mathematical literature) comes from 
the fact (to be proven shortly) that the spectra of systems with finite (but 
not vanishing) Aj approach to the Neumann spectrum as one looks higher 
up in the spectrum. A similar situation is well known for the spectra of 
Schrodinger operators on domains with boundaries, where the boundary 



\/ l<b< B : Hbi^b = k^^h 



(2.7) 



(2.8) 




(2.9) 
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conditions intermediate between Dirichlet and Neumann are studied |224j . 
Note that the case studied in the example above of a vertex with f = 2 is 
trivial under the Neumann condition. Indeed, (j2.9p impHes that the wave 
function and its first derivative are continuous so that the point x = be- 
comes an ordinary point on the interval. Thus, a Neumann vertex with 
V = 2 can simply be erased from the graph, without any effect on the spec- 
trum or wave functions. This property can be used in the reverse direction 
as well. In particular, every non simple quantum graph can be turned into 
an equivalent (with respect to spectra and wave functions) quantum graph 
of simple topology (no loops, no multiple connections) by adding two Neu- 
mann vertices on each loop and one Neumann vertex on each bond which 
is responsible for multiple connectivity. Because of this reason, and un- 
less stated otherwise, we shall henceforth consider only graphs with simple 
topology: 

dj E {0, 1} and d^i = . (2.10) 

This significantly simplifies the notation while the transition back to the 
original connectivity is straight forward. 

The Dirichlet boundary conditions imply that the value of the wave 
function vanishes at all the vertices. This isolates the various bonds and 
the spectrum of the Schrodinger of the graph reduces to the union of the 
bond spectra, with Dirichlet boundary conditions on each of the ends of the 
bonds. The wave function corresponding to a nondegenerate eigenvalue k"^ 
is identically zero on all bonds but one and the spectrum is given by the 
union of the independent spectra on the bonds 



^Dirichlet — 




(2.11) 



2.2 The secular equation 

In the previous section the Schrodinger operator and the accompanying 
boundary conditions were defined and discussed. Here, we shall derive the 
secular function - the real valued function whose zeros (the values of the 
argument where the function vanishes) stand in one to one correspondence 
with the spectrum of the graph. The spectrum of Dirichlet graphs was 
computed previously and we shall exclude this case in the following. 

We shall first assume that the bond potentials vanish, wi,{xb) = 0. We 
will comment on non-vanishing bond potentials at the end of this section. 

The starting point is the potential-free solution (j2.3|) . On a bond b = 
with i < j, the general solution of the bond Schrodinger equation is a 
superposition of the two solutions (j2.3|) . which can be written as 

Mxb]k) = ———{^ je'^^^^ sin kxb + (pi sin k{Lb-Xb)) . (2.12) 
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At the two endpoints X{, = and Xf, = L^, the wave function ijjh{xb;k) 
assumes the values (pi and cpj, respectively. These same values appear in 
the wave functions on bonds which connect i or j to other vertices, and 
therefore, the resulting graph wave function is continuous. Thus, the most 
general form of ^' G T>h is completely determined (up to a scalar factor) by V 
complex parameters = (</>i, • • • , (pv)- The appearance in the denominator 
of sin kLh for all b will eventually lead to poles of the secular function at the 
Dirichlet spectrum. 

The eigenfunctions and eigenvalues of the Schrodinger operator are found 
by demanding that the B functions (j2.12p satisfy the boundary conditions 
(|2.6p . Upon substitution, we obtain a set of V homogeneous and linear 
equations for the vertex wave functions 0, 

V 

Y,(Pjhj4k) = yi = l,...,V, (2.13) 

J=l 

where, 

hi,iik) =^ + ^ cot kLh 

{b,u)GS^;^ (2.14) 

gi uiAbLt 

hj,i{k) = - Cj^i . , {h,uj) = . 

Sin ixiIj^ 

Note that the matrix h is hermitian because the indicators lo change their 
sign when i and j are interchanged. Equations ()2.13p have a non-trivial 
solution if and only if 

Ch{k) = det h{k) = 0. (2.15) 

The function C,h{k) is the secular function. Because of the symmetry 
'tpb{x;k) = —ip[i{x,—k), the zeros of the secular function appear symmet- 
rically on the negative and positive k half-lines at =b/i;„. The Schrodinger 
spectrum is A;^. In Chapter O an alternative secular function will be de- 
rived. The suffix h in C,h{k) is added to distinguish between this and the 
other secular function. 

In the previous section we commented (without proof) that in the limit 
of large eigenvalues, the spectra of the Schrodinger operators with Ab 7^ 
converge to the Neumann spectrum. This is immediately apparent from 
(I2.14p . since the parameters A;, appear in the combination ^ which vanishes 
in the limit of large k. 

Finally, it will be instructive to write explicitly the secular function for 
a simple example. We do it for a star graph which was defined in Section 
11.11 We denote the central vertex by the index 0, and label the emanating 
bonds as well as their exterior vertices by i, i = 1, • • • The only non 
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vanishing matrix elements of h are 



^0,0 



- — h / cot kLi 




Aj 

-p- + cot kLj 

k 



(2.16) 




sin kLj 



The secular function can be computed directly, 



Ch{k) = ho,o - J2 




11": 



B 




(2.17) 



It is worth while noticing that the magnetic fluxes Aj do not appear in 
the secular function. The reason for this is that the star graph is simply 
connected. Moreover, canonlyvanishifthefirstfactor ontherighthandside- 
ofequation (I2.17p vanishes. The secular equation reduces to 



To end this section we shall indicate how the computation of the secular 
equation can be generalised to include non- vanishing potentials Wb{xb) on the 
bonds. In the general case one cannot write an explicit solution of the bond 
wave equation as in ()2.12p . However, since it is of the Sturm-Liouville type, 
it is possible to find two independent solutions, which satisfy independent 
initial conditions. One function vanishes at the vertex i and its derivative is 
set to 1, and the other vanishes at the vertex j and its derivative is 1. The 
most general solution which takes the values and (j)j at both ends of the 
bond can be written as a linear combination of the two solutions, and from 
here on one can proceed as in the potential free case. Use should be made 
of the Wronskian relation, which in the presence of the magnetic flux takes 
the form 





(2.19) 



Chapter 3 

The Quantum Evolution 
Map 

It is quite common in theoretical physics that a given subject can be for- 
mulated and studied from several points of view which use different con- 
cepts and tools. Consider for instance, the Lagrangian and Hamiltonian 
approaches to classical mechanics. The phase-space description reveals the 
symplectic structure of classical dynamics in a natural way, which is only 
implicit in the configuration-space description. Following this example, we 
shall present in this chapter another formulation of the spectral theory of 
graphs which uses other concepts and structures then those used in the pre- 
vious chapter. The main observation and intuition which underlays the new 
formulation is that a wave function on the graph can be written as a super- 
position of waves travelling in opposite directions on the bonds. The waves 
which propagate towards a given vertex are scattered from the vertex and 
emerge as outgoing waves, which scatter again and again. A wave function 
is an eigenfunction if it is stationary under the multiple scattering scenario 
described above, and this, in turn can happen only if the wave number k 
of the propagating wave is correctly selected. Thus, a new form of the sec- 
ular equation can be obtained, based on this multiple scattering approach 
[160j. The advantages gained by developing the alternative theory are both 
conceptual and technical, and they will be revealed as the theory is system- 
atically unfolded in the following sections. These advantages are gained at 
a price: the number of parameters necessary to represent a wave function 
increases from V to 2B, which is the number of amplitudes of waves which 
propagate back and forth on the bonds. 

Note: In the present chapter we consider the "bare" graphs so that on 
all the bonds Wf, = 0. We also want to remind that we assume a simple 
topology without loss of generality (see Section f2.ip . 
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3.1 Vertices as scattering centres 

Given a graph Q, the bond lengths Lj,, magnetic fluxes Af, and the boundary 
condition parameters Af,, one can write a general solution of the Schrodinger 
equation ()2.12p (for i < j) 

ipbixb] k) = . Uje' ^ sin kxb + sin k{Lb - Xb)) , , 
svnkLb (3.1) 

where, aj^^j = a" and a°j^*j = a|^"* are the (complex valued) amplitudes of 
the counter propagating waves (incoming and outgoing at vertex i) along 
the bond 6, and they can be readily written in terms of the parameters 
and In Section [LD we defined the outgoing (incoming) stars (S'i*'') as 
the set of directed bonds which point away (point at) the vertex i which are 
convenient for the following discussion. Consider now a single vertex i and 
the bonds h G 5^*^ For the present discussion we also use the convention 
that the vertex i corresponds to the point Xf, = for all the bonds in the 
star 5^. From we see that (/>j can be expressed as 

'A. = C + 4%, (3.2) 

on all the bonds h G 5^*^. These relations provide Vi — \ homogeneous 
equations 

C + = <t + for b, b' G 5« . (3.3) 



Another homogeneous expression is derived by substitute both (13. ip and 
()3.2p in the boundary condition at vertex i (j2.6p . 



These Vi relations enable us to write the outgoing amplitudes a°"*, a G S^^ 
in terms of the incoming amplitudes a™, /3 G S*^^ in the form 



where, 



and, 



-a = - (3-6) 
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This can be written shortly as a'^*^ = ^,^^^|^. IE — 1 where IE is a full matrix 
with unit entries. 

It is not difficult to prove that the Vi x Vi symmetric matrix a^^^ is 
unitary, which implies that the total outgoing probability current equals the 
total incoming probability current, as stated (without proof) in a previous 
chapter. The matrices cr^*-' will be referred to as the vertex scattering 
matrices. The simple topology allows us to use the previous and the next 
vertex as indices of the vertex scattering matrix instead of the corresponding 
directed bonds, 

for all j,j' in the neighbourhood of i. We shall use both notations at con- 
venience. 

Consider as an example a Neumann vertex with v = 2 and Sf' = {a, /?}. 
It follows from (13. 6p that aa,a = '^/3/3 = 0; t)ut 17^^ = cr/3,a = 1- In other 
words, a V = 2 Neumann vertex transmits waves without any reflection. 
This is consistent with a previous remark that such vertices have no effect 
on the quantum dynamics on graphs, and can be added or removed at will. 

The vertex scattering matrices are the building blocks of much of the 
subsequent theory, and they deserve some further discussion. The vertex 
scattering matrices presented above are but one of many possible examples, 
and they can be constructed to model various physical systems. For example, 
an engineer considers the graph as a model of a network of wave-guides 
connected by junctions. In the ideal situation, there is no dissipation at 
the junctions, and therefore the transmission and reflection at the junctions 
are characterised by unitary scattering matrices. Their details (such as 
their dependence on the frequency), reflect the particular way by which the 
junction is designed and constructed. Also, the quantisation procedure of 
Chapter [5] can be generalised in several ways. We shall dedicate the special 
Section 13.31 to present and discuss the most commonly used classes of vertex 
matrices beyond the present quantisation scheme. We place this section at 
the end of the chapter not to interrupt the flow of the exposition. 



3.2 The quantum evolution map 

The graph is a network of connected vertices, and the waves are scattered 
between them along the connecting bonds. When the graph is assembled 
from all its vertices and bonds, it is important to remember that given two 

connected vertices i and j, a directed bond a = is outgoing from j, 

(i) (i) 
a £ but it is incoming to i, a G S_ . That is, when Cij = 1 

[i,j] = n [j, ^] = n 5? . (3.9) 

The wave function (j3.ip on the bond b = {i, j) was written down by adopting 
the convention that x;, = at the vertex i - this choice emphasises the vertex 
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i as the scattering centre. However, we can write the same wave function in 
terms of the coordinate x^^ = Lb — Xh, so that = marks the vertex j, 

Again, the wave function is expressed in terms of counter propagating waves 
with incoming and outgoing amphtudes with respect to vertex j. However, 
the magnetic flux appears now with a different sign, as befits its role in the 
theory as the element which breaks time reversal invariance. It is therefore 
only natural to consider the magnetic flux as a quantity which is associated 
with a directed bond, with = —A^, in contrast with the bond length 
which does not depend on orientation. Comparing the two expressions (13. ip 
and (j3.1U|) for the wavefunction on the bond b one can read off that the 
outgoing amplitude at the starting vertex of a directed bond a = {b, uj) and 
the incoming amplitude at the next vertex differ by a phase factor 

which the wave aquires from one end to the other of the bond. 

Equipped with all the above notations and remarks, we can now de- 
mand that the bond wave functions ()3.ip or (jS.lOp satisfy the continuity 
and the boundary conditions (|2.6p at all the vertices. This results in 2B ho- 
mogeneous linear equations for the 2B coefficients which can be written 
explicitly as 

V a : aL" = Y,^B{k)a,i3a}^ , (3.12) 

/3 

where the sum over /? is over all 2B directed bonds and the matrix UB{k) 
is written as 

Ueik) = T{k)S{k) , (3.13) 

with T(/c), the bond propagation matrix, which is the diagonal 2B x 2B 
matrix 

T{k)^b,uUb',.') = <566'5..'e'('=+"^'')^'' . (3.14) 
The 2B X 2i?-matrix S{k) contains the vertex scattering coefficients 

c - j^a^ a if follows a at vertex i, a' G J^ia, 

Ja',a[k) ~ 1 Q ' ^jg^ [6. Lb) 

We will refer to S as the graph scattering matrix. A non trivial solution 
of (|3.12p exists for the wave number A; > for which 

(sik) = det(l - Usik)) = . (3.16) 



Csik) is another secular function, whose zeros define the spectrum of the 
graph. In contrast with Chik), it involves a determinant of a matrix of larger 
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dimension, yet it has many advantages, the most outstanding is that Cb(^) 
has no poles on the real k axis, in contrast to Chik) which has poles at the 
Dirichlet spectrum of the graph. The Cs-function can be thought of as the 
characteristic polynomial det (Al — UB{k)) = Yl'j=o^j^''- ^ consequence of 
the unitarity of Usik) is 



This identity is the analogue of the Riemann-Siegel 'look-alike' symmetry 
which holds only approximately for spectral ^-functions of chaotic quan- 
tum systems. Other properties of Cb(^) will be discussed as the theory is 
developed. 

The matrix Z//s(^) will be referred to as the quantum evolution map 
because of the following reasons. Being an involution of two unitary matri- 
ces, it is unitary - which is a basic requirement for a quantum map. The 
action of the map is a composition of two successive operations: scattering 
followed by a propagation along the bonds: The scattering map operates on 
incoming amplitudes at all the vertices and produces the corresponding out- 
going amplitudes. For finite Aj > the graph scattering matrix depends on 
the wavenumber - for Aj = and for some of the generalisations to discuss 
later it becomes independent. The matrix T(k) propagates the outgoing 
waves along the bonds and provides the correct phase for the next scatter- 
ing event. Starting with an arbitrary distribution of amplitudes a™ of waves 
with a wave number k, one can examine the evolution of the wave pattern 
as a function of the number n of scattering events by applying UB^{k) to 
the initial distribution. The consistency condition ()3.12p can be interpreted 
as a requirement that at an eigenvalue of the graph, the wave function is 
stationary with respect to the quantum evolution map - a very natural re- 
quirement. More care has to be taken for vanishing wave number k = 
where e^''^^'' — > 1 is not propagating. That is, only the strictly positive part 
of the spectrum A:„ > is equivalent to positive zeros of the secular func- 
tion while this one-to-one correspondence may fail at /c = 0. The simplest 
example for such a failure is the line of length L with Dirichlet boundary 
conditions, for which Csik) = 1 — e'"^^^. Obviously C-b(O) = in this case, 
but = is not in the spectrum. Note, that for Neumann boundary con- 
ditions C-b(^) is not changed but now /c = is in the spectrum. Such "false 
zeros" of the secular function at A; = appear quite often. For instance, 
quantum graphs with Neumann boundary conditions have a non-degenerate 
ground state at A; = which is just the constant function on the graph 
and the secular function vanishes Csik = 0) = 0. In most topologies this 
is however a highly degenerate zero of the secular function (for instance in 




(3.17) 



which implies 




(3.18) 
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all cases where all valencies are larger than two). In spite of this possible 
failure of the secular function at A; = the graph scattering matrix S com- 
pletely determines the boundary conditions at all vertices of the graph and 
the complete evolution operator UB{k) determines the complete spectrum if 
the limit A; — > is considered correctlj0. 

So far we have shown that the condition (j3.16|) is equivalent to the bound- 
ary conditions p.6p . However, with the above interpretation, we can obtain 
the spectrum of a wider class of graphs, whose vertex scattering matrices 
are unitary matrices derived or postulated in other ways, as discussed in the 
next section. 



3.3 Examples of generalised vertex scattering ma- 
trices 

In the quantisation procedure of Chapter [2] we have considered only bound- 
ary conditions for which the wave function is continuous through each ver- 
tex and we have given the most general account of them. The continuity 
condition may however be relaxed, and more general boundary conditions 
corresponding to a self-adjoint operator on a metric graph have been de- 
rived for this case by Kostrykin and Schrader |152| . Generalised boundary 
conditions at the vertex i can be written in the form 

Ecj']^(..)(4.)) + ^IJd^<^(..)(4.)) = (3-19) 

where the sum runs over all vertices j in the neighbourhood of i. The {k 
independent) Vi x Vi matrices C*-*) and D^^^ can be chosen arbitrarily for 
each vertex, restricted only by the conditions that (i) the Vi x 2vi matrix 
(C^*-*, Z)^*^) has maximal rank and (ii) C^^^D^^^ is hermitian. It then follows 
that the Schrodinger operator on the graph with these boundary conditions 
is self-adjoint. A small exercise gives the corresponding vertex scattering 
matrices. 

a«(/fc) = (i/tl^W + (:7«)-i(ifcD» - C«) . (3.20) 

The vertex scattering matrices defined in this way are unitary and generally 
fc-dependent. With the given restrictions on C^*^ and -D^*-* the matrices 
ifc-D^*-* zbC^*-* are always invertible |152j . However the matrices Z?*-*^ and C^*-* 
are in general not invertible and the limits A; ^ and /c ^ cxd in (|3.2U|) are 
non-trivial. 

^This limit has to take into account, that for k = the general solution of the one- 
dimensional Schrodinger equation (without magnetic field) is a linear function a + bx. 
The rest is straight forward. In the following we will always imply a full knowledge of the 
quantum map when referring to the condition p.l6p such that the oddities at fc = can 
always be removed by the correct limit. We will only come back to this point if it has 
some non-trivial consequences. 
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As an example, let us show that the boundary conditions from Chapter 
E] are a subset of this generalised approach. Choosing the Vi x Vi matrices 

C« = -Ail L»«=E (3.21) 

where Ejj/ = 1 for each matrix element is equivalent to the boundary con- 
ditions ()2.6p and satisfies the conditions on C*-*^ and D^'^\ The vertex scat- 
tering matrix cr^*) = — 1 + ^.^^_^|^, E given by (j3.6p is easily seen to satisfy 
(ifcL>(*) + C(*))(t(*) = i/cZ)(*) - C(*) using the property = u^E. In the limit 
Aj ^ our choice of C^*-* and D*-*^ does not satisfy the required condition 
that (C'-^^D^*^) has maximal rank. However, the choice of C^*^ and D^^^ 
is not unique and it is easy to find a choice that works for Aj = 0. To be 
explicit, for Vi = 3 choose 



C7» = 1 -1 D« = . (3.22) 





Here C'^'-'-D^*-'^ = which is self-adjoint and it is obvious how to extend it 
to higher valencies Vi > 3. 

Instead of defining a self-adjoint Schrodinger operator one may as well 
quantise a metric graph by requiring that a unitary quantum evolution oper- 
ator of the formZ^B(/c) = T{k)S{k) exists. This is equivalent to assigning an 
arbitrary unitary vertex scattering matrix a^^\k) to each vertex i. Such an 
arbitrary choice is in general not equivalent to a self-adjoint Schrodinger op- 
erator on a metric graph with boundary conditions specified by (j3.19p . With 
little abuse of language physicists nonetheless speak of boundary conditions 
which are defined by such arbitrary vertex scattering matrices - and so will 
we in the sequel. While mathematically less satisfactory this is a physically 
well motivated and well defined generalisation. Physically, a vertex is some 
scattering centre which may have some internal structure. While the details 
of this internal structure are not relevant one just describes it through its 
unitary scattering matrix which ensures probability conservation (unitarity 
of the quantum evolution map). The scattering at the vertex may or may 
not depend on the wave-number k depending on the internal structure of the 
vertex. In most physical applications they are just chosen to be constant. 
In the remaining chapters we will stick to an arbitrary but /^-independent 
choice just for simplicity. Note, that a ^-independent choice is also mathe- 
matically satisfying. It has been shown by Carlson [64j that such a choice 
defines a self-adjoint Schrodinger operator on a (directed) metric graph. An 
example of a widely used vertex scattering matrix is the Discrete Fourier 
Transform (DFT) matrix for which 

4'°^^ = ^f'"-^ (3.23) 
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where n(j) is a one-to-one mapping of the Vi neighbouring vertices j to the 
numbers 0, 1, . . . , — 1. In contrast to a vertex with Neumann boundary- 
conditions where backscattering is favoured for large valency v > 4 the 
scattering amplitudes of a vertex with a DFT vertex scattering matrix has 
equal absolute value for all incoming and outgoing bonds. 

In other applications one may have a specific scattering system in mind 
(e.g. the bonds that are connected to a vertex are really channels that couple 
to a quantum billiard). Quite general vertex scattering matrices may also 
arise starting from a quantum graph with Neumann boundary conditions. In 
that case one may define a composite vertex which combines some subset 
of the vertices of the graph and combines them to a single one. This may 
simplify the topology (or calculations) considerably. For the new composite 
vertex one may derive a vertex scattering matrix by eliminating all bonds 
that now belong to the internal structure. As a simple example let us look 
at a composite vertex defined by replacing 



where one bond of length a is eliminated. For Neumann boundary conditions 
a^^ = 1 and cr^') = |E3 — I3 one gets the fc-dependent composite vertex 
scattering matrix 

^comp 1 / e'^"'^ - 1 2(e'2»'= + 1) 
3 + Q>2ak \^2(e'2«^ + 1) e'^"'^ - 1 

Here, we have written the last line in such a way that an interpretation in 
terms of trajectories is apparent. The first two terms — 12 + |lE2 describe 
direct processes for which the particle does not visit the eliminated bond. 
In the third term the two factors 2/3 are the amplitudes for entering and 
exiting the eliminated bond and the factor 

„\2ak °° / 1 \ " 

e _ J2afe^/_igi2afc\ (3 35) 



1 + ie'2»'= n V - / 

o n=0 

is a geometric sum over trajectories for which the particle is scattered back 
to the eliminated bond n times. 

Eliminating bonds by defining composite vertices is very useful in nu- 
merical calculations (it reduces the matrix dimension) , and also in analytical 
approaches. For instance, if there is a length scale separation between small 
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bond lengths inside the composite vertex and large bond lengths in the rest 
of the graph like in 




the vertex scattering matrix of the composite vertex will weakly depend on 
the wavenumber k. 



Chapter 4 

Classical evolution on graphs 



So far we discussed graph dynamics from a quantum mechanical point of 
view. At the present stage, we would like to study graphs from a different 
point of view, which provides the classical counterpart of the quantum theory 
|160t[26lll94j . Usually the connection between the quantum and the classical 
description is provided by quantising a classical system. Here we take the 
process in the reverse direction, for reasons which will be explained below. 

4.1 Classical phase-space and transition probabil- 
ities 

Given a graph, we consider a classical particle which moves freely as long 
as it is on any of the bonds. The vertices are singular points, and it is not 
possible to write down the analogue of the Newton equations at the vertices. 
To circumvent this intrinsic difficulty we employ a Liouvillian approach, 
based on the study of the evolution of phase-space densities. The phase 
space evolution operator assigns transition probabilities between phase space 
points, and it is the classical analogue of the quantum evolution operator. 
We shall employ this analogy to construct the classical dynamics. For this 
purpose we should first establish what is the graph classical phase-space, and 
second, construct the classical analogue of the quantum evolution operator. 

The phase-space description will be constructed on a Poincare section 
which is defined in the following way. Crossing of the section is registered as 
the particle encounters a vertex, thus the "coordinate" on the section is the 
vertex label. The corresponding "momentum" is the direction at which the 
particle emerges from the vertex. This is completely specified by the label 
of the next vertex to be encountered. In other words. 



The set of all possible vertices and directions is equivalent to the set of 
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2B directed bonds. Thus, the classical phase space densities are defined on 
the same space as the corresponding quantum evolution map. 

The evolution on the Poincare section is described in discrete time steps 
(topological time) and is well defined once we postulate the transition prob- 
abilities P^^) > between the directed bonds a and a' where a' follows a 

a < — a — 

at the vertex i, a' G J^ia. Probability conservation then requires 
for all a. 

For a general classical dynamics (Markov process) on a graph one may 
postulate any transition probabilities Pj^}^^ > that satisfies ()4.2p . To con- 
struct the classical analogue of the quantum graphs, we choose the classical 
transition probabilities to be equal to the quantum transition probabilities, 
expressed as the absolute squares of the U b matrix elements 

(4.3) 

In addition to the general condition ()4.2p the transition probabilities on a 
quantum graph also satisfy 

for all a' because the u^*^ are unitary. 

As examples, we quote explicitly the transition probability which corre- 
spond to a DFT, Neumann and Dirichlet vertex scattering matrices. 

— for a DFT vertex scattering matrix, 

^1 — Sa,a' + ^ for Neumann boundary conditions, 
6a a' for Dirichlet boundary conditions. 

(4.5) 

(i) 

In the DFT case a particle on a directed bond a £ Si_ is scattered 

(i) 

with equal probability to all outgoing bond a G Sy . In contrast, for Neu- 
mann boundary conditions at a vertex of high valency w j > 4 backscattering 
a' = a is strongly favoured. One may thus expect that the classical decay 
of correlations is much faster in a graph with DFT vertex scattering ma- 
trices compared to a Neumann graph. The transition probability for the 
Dirichlet case, admits the following physical interpretation. The particle is 
confined to the bond where it started and thus the phase space is parti- 
tioned to non-overlapping components ("tori" - since the dynamics on each 



(4.2) 



p(0 



a 



(4.4) 
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bond is periodic). Thus, the Dirichlet boundary conditions correspond to 
an integrable classical analogue. 

In general, a graph is dynamically connected if the directed bonds 
cannot be split into two nonempty subsets such that all transition proba- 
bilities to go from one set to the other vanish. Dynamical connectivity for 
quantum graphs may equivalently be defined via the scattering amplitudes 



4.2 The classical evolution operator 

(i) 

The transition probabilities P^/^^ can be combined into a 2B x 2B matrix in 
the same way that the vertex scattering matrices have been used to construct 
the graph scattering matrix. This defines the classical evolution operator 
(also called Frobenius- Perron operator) 

-M«',« = \S^'A' = pB.'Ak)\' = llT"^'"'' (4.6) 

10 else. 

M does not involve any metric information on the graph, and for Dirichlet 
or Neumann boundary conditions it is independent of k. The evolution is a 
discrete Markov process with transition probabilities Ma', a- The unitarity 
of the quantum evolution map Usik) guarantees that is a bistochastic 
matrix (also known as doubly stochastic matrix) which means 

'^Ma',a = '^Ma',a = '^ and < A^^/^a < 1 . (4.7) 

a a' 

Bistochastic matrices that are defined by a unitary matrix in the form (14. 6p 
are also called unistochistic matriceqj. These properties of the classical 
evolution map are equivalent to the condition ()4.2p and ()4.4p which ensure 
probability conservation. 

If Pain) > denotes the probability to occupy the directed bond a at the 
(topological) time n, then we can write down a Markovian Master equation 
for the classical density: 

p{n + l) =Mp{n) , (4.8) 

with the 2i?-dimensional vector p = (pi, • • • , P2b)- The total probability to 
be anywhere on the graph is unity thus 

llp|| = 5^Pa = l- (4.9) 

a 

^In general unistochastic matrices are a subset of bistochastic matrices. A unistochastic 
matrix defines a Markov process which can be quantised (though not uniquely) [351 12581 

BO]. 
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From the bistochastic property of M follows that the uniformly dis- 
tributed density on the directed bonds 

pr = ^ (4.10) 

is invariant under the classical evolution 

pi°^ = Mp'""^ . (4.11) 

In most cases there exists only one invariant probability distribution. It 
then follows that the Markov process is ergodic which means that the time 
average of any classical observable f £ (that is, fa is the value of the 
observable on the directed bond a) is equal to an average over the invariant 
uniform probability distribution 



N 

lim V / • 7W"p(0) = / • p'"^^ = y ^ (4.12) 



?i=0 a 



for every initial distribution p(0). This is equivalent to the statement that 
the time averaged occupation probability on a directed bond a is uniformly 
distributed over a 



1 ^ 1 



AT^OO + _ 

n=0 

This is not a very strong statement on the classical dynamics of a graph. It 
is known that the classical dynamics for every dynamically connected graph 
is ergodic [26j . A Markov process may have the stronger dynamical property 
of being mixing which is defined by 

lim 7W"p(0) = p'""^ (4.14) 

n^oo 

for any initial probability distribution p(0). 

The properties of the graph which determine whether it is ergodic or 
mixing are encoded in the classical spectrum , £=!,•••, 2B which is 
defined by 

Mxi = viXi ■ (4.15) 

The classical spectrum is restricted to the interior of the unit circle jf^l < 
1 and v\ = \ corresponds to the uniform distribution Xx = p™^. Any 
probability distribution can be written in the form 

2B 



J^a^X^ (4.16) 
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and ai = 1. For f ^ / 1 one has = J2a Xa,e = due to the probability 

conservation of the dynamics. 

Ergodicity implies that z^i = 1 is the only unit eigenvector. Thus there 
is a spectral gap 

Ao = min \1 - uA = ^— (4.17) 
^ e=2,...,2B ' ' n'^'^s ^ ^ 

which determines the time scale n^^^ on which the left hand side of (|4.13|) 
decays to the uniform distribution. Thus, ergodicity allows for eigenvalues 
on the unit circle vi = e'^ which are damped by the time average. In practice, 
such eigenvalues are non-generic - however there are some example, such as 
star graphs which have an eigenvalue = — lH. 

If there is no eigenvalue on the unit circle apart from vi = 1 all non- 
uniform modes of the probability distribution decay 

2B 

M^p{^) = p''^^ + u^a.xi P'"" (4.18) 

1=2 

and the dynamics is mixing. The rate at which equilibrium is approached 
is determined by the gap 

Kg = min 1 - = (4.19) 

^ e=2,...,2B n™"" 

between the next largest eigenvalue and 1. In general ergodicity sets in 
faster than mixing n'^^^ < n™'^. 

We shall end this section by showing how the trajectories which we 
defined formally in Section Fl . 1 1 acquire a dynamical significance, and emerge 
naturally in the classical framework developed above. 

The classical probability to make a transition from a directed bond a to 
a directed bond (3 after n steps equals 

ter„(/3,o) 

where the sum extends over the set a) of trajectories which start at a 

and get to (3 in n steps, and 

n-l 



^Aa(t)=n-^°=+i."^ (4-21) 



s=l 



^Star graphs are an example of bipartite graphs for which the directed bonds can 
be spht into two nonempty sets and the transition probabihties within each of the two 
sets vanish. For star graphs due to its topology there is no transition from any outgoing 
(incoming) bond to another outgoing (incoming) bond. One may easily construct other 
examples which rely on the dynamical connectivity as well. Dynamically connected star 
graphs are generally mixing if the classical dynamics is defined on bonds instead of directed 
bonds as in Chapters |8] and [O] 
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are the classical weights or probabilities contributed by the trajectory t to 
the transition probability. This form of the transition probability emphasises 
the similarity between the classical dynamics on a graph and the evolution 
under Hamiltonian maps. 

One of the characteristic features of chaotic (mixing) classical dynamics 
is the exponential proliferation of the number of classical trajectories which 
connect the same initial and final points as the transition time increases. 
This property is also shared by the trajectories on the graphs. This follows 
immediately from the fact that we have a natural code - the sequence of 
directed bonds - which associates a string of numbers with a trajectory 
t = (ai, . . . ,a„) (by numbering the directed bonds from 1 to 2B), and a 
connectivity matrix which established a Markovian grammar on the codes. 
In the language used in the theory of dynamical systems this is a Bernoulli 
code, which in Hamiltonian systems guarantees chaotic dynamics (see \120\ 
[79] for more details). On a simple graph, the sequence of vertices is often 
used to define an alternative (equivalent) Bernoulli code. 

4.3 The return probability and a classical sum rule 

Of prime importance in the discussion of the relation between the classical 
and the quantum dynamics are the traces 

2B 

= tr>l" = ^ (>!"),,, . (4.22) 

a=l 

This is the return probability - the classical probability to perform n- 
periodic motion. The mean probability to return is Un/2B. For mixing 
dynamics where only one eigenvalue v is on the unit circle, one has 

2B 

Un = Y, 1 • (4-23) 

Following the discussion at the end of the previous section, Un can be 
written also as a classical sum over closed trajectories of length n. The 
classical weight of a closed trajectory tc with a code (ai, ai) is the 

product 

n 

Vt^ = llM^^^,,a. (4.24) 

s=l 

where a„+i = ai. The weight can be interpreted as the probability to 
remain on the trajectory. Any cyclic permutation of the bond indices in the 
code of a closed trajectory is also a closed trajectory with identical classical 
weight. Hence, instead of summing over closed trajectories, we may express 
the return probability as a sum over periodic orbits . A code of an n-periodic 
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orbit can be written as a repetition of primitive codes whose lengths Up 
divides n. The set of all the primitive orbits which build n periodic orbits 
will be denoted by V{n). Denote the weight of a primitive orbit p by Wp. 
Then, each primitive orbit contributes to the return probability a term which 
consists of its weight taken to its r power, multiplied by Up = n/r, 



This is a basic formula in all our subsequent discussion and it deserves a few 
comments. 

The weights Wp can also be defined for codes which violate the dynamical 
connectivity, by assigning them the value Wp = 0. Sums over periodic orbits 
that involve the weight as a factor can immediately be replaced by sums over 
all arbitrary cyclic codes for this reason. 

The code gives rise to a variety of classifications of periodic orbits. For 
instance, a periodic orbit is either reducible or irreducible. An irreducible 
periodic orbits visits each directed bond at most once. All directed bonds 
appearing in its code are different ai ^ aj. The maximal period of an irre- 
ducible periodic orbit is rimax = 2i?. All other periodic orbits are reducible, 
that is at least one directed bond is visited more than once. All periodic 
orbits of period larger 2B are reducibl^. Reducible orbits have the remark- 
able property that their weights can be written as products of weight of 
shorter orbits Wp = Wp^Wp^. 

We are now equipped with the background to obtain a classical sum- 
rule for mixing dynamics by substituting a periodic-orbit expansion of 



^The definition of reducible and irreducible periodic orbits depends on the code 
wliicli is used to define a periodic orbit unambiguously. For a simple graph one may 
use the the sequence of vertices p = iii2...in as a code (symbolic dynamics) which 
uniquely defines a periodic orbit. In the directed bond code the same orbit is given by 
p ~ [i2, ii], [is, 12], ■ • • , [ii,in]- Irreducible periodic orbits with respect to the vertex code 
{vertex-irreducible periodic orbits) do not intersect at any vertex such that all ij are 
different. Our definition which is based on the directed bond code is not equivalent to 
the one based on the vertex code. In general there are bond-irreducible orbits of period 
n > V and the latter cannot be vertex-irreducible. Both codes have their advantages and 
disadvantages. The vertex code has the advantage of a smaller number of symbols and 
a grammar which is readily encoded in the connectivity matrix. In the present context 
the grammar is less relevant since a code which does not correspond to a periodic orbit 
has zero weight. The advantage of using the bond code to define reducible orbits, is that 
weights of reducible orbits are products of weights of non-irreducible orbits. 




(4.25) 



p6P{n) 
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Combining (|i:23]) and (f05|) we ge10 

un= ^WpY^^ 1- (4-27) 

This sum rule is analogous to the Sinai-Bowen-Ruelle \2'27\ [BOl 1210] 

sum rule for classically chaotic dynamics (in the physics literature also 
known as Hannay-Ozorio de Almeida sum rule 1137) 1. The weights Wp 
are the counterparts of the stability weights | det(/ — Mp)\~^ for hyperbolic 
periodic orbits in Hamiltonian systems, where Mp is the monodromy ma- 
trix. Graphs, however, are one dimensional and the motion on the bonds 
is simple and stable. Ergodic or mixing dynamics is generated because at 
each vertex a (Markovian) choice of one out of v directions is made. Thus, 
chaos on graphs originates from the multiple connectivity of the (otherwise 
linear) system. 

Using the expression ()4.27p for u„ one can easily write down the complete 
thermodynamic formalism for the graph. Here, we shall only quote the 
periodic orbit expression for the Ruelle ^-function 

Cr{z) = det{l - zM)-^ =e-trin(i-^M) ^gE„^«n 

^TT 1 (4.28) 

p (1 - z"'" exp(-np7p)) 

where the product extends over all primitive periodic orbits p and we have 
set Wp = . The Ruelle (^-function is an important tool to analyse 

dynamical systems |1201 [79] . 



* In the non-mixing case a dynamically connected graph is ergodic. With an additional 
time average over an interval An ^ n one then gets the sum rule 

n + An— 1 

hT. -^-i^ E -.^^ 1 • (4-26) 



Chapter 5 

Spectral theory for quantum 
graphs 

Trace formulae are an important tool in spectral theory. They express spec- 
tral sums in terms of sums over periodic orbits of the underlying classical 
dynamics. In this chapter our main goal is to introduce the trace formula 
for the density of states of a quantum graph. The first trace formula for a 
quantum graph was derived by Roth [2061 1207j who showed that a spectral 
determinant (which can be interpreted as a variant of the secular function 
C-b(^)) can be written as a sum over periodic orbits. The trace formula for 
the density of states goes back to Kottos and Smilansky [159^ I160j and is for- 
mulated in analogy to the semi-classical Gutzwiller trace formula [l'S'2\ I133| 
for chaotic (hyperbolic) Hamiltonian quantum systems. In contrast to the 
latter, the trace formula for quantum graphs is exact. Apart from the den- 
sity of states we will also discuss how other spectral functions are related to 
sums over periodic orbits. 

The starting point for the derivation is the secular function Csik) which 
is expressed in terms of the quantum evolution map Usik) (see Chapter 
[3]). Thus, the boundary conditions at the vertices are defined by a given 
set of vertex scattering matrices and the wave functions are not necessarily 
continuous across a vertex. Without loss of generality we will restrict the 
discussion to dynamically connected quantum graphs of simple topology. 
If the graph is not dynamically connected (e.g. a graph with Dirichlet 
boundary conditions) the graph can be divided into subgraphs which have 
independent spectra and wave functions - that is, the complete spectrum 
is the superposition of the independent spectra. For simplicity, we will also 
assume throughout this chapter that the vertex scattering matrices do not 
depend on the wave number k (the generalisation to k dependent scattering 
matrices follows in a straight forward way). 
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5.1 The density of states and the counting func- 
tion 

If {kn} (n G N) is the spectrum of a quantum graph (degenerate eigenvalues 
appear according to their multiphcity) we define the density of states as 

oo 

d{k) = ^6{k-kn) . (5.1) 

n=l 

The spectral counting function which provides the number of eigenvalues 
kn which are smaller than k is given by the integral 

N{k) = 9{k)No + lim / dk'd{k') = V ^(yfc - kn) , (5.2) 

n=l 



where A'^o is the number of eigenvalues at A; = and 

0{k)-- 



for < 

1/2 forA: = (5.3) 



^1 for /c > 

is the Heaviside step function. The limit 5 — > 0"*" is needed only when N{k) 
is computed at /c ^ 0, and it ensures that limfc^o-^(^) = -^o- When the 
secular function Cb(^) does not have a zero at A; = and Nq = one can 
set 5 = from the start in the following discussion. 

We will now express the counting functions N{k) in terms of the secular 
function C-b(^) = det(l — UB{k)). To this end, we introduce the eigen- 
values {e'^t('')} (£ = 1,2,..., 2B) of the quantum evolution map Usik) at 
fixed wave number k and first discuss some of their properties. The phases 
(j)i{k) are important in the present context since the quantisation condition 
C-b(^) = is equivalent to the requirement 

(l)i{k) = 2'Kz with z = 0,1,2,... (5.4) 

for one of the (peik)- Each of the quantisation conditions (15. 4|) yields a 
discrete subset {kn^} C {kn} of the complete spectrum and this subset is 
free of degeneracies since ^ > 0. Indeed, a direct computation reveals that 

^ = Ei.(i«!t,wi^ + K?-„,wr) > «. (5.5) 

b=l 

it) 

where o}{p^^ the components of the I eigenvector. This proof of the mono- 
tonic increase of the eigenphases depends crucially on the assumption that 
the vertex scattering matrices are independent of k. It follows immediately 
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that the maximal degeneracy in a quantum graph is 2B. Direct inspec- 
tion shows Ljnin < ^ < -^max where Lmin (-^max) IS the Smallest (largest) 
bond length, and Yl'e=i ^ = 2BL where L = ^ J2b=i mean bond 

length. For Neumann graphs with rationally independent bond lengths a 
much stronger statement has been proven: their spectra are generically de- 
void of degeneracies [116j . 

We also note that Usik + ie) is subunitarity for any finite e > 0. 
That is, I e''^'^ I < 1. This is an immediate consequence of the obvious 
subunitarity of the bond propagator T{k + ie). 

A counting function for the subspectrum {kn^} is obtained through the 
correspondence to the zeros of the real function sin^^^. A sum over step 
functions at the (strictly positive part of the) subspectrum can thus be 
written as 

Y: 0{k - feW) = Im An sin _ in sin (5.6) 



where the limit e — > 0"*" is always implied. The e-shift ensures that the 

2 



complex number sin '^"(^+'^) rotates counterclockwise around the origin of 



the complex plane. The term In sin ^i^^ ensures that one starts to count at 
zero and e has been set to zero to indicate that the limit e ^ 0"*" has to be 
performed before the limit 5^0^. 

The full counting function is obtained by summing over all 2B subspec- 



tra. Using sin = 2 (1 _ e '2 ) we get, 



2 2 



Nik) =No + ^Imlndet^^i^ - llmln^^^^ + 



2vr_ Ub{S) vr Cb{S) 

BL 1 1 

=iVo H k+ -ImlnCB((5) ImlnCB(fc+ ie) • 

vr vr vr 

We have used the assumption that the vertex scattering matrices do not 
depend on k in the second line. 

The counting function can be decomposed into a smooth term (the so- 
called Weyl term) and an oscillatory term 

N{k) = A^Weyl(^) + Ar°«=(A;) . (5.8) 
The smooth term is given by 

vr 

where the first term —k describes the linear increase of the counting func- 
tion and 

iVWeyi(o) = iVo + ilmlnCij(5). (5.10) 
vr 
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is a constant that depends on the boundary conditions at the vertices. This 
constant ensures that the complete counting reduces to A'q at A; = 6. When 
one wants to calculate N^^^^{0) from (jS.lOp one should be aware that the 
logarithm is a multivalued function and it is not always obvious which is 
the correct sheet. The correct sheet can be identified by the condition 
^J^ dk (A^Weyi(^) _ ^ which just means that the difference is 

oscillating around zero with a bounded amplitude. For general Z^b(O) the 
shift N^^^^{0) can take real values. In many cases (e.g. for Neumann 
boundary conditions) Ub{0) = 5" is a real unitary matrix. Then, A^^^^'(O) 
is some half-integer number. This can be calculated directly by observing 
that all eigenvalues of S with non- vanishing imaginary part appear in com- 
plex conjugated pairs such that their contributions cancel each other. Real 
eigenvalues take the values ±1. —1 does not contribute since the imaginary 
part of In (l — (— e'^)) vanishes for 6^0. Only eigenvalues +1 contribute 
ilmln (1 - e'-^) ~ ^Imln e-'"'/'^6 = and iVWeyi^o) = iVo - f where z is 
the number of unit eigenvalues. 

The linear dependence of N'^'^y^{k) on /c is a consequence of the fact 
that the graph is a one dimensional object. The fact that it is not simply 
connected appears only in the oscillatory part of the counting function - the 
second term in (15.80 - 

iV°s=(A:) = --ImlnCB(A; + ie) . (5.11) 
vr 

The oscillations are bounded by the number of bonds \N°^^{k + ie)| < B. 
Equation (j5.11|) will be the starting point for our derivation of the trace 
formula in Section [5.21 

It is now straight forward to derive the density of states by differentiation 
with respect to A: as a sum 

d{k) = d'^^y^k) + d°'%k) (5.12) 

of a smooth Weyl term 

^Weyl(^) = EL (5.13) 

vr 

and an oscillatory part 

d°'%k) = ---^ImlnCB(fc + ie) . (5.14) 
vr dA; 

5.1.1 The spectrum of the quantum evolution map 

The 2B eigenvalues e"'^^'-'^-' of the quantum evolution map lAsik) at fixed k 
can also be used to define a density and a corresponding counting function 
which will turn out to be interesting objects in their own right. We define 
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the density of eigenphases as 

2B 



d{ct>;k) = Y,52A(l>-Mk)) (5.15) 



and the eigenphase counting function as 

N oo 

64>'d{4>'; k) = 

Is 



„0 N oo 

NicP; k) = No e{(l)) + / d,/.'(iV; A:) = ^ ^((/> - </'n(A:) + 2t:u) . (5.16) 

J8 1 n 



n=l v=0 

In (j5.15p 527r(0) = YlV=-oo ^{4> + ^iriy) denotes the 27r-periodic (5-function. 

In complete analogy to the derivation for the spectral counting function 
above one may write 

N{(j); k) = iV^^^y^^; k) + iV°^'=(</.; it) (5.17) 

and express the smooth and oscillatory parts in terms of the (/)-dependent 
secular function 

Cb(</'; k) = det (1 - e-'^Usik)) . (5.18) 
As a result one obtains 

iVW<'yi(0; k) = -(/) + No- -ImlnCB((^, k) (5.19) 

vr TT 

for the smooth part, and 

iV°'^=((/,;/fc) = ilmlnCB((/'- ie,fc) (5.20) 

TT 

for the oscillating part. These expressions for the eigenphase counting func- 
tion can alternatively be derived using Poisson's formula. Note that the 
second term of the Weyl part — ^Im In (^b i^-, k) is independent of (j), however, 
as a function of k it can oscillate with a maximal amplitude of B. 
The smooth Weyl part of the density of eigenphases 

d^^y\cj)-k) = - (5.21) 
vr 

is consistent with having 2B eigenvalues distributed in an interval of length 
2vr. Finally, the oscillating part of the density of eigenphases is given by 

1 rl 

d"'%(p;k) = -— ImlnCB((A- ie,A;) . (5.22) 
vr 0(p 
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5.2 Periodic orbits and the trace formula 

We will now consider the oscillating parts of the spectral functions N(k) and 
d{k) from a periodic-orbit perspective. One of the reasons for the interest 
in quantum graphs in the quantum chaos community is the analogy of an 
exact trace formula for quantum graphs to Gutzwiller's semiclassical trace 
formula for the density of states of chaotic Hamiltonian systems. While 
there is no underlying deterministic classical dynamics for graphs they still 
display the generic behaviour of chaotic Hamiltonian systems. At the same 
time they are less resistant to either rigorous or numerical approaches. 

The logarithm of the secular equation appearing in (15. lip can directly 
be written in terms of traces of powers of the quantum evolution map 

oo 

IndeUl - UBik)) =tvln(l- UB(k)) = -tr Kb" (k) . (5.23) 

n 

n=l 

We can now follow the path we used in Section 14.31 to compute the classical 
return probability u„ = tr in terms of periodic orbits. The result is 
again a sum over primitive periodic orbits [16Dj 



oo 

Ar°«^(A;) =-Im V -tvUB^'ik + ie) 



vr ^ — ' n 

n=l 

oo ^ 

=Im Yl — ^;e'''^''(^+'^)+''-*^ (5.24) 

p r=l 
oo ^ 

= Y1 J^— l-4p^sinr(Lp/c + $p + /Up)e-^^'^ 

p r=l 

In the second line we have introduced the sum over all primitive periodic 
orbits p of the graph and their repetitions r. Here, Lp = Yl^Zi is the 
length of the primitive periodic orbit p = ai, . . . , a„p, 

Ap = |.Ap|e ^ = SaiOnpSan^Onp-i • ■ ■ S a2ai (5.25) 

is the quantum amplitude of the primitive orbit defined as the product of 
all scattering amplitudes along the orbit, and 



^p = Y,^c.,=Y.^iA,, (5.26) 

1=1 1=1 

is the overall magnetic flux through the periodic orbit. The step from the 
first line of (j5.24p to the second is performed by collecting all contributions 
to the n-th trace that stem from the same primitive periodic orbit such that 
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n = rup. After a resummation over the repetitions r in the second hne of 
(I5.24P the oscillatory part of the counting function can be expressed as 

N^'^ik) = Im i ^ In (l - ^pe'^''('=+'^)+'*'') . (5.27) 
p 

It should be noted, that these period-orbit sums are not absolutely con- 
vergent due to the exponential proliferation of periodic orbits. That is, 
the number of periodic orbits of period n grows exponentially in n. Abso- 
lute convergence is only obtained for complex k beyond an entropy barrier 
Imk > Scrit- For real k the limit e ^ implies that the infinite sum over 
periodic orbits is ordered with respect to the length rLp of the periodic 
orbit. 

The derivative of (j5.24p with respect to k leads to the, equally exact, 
trace formula 

OO 

d°''ik) =Re Yl ^A;e"-^^'^^+'"^+'""^'' 

p r=l 

OO 

= YY1 —l-^pT COS r{Lpk + % + ap)e~^p' (5.28) 

p r=l 

=Re Y 



TT 1 - ^„e'^!'(^+'^)+'*p 
p ^ 

for the oscillatory part of the density of states. Instead of the sum over 
primitive periodic orbits and its repetitions, one may combine the two sums 
in a single sum over all periodic orbits d°^'^{k) = Re Yip In that 

case Ap is the amplitude of the full orbit, and Lp its full length while the 
primitive length is Lp/rp where is the repetition number of the periodic 
orbit. Note, that the classical probability to stay on a periodic orbit (j4.24p 
is just the absolute square of the quantum amplitude 

Wp = \Ap\'^ . (5.29) 

Comparing the trace formula (j5.28p to the Gutzwiller trace formula reveals 
a complete analogy. 



5.2.1 Trace formulae for the eigenphase spectrum 

The oscillatory parts of the spectral functions A'^(<^; k) and d{(j); k) for the 
eigenphase spectrum of quantum evolution map at fixed wavenumber k can 
be written in terms of primitive periodic orbits in an analogous way. This 
leads to the trace formulae 

OO 

iV°'^^((/,; k) = -Im Y Y — ^;;e''^(^f'=+*f-™p('^-i^)) (5.30) 

p r=l 
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and 



oo 

rir. 



d°'^=(<^; fc) = Re ^ ^ ^_4rgir(Lpfc+$p-np(</,-i.)) _ (5 3^^ 

p r=l 



5.3 The length spectrum and 

the quantum-to-classical duality 

The trace formula (|5.28|) for the density of states may be written in the form 

/oo 
dL A{Ly^^ (5.32) 

where 5^0 will always be implied, and 

bT 1 °° 

A{L) = 5{L) + -YY LpAle"''^'^6{L - rLp) (5.33) 

p r=l 

is a sum over primitive periodic orbits and their repetitions plus a zero length 
contribution (related to the Weyl term) that is proportional to the length 
BL of the graph. The complex function A{L) is a weighted sum over 6- 
functions located at the length spectrum {Lp} (and their repetitions). This 
reveals a quantum-to-classical duality between the length spectrum {Lp} 
and the spectrum of the quantum graph. The relation may also be 
inverted 

1 fOO 1 

A{L) = - dk d{k)e-'^^ = - y e~'^"^ . (5.34) 

The length spectrum is highly degenerate. This follows from the fact 
that the length of any periodic orbit of period n is an integer combination 
of the bond lenths 

B 

Lp = '^qbLb 96 = 0,1,2,..., (5.35) 
6=1 

with X^^i Qb = B. Periodic orbits with the same length differ only in the 
order by which the bonds are traversed. Note, that not every combination 
of integers qb is consistent with the connectivity of the graph. The degen- 
eracy class will be denoted by the vector {q} where q = {qi, . . . , qs)- For 
star graphs the degeneracies are maximal, only even q^, are allowed by the 
connectivity without any further restriction. 

The presence of the exact 'quantum-to-classical' duality imprinted by 
the relation between the length spectrum {Lp} and the quantum spec- 
trum {kn} has been used efficiently in the context of isospectrality. In a 
variant of Kac' famous question 'Can one hear the shape of a drum?' jl45j 
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one may ask the question if there are topolgically or metrically different 
quantum graphs which have the same spectrum {kn} [H]. Such graphs are 
called isospectral. Kac' question can be reformulated as 'Can one hear the 
shape of a graph?'. Exploiting the quantum-to-classical duality an affirma- 
tive answer has been proven in if the class of graphs is restricted to 
simple topologies and rationally independent bond lengths. In contrast, if 
the restriction of rationally independent bond lengths is dropped, a lot of 
examples of isospectral graphs are known today. Still, isospectrality is not 
understood in sufficient detail. Recently, it has been asked [2211 how one 
may resolve isopectrality for quantum graphs by adding non-trivial addi- 
tional information such as nodal counting |124j . 

5.4 The secular function is a (^-function 

Let us set tp = Ape'^^p for a primitive orbit p. For a fixed primitive orbit 
p the sum over repetitions r in the trace formula for the counting function 
(I5.25P can easily be performed X^^i ^p/^ ~ ~^^{^ ~ *p)- Summing over all 
primitive periodic orbits and comparing to (|5.24p then leads to the exact 
expression 



for the secular function, a truly remarkable identity. The expression on the 
left hand side is defined as the determinant of 1 —UBB{k). As such, when 
expanded in periodic orbits it has contributions from irreducible periodic 
orbits only [7]. A superficial look on the right hand side could lead to the 
conclusion that it involves all possible multiples of periodic orbits and their 
repetions. This apparent paradox is resolved by observing that there is an 
exact cancellation mechanism which is responsible for the identity of the 
left and right hand sides. From this follows that a truncation of the infinite 
product or an approximation of any of the terms tp leads to systematic errors 
and therefore (j5.36|) has no practical advantage for graphs. This contrasts 
the situation in the analogous semiclassical trace formulae where shadow- 
ing of periodic orbits is the analogue of the exact cancellations mentioned 
above. Nevertheless, for pedagogical reasons mainly, we shall investigate 
further. 

The similarity to the Euler product formula for the Riemann C-function 



CRiemann (*) = Dprimes: n ( ^ ~ ) justifies the letter C by which the secu- 



lar function is denoted. The product over primitive periodic orbits in the 
secular function corresponds to a product over primes in the Riemann C- 
function. Like the periodic-orbit sum the infinite product on the right hand 
side of ()5.36p is absolutely converging only for complex k beyond the en- 
tropy barrier. For real k the product is defined by analytic continuation. 
Note that a complex k with positive imaginary part, orders the product 




(5.36) 



p 
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by the length of the periodic orbit. Thus this ordering is imphed in the 
analytic continuation. Ordering the contributions by the period n instead 
of the length is essentially equivalent. Other non-equivalent orderings gen- 
erally diverge, or worse: they converge to a wrong result (with false zeros 
|147j ). Note, that generally \tp\ < 1 for a periodic orbit (if we exclude the 
line, that is the trivial graphs with B = 1). Thus the zeros of the secular 
function which define the spectrum of the graph are not equivalent to zeros 
of the factors 1 — tp{k) (which have zeros only for complex k with negative 
imaginary part). Instead the zeros of the secular equation for real k are due 
to the infinite product over primitive periodic orbits. 

Let us now come back to the exact cancellations on the right hand side 
of (j5.36p which eventually reduce it to the finite polynomial on the left hand 
side. One principle behind these cancellations which is known as shadowing 
may be understood from a closer look at the the expression tp for reducible 
primitive orbits. Let pa and pb be two irreducible primitive orbits with codes 
a = aPa and h = aPb which have exactly one directed bond in common. 
Here, Pa and Ph are the codes for two paths that close the periodic orbits 
such that Pa and Pb never visit the same directed bond. For the reducible 
primitive periodic orbit pab with code ab = aPadPb one gets the identity 
tab = tah [7J. Longer reducible primitive orbits may be obtained from pa 
and Pb- With an obvious short-hand notation for the codes, aab and abb 
are the only possibilities to create reducible primitive orbits composed of 
three irreducible orbits (with four irreducible orbits one has aaab, aabb, and 
abbb; for five aaaab, aaabb, aabab, aabbb, ababb, and abbbb). Performing the 
product over the primitive orbits p that are composed of pa and pb one gets 

^(1 - tp) =(1 - ta){l - tb){l - tab){l - taab){l " Ufe) • • • 
P 

= (1 -ta-tb + tatb){l - tatb){l " t^4)(l " tatfo) • • • ^^'^^'^ 
= (1 -ta-tb + tltb + tail + tltl){l - tltb){l - tatl) . . . 

In the third line all contributions of second order in ta and tb have cancelled. 
It is obvious at this stage that the third order terms will cancel after the 
next two factors have been multiplied. Eventually all higher order terms 
cancel and one gets the simple result Wp{^ ~ ip) = ^ — ta — tb which is a 
finite polynomial that only depends on reducible primitive orbits (and its 
zeros are not zeros of one factor 1 — tpj^. If more irreducible orbits are taken 
into account a direct multiplication gets much more involved but the exact 
shadowing of reducible periodic orbits remains valid. Eventually, when all 

^One may prove, that each finite product over factors 1 — tp equals 1 ~ ta — tt + R 
with R = 0{t") if all contributions from periodic orbits that are composed of less than 
n irreducible periodic orbits appear in the finite product. It is a signature of the bad 
convergence of the infinite product that the modulus of the remainder for real k in 
generally does not go to zero in the limit of large n. 
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factors 1 — tp are expanded and the shadowing principle is used to replace 
reducible orbits by products of irreducible orbits the final result only depends 
on the finite set of irreducible orbits and their amplitudes tp. Note, that all 
these orbits have a period smaller than (or equal to) the dimension 2B of 
the quantum mapH. A more systematic treatment of approximations to the 
secular function for calculating the spectrum from a finite set of primitive 
periodic orbits is obtained within the cycle (or curvature) expansion [79l 
I12U) . The cycle expansion also holds for the Ruelle (^-function Cr{z) (|4.28|) 
that describes the classical dynamics on the graph. 

Cycle expansions are also used efficiently for semiclassical expansions of 
(^-functions in Hamiltonian quantum systems [52]. As soon as semiclassical 
approximations are used the analytical properties of the (^-function are de- 
stroyed. Also, replacing long orbits by products of short orbits (also called 
pseudo-orbits in this context) is a semiclassical approximation in these sys- 
tems. For quantum graphs the expansion of the secular function into prod- 
ucts of irreducible periodic orbits is exact due to the finite dimension of the 
quantum evolution map - that is, due to the discreteness of phase space. 



^The final result can also be obtained with Newton's formulas which express the coef- 
ficients a„ of a secular determinant det{z — A) = X]^=o through the first A'^ traces 
trA" in a recursive way. 



Chapter 6 

Spectral statistics 



Any respectable statistical theory starts by specifying the space of variables 
over which the probability distribution (measure) is defined, and where the 
statistical tests and conclusions are drawn. Here we would like to study the 
spectrum of a single graph or of an ensemble of graphs, from a statistical 
point of view. To remain within the range of acceptable respectability we 
should at least state what is our space of variables. In other words, what 
are we averaging over to get a statistical (probabilistic) description. The 
statistical approach to the spectra of graphs can be developed at several 
levels. The conceptually simplest level is to consider the ensemble of all 
the graphs which share the same connectivity but whose bond lengths, and 
possibly also vertex scattering matrices are drawn randomly according to 
prescribed probability distribution functions. This approach is very similar 
to disorder averaging common e.g., in the study of mesoscopic systems. The 
less obvious approach is the study of the spectral statistics of a single graph. 
Here, one can consider non overlapping spectral domains and study the 
distribution of the spectral intervals within each of the domains, assigning 
to the domains equal probabilities. This approach, which is reffered to as 
spectral averaging, makes sense only if one can prove that the probabilities 
converge to a well defined limit distribution as the spectral domains and 
their number increase. 



6.1 Spectral averages 

The spectral average over some interval ko < k < ko + K of a spectral 
function f{k) of a quantum graph is defined by 




(6.1) 
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In the simplest case one chooses f{k) = d{k) and calculates the mean density 
of states 

{d{k))u = d'^'y' + -^(iV°^^(A;o + K)- N-'-{ko)). (6.2) 

Since |A^°^'^(A;)| < B one may neglect the second term \i K =. Thus the 
spectral average of the density of states equals the Weyl term if the spectral 
interval is chosen larger than the inverse mean bond length (equivalently, 
when the number of states in that interval is large compared to the dimen- 
sion of the quantum evolution map). In most cases we will take an infinite 
spectral interval K ^ co for averaging. For finite K the result explicitly 
depends on k^. For small K < l/d^^^^^ = A where A is the mean level 
spacing the 5-peaks are broadened but generally do not overlap. More in- 
teresting is the regime A ^ < = where pronounced oscillations modulate 
the Weyl term. Replacing the periodic orbit sum for A^°'^'^ one sees that the 
most pronounced of theses oscillations ~ e'^p'^" stem from the short periodic 
orbits with length Lp. 

Some spectral averages remain very spiky functions even after the spec- 
tral averaging. A simple example for this is the correlator 



{d{k)d{k + k'))k =025{k- kn)S{k' + kn-kr, 

\m,n I ^ 

<K 

=dW^y'<5(A:') + lim -- y 6{k' + kn-k„ 



(6.3) 



where we have taken the average over the complete spectrum and assumed 
that there are no degeneracies. In the limit K ^ co the second part be- 
comes an infinite sum over (^-functions with a weight ~ \/K which is often 
equivalent (in a distributional sense) to some well-behaved function. Both 
in numerical and in analytical approaches an additional average over a short 
interval in k' is usually needed to make a well-behaved function apparent. 
In chaotic systems the interesting features of such correlators occur on the 
scale of the mean level spacing A = ^^ws^ = t^/BL and the additional aver- 
age over k' has to be on a much shorter scale. Equivalently, one may replace 
the 5- functions by a Lorentzian or Gaussian function of width AA:'. The 
Lorentzian is natural when the trace formula is used. The addition of an 
imaginary part ie to the wavenumber results in replacing the (5-functions in 
the density of states by Lorentzians of width e. 

In practical (numerical) calculations one can circumvent the difficulties 
encountered in direct application of equation ()6.3p by dividing the spec- 
trum of length K into M subsequent intervals. For each of the intervals 
one can use (|6.3p replacing the (5-functions by Lorentzians or Gaussians of 
finite width. The average over all the intervals is a well-behaved numerical 
approximation. As a rule of thumb one should use 
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6.2 Disorder averages 

In the study of spectral statistics averaging over an ensemble of systems, 
usually referred to as 'disorder average', is a usefull theoretical tool which is 
naturally called for in many applications. Disorder averaging for graphs can 
be implemented in several ways. The most obvious one is to take a family of 
graphs with the same connectivity but with a random distribution of bond 
lengths. Another possible variable to average over are the vertex scattering 
matrices which can be picked at random from some well-defined ensemble. 
Both methods were used in past works and we illustrate this method in 
Section 18.51 For graphs there is an intimate connection between disorder 
and spectral averages which will be explained in the next section. 



6.3 Phase ergodicity for incommensurate graphs 

All spectral functions f{k) that we are going to consider depend on k via 
the B phase factors e'^^*" in the bond propagator T{k). They are thus 
quasiperiodic functions of k. Let us again assume that the bond lengths 
are incommensurate. Under this condition the flow 

covers the 5-dimensional torus ergodically. With = fcLf, and 

modest conditions on the function f{k) = f{e^^^^^\...^&^'^^^^'>) we have 
thus the exact equality |144l [25] 

{f{k))u = lim 1 r 6k /(e'-^if^), . . . , e'^^^'^)) 

^(/(e'<^\...,e'<^«))^ . 

Thus the spectral average {■)k is equivalent to a phase (or torus) average 
{■)^. We will make extensive use of this equivalence in the following section 
and in Chapter [9l 

If only a subset of bond lengths is incommensurate one may deflne er- 
godic flows on tori of smaller dimension and establish an equivalence of the 
phase average with phase averages over these tori. We refer to the literature 
|25j for details and the application to quantum graphs. 

The equivalence (|6.18|) between spectral statistics and eigenphase statis- 
tics in large graphs on one hand and the equivalence between a spectral 
average and a phase average (j6.5p on the other hand also gives a connection 
between the spectral statistics of an individual quantum system (with the 
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spectrum {kn}) and the spectral statistics of an ensemble of unitary ma- 
trices Usik) = ST{k) 1-^ ST{(p) parameterised by B parameters cpb- Thus 
the spectral statistics of a large graph (with incommensurable bond lengths 
and moderate bond length fluctuations) is equivalent to that of a (special) 
ensemble of unitary random matrices. 

6.4 Spectral correlation functions 

Let us now define spectral correlation functions. We have already shown 
that the spectral average of the density of states reduces to the Weyl term. 
Spectral correlations thus appear in the oscillating part. We define the n- 
point correlation function as [176] 



Rnisi, Sn-i) = A"(d°^^(A;)d°^^(A; + siA) . . . d°^=(A; + s„-iA))fc (6.6) 



where the average is over the whole spectrum < A; < oo. In this defi- 
nition we anticipated that we will be interested mainly in correlations on 
the scale of the mean level spacing A. The corresponding length scale is 
the Heisenberg length Lh = = = 2BL and the corresponding period, 
the Heisenberg period is nn = 2B which is just the dimension of the 
quantum evolution map. 

The two-point correlation function R2{s) and its Fourier transform, the 
spectral form factor 



will be at the centre of much of the discussion in chapters [8] and [H Though 
r and s are dimensionless quantities that have been obtained from rescaling 
a length and a wave number we will often follow the convention to call them 
"time" and "energy" - r = 1 corresponds to the Heisenberg time and s = 1 
to the mean level spacing. 

As mentioned above, spectral correlation functions (and their Fourier 
transforms) are not self- averaging quantities. With the original definition 
of the density of states it is obvious that they are usually a sum over 5- 
functions. An additional average over a small window As ^ 1 in the cor- 
relators or Ar ^ 1 in the form factor has to be performed to smooth the 
function. A combination of a spectral average and a subsequent average 
over r or s will be denoted by {■)s or (•),-. 

Instead of correlations in the spectrum {k^} of the graph one may also 
be interested in the correlations in the eigenphases of the quantum map 
lAsik) averaged over k. The corresponding correlation functions 




(6.7) 




{d°''{<i>; k)d°''{<P + fc) . . . d°'^{ct> + A:))^ (6.8) 
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are defined analogousl}0. For the spectral form factor of the quantum map 
this reduces to 

^ oo 

^(-) = ^EH\-\-^)^n (6.9) 



2B ^ V ' 2BJ 

n=l 

where 

Kn = ^ {\tvl(B{kr\'')k . (6.10) 



The (5-functions in ()6.9p can be smoothed by an additional time average over 
an interval At > 

n \ 1 r+^^/^ 



— }2 = 



2BAt 

n':\n'-n\<BAT 



6.5 Equivalence of spectral correlators and eigen- 
phase correlators in large graphs 

For large graphs with moderate bond length fluctuations spectral and eigen- 
phase correlation functions are equivalent. We will now discuss this state- 
ment in a periodic-orbit approach to the two-point correlation function of 
a quantum graph with incommensurate bond lengths. Replacing cP^'^ by 
the sum over periodic orbits (|5.28|) and performing the spectral average one 
obtains 

p,p' r,r'=l 

(6.12) 

Here, the Kronecker 6{rpLp,rp'Lp') restricts the sum over pairs of periodic 
orbits to pairs of the same (full) length VpLp = Vp/Lp/. This result does 
not depend on the assumption of incommensurate bond lengths. The latter 
has the additional implication that each bond is visited the same number of 
times by a pair of periodic orbits with the same length (thus the full periods 
rpUp and rpiUpi are also the same). Using the notation introduced in Section 
I5.3l the two-point correlator of a quantum graph with incommensurate bond 
lengths can be written as 



The fc-average destroys any dependence on the value which can thus be set to = 
without loss of generality. 
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where {q} is the degeneracy class of periodic orbits of length Lq = X^^^j^ QbLt- 
Interference only takes place between the amplitudes of periodic orbits in 
the same degeneracy class. 

The spectral form factor (for r > 0) 

{<?} ^ ^ P,P'e{q} P P 

is a sum over (5-functions located at the lengths of the periodic orbits - it is 
thus directly related to the length spectrum discussed in Section 15.31 
Finally, 

is the periodic orbit expression for the discrete time form factor obtained 
from the Fourier transform of the two-point eigenphase correlator. Here, the 
sum includes only pairs of orbits of full period n = VpUp = rpiUpi = Ylih=i lb- 
An additional time average over a small interval reveals the equivalence 
of the discrete time form factor Kn and the spectral form factor K(t) for 
large graphs. To this end, let us write the spectral form factor as a sum 

oo 

K{T) = Y,Kn{r) (6.16) 

n=l 

where Kn{T) contains only those periodic orbits of full period n = Upr = 
Upir' . Next, let us discuss the distribution of the quantity Tp = ^rpLpj^Tx 
(this is the metric length in units of the Heisenberg length Lu = 27r/A, 
or the traversal "time") for periodic orbits with period n. If n ^ 1 one 
may invoke the central limit theorem and set Tp = nAL/(27r) + 6Tp where 
nAL/(27r) is the mean traversal time. The fluctuations 6Tp vanish in the 
mean over all periodic orbits and have the variance {6Tp)p = nA^ AL^/ (27r)^ 
where AL^ = {Yli,{Lf, — L)"^) / B is the variance of bond lengths. We consider 
periodic orbits of large period in a large graph n, B ^ oo where ^ is 
constant. Additionally we assume moderate bond length fluctuations in the 
sense that remains constant when B ^ oo. Thus Kn{T) contributes to 
K{t) in a small time window of width ~ ^/njB around r = n/2B. The 
overall contribution of Kn{T) to K{t) can be calculated by the integral 

{Q}-T.b1b=n P,p'e{q} (6 17) 
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As a consequence, with a suitable time average over an interval Ar ~ B~^/'^ 
the spectral form factor is well approximated by the discrete time form 
factor (averaged over periods which correspond to the time interval) for 
large graphs and r ^ 1/2B 

{K{T))r = {Kn)n. (6.18) 

This proves the above statement that the statistics of the spectrum {kn} 
and of the spectrum of eigenphases of the quantum map are equivalent for 
large graphs (if bond length fluctuations are moderate). The discrete time 
form factor does not explicitly depend on the bond lengths. It only 
depends on the degeneracy classes of periodic orbits which have the same 
lengths. All incommensurate choices of the bond lengths lead to the same 
degeneracy classes {q} and thus to the same value for K^. As a consequence 
the spectral form factors {K{T))r for different incommensurate choices of 
the bond lengths are equivelant for large graphs. This is a first signature of 
universality: all large graphs with the same vertex scattering matrices but 
different incommensurate choices of bond lengths share the same spectral 
statistics. 



6.6 The level spacing distribution 
for quantum graphs 

The level-spacing distribution P{s) is one of the most well-known signatures 
of quantum chaos. For a graph with the ordered spectrum {kn}, kn+i > kn 
it is defined by 

A ^ - 

P{s) = lim - V (As - {kn+i - kn)) (6.19) 

A* — s-oo iV ^ — ^ 
n=l 

where 5e{x) is some continuous approximation to a (5-function with width e. 
The limit e — > is implied at the end after N ^ oo and we will omit e in 
the following. 

The level-spacing distribution depends on spectral correlators of all or- 
ders. For this reason it is desirable to express the level spacing distribution 
of a quantum graph directly. Such an expression has been derived by Barra 
and Gaspard [25j and we will outline their approach in this section. 

The starting point is the observation that the spectrum is given by the 
(positive) zeros of a quasi-periodic function f{k). For every graph there 
are many such functions. For example, the secular function f{k) = Cb(^) 
will do the job but C,B{k)a{k) where |a(A;)| > will equally be valid. We 
will not specify the choice of f{k) in the following and the approach can be 
generalised beyond quantum graphs to any spectrum which is given by such 
a function. By definition quasi-periodicity implies that there is a function 
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/((/)) = /((/>!, . . . , 0_b) of B variables (p = i'Pij • • • > 'Pb) which is periodic in 
each argument 

f{<Pi, ...,<Pb + 27r,...,<PB) = f{(pu ...Ab,...,(pB) (6.20) 

such that 

fik) = = Lik, ...,^B = Lek) (6.21) 

where Li, . . . , is a set of rationahy independent frequencies. In the case 
of quantum graphs these frequencies are the bond length^. For the secular 
function Csik) of a graph quasiperiodicity is easily seen since it depends 
on k only through the bond propagator T{k). Replacing Lf,k in 
T{k) I— > T{(f)i^ . . . , (^b) leads to a periodic function C-b(01i • • • > 'Pb) which 
shows that f{k) = Qsik) is a quasiperiodic function. We have already 
encountered this in Section 16.31 The variables (pb live on a B-dimensional 
torus and 

^\(P) = {<Pi + kLu...,(pB + kLk) (6.22) 

defines an ergodic flow in 'time' k on this torus. When the initial condition 
of the flow is set to (j)^'' = ■ ■ ■ = = one has 

/(fc) = /(cI>^(0W)). (6.23) 
Thus the intersections of the flow with the Hyperplane E deflned by 

/(0) = (6.24) 

give the spectrum. Denoting a point on the B — 1-dimensional hypersurface 
S by X the flow defines a Poincare map 

X.+l =HXn) ^g_25) 
kn+l =kn + TRiXn) 

where Xn is a point on S that is mapped by the flow to Xn+i on S. The two 
points are the successive intersections of the trajectory = kLfj with S at 
times kfi and kn^i. The difference T^(Xn) = ^n+i — is the first return 
time to the surface of section. 

We can now write the level spacing distribution as 

A ^ _ 
P(s)= hm ^5{As-TRUn)) 

N~*oo I\ ^ — ' 

n=l 



A ^ 



(6.26) 



N 

n=l 



^ If not all bond lengths of the quantum graphs are rationally independent the number 
of variables has to be reduced to some number smaller than B. The generalisation of the 
following argument is straight forward. However, we will keep a notation that implies 
incommensurability of all bond lengths. 
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As a consequence of ergodicity this distribution is independent of the initial 
conditions (j)^^ = of the flow. Thus, the initial conditions can be taken 
anywhere on the torus without changing the limiting distribution for N 
oo. Moreover, ergodicity implies that the initial condition can be taken 
anywhere on the surface S and that there exists an appropriate measure 
6^~^x ^{x) such that 

P{s)= f d''-^xu{x)A5{As-TR{x{x))) , (6.27) 
Jt. 

where x = (xi, . . . ,xb~i) is a set of i? — 1 variables that parameterise 
Xix) G E. 

We now turn to the calculation of the function i^ix) that determines the 
invariant measure on the hypersurface. A good starting point is the known 
invariant measure on the complete torus. That is for every (measurable) 
function g{(j)i, . . . , 0_b) on the torus ergodicity implies 

^lirn^-i^'' dkgi<^>\<p('^^)) = ^ J d^<Pg{<P) . (6.28) 

where (2^p"d'^</' is the invariant measure on the torus implied by the homo- 
geneous ergodic flow. Defining r[$'^((/)('^^)] as the time of flight after the last 
intersection with the surface of section S (thus T[^^{<j3^^^)] = k — kn the 
last intersection happened at time /c„) and setting (j)^^'' = xJl to the first 
point of intersection we replace the function g{(j)) by 

oo 
n=l 

Let us first show that the left hand side of (j6.28p is the cumulative function 
of the level spacing distribution, that is 

I(s) = / ds' P{s' 
Jo 

= hm -/ dkQ{As-T[^^-^^{xi)])y^5{k-h 

We assume that there are intersections at times k2, ■ ■ ■ , A^at+i of the tra- 
jectory 0*^(xi) with S in the interval < k < K. The integral over k then 
gives 

rK N 



K oo (6.30) 



/ dkg[<f''-''^{xi)] = V Ae(As - (A:„+i - kn)) 



n=l 

N 



(6.31) 



:^Ae(As-T^,[|."-Hxi)]) . 



n=l 



'^This implies a shift fc i-^ fc — fci without changing the time difference between the last 
intersection. Thus, t[$''"'°i (xi)] = fc — fc„. 
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For large values one may replace K = NA and arrives at 



I(s)= lim le{As-TR[^^-\xi)]) 

N^oo iV 



(6.32) 



which shows that /(s) is indeed the cumulative function of the level spacing 
distribution. 

For the calculation of the right hand side of ()6.28p one first needs to 
rewrite the sum over (5-functions such that it depends explicitly on the point 
(p on the torus. Using that the values /c„ are given as the zeros of f{k) and 
thus 

dfik) 



6{k — k„) 



n=l 



dk 



(6.33) 



Now f{k) = /(^'^-^^(xi)) and % = ^Li M-M lead to 



5[$'"'Hxi)] =AQ{As - t[^^-^\xi)])^ 

1)1. 



6=1 



5(/(^'-'Hxi))) 



(6.34) 



As a next step we change the variables (j) ^ (3;, t) where the B — 1 variables 
X = {xi, . . . ,xb-i) parameterise the surface of section S for r = and r is 
the time of flight of the flow since the last intersection with S. That is 



LbT + Sbix) 



(6.35) 



where the B functions satisfy /(si(x), . . . , sb(x)) = such that xi^) ^ 
{si{x), . . . , sb{x)) G S In the new variables we have f[(^^~^^ (xi)) = f{^^ t) 
and a measure 

i^^-'^B-i^^^ Jr^^ (6.36) 
with the Jacobean determinant 



J(x) 





^xdr J{x) 


Li 


... Lb 


dsi 


dsB 


dxi 


dxi 


(9s 1 


(9s s 


dxB-i 


dxB-i 



(6.37) 



The cumulative function becomes 
A 



(27r)^ 
B 



aB~1 



x J{x) 



dre(As -r)x 



(6.38) 



6. Spectral statistics 



57 



The integral over r can be performed explicitly by introducing the new 
variable 

n(T) = />,r) (6.39) 
where x is fixed. With ^ = X^^j^ arrives at 

I{s) = J d^-^x J{x) j due(As - t{u))5{u) . (6.40) 

The integral over u picks up the value of t{u) at u = 0. This has the two 
solutions r = and r = tr(x(x)). Since r is the time of flight since the last 
intersection only the second solution is valid and we finally arrive at 

=7^ / d''-'xJ{x)e(A)e(As-TRixix))) 

_ _ (6.41) 

''-'xJ{x)e{A)e{As-TR{x{x))) 



_ 1 

-/^d^-ixJ(x) is 
where the last line follows from 

A 



(27r)^ 



/ d^-^xJ(x) = 1 . (6.42) 



Equation (j6.42p is consistent with the normalisation of the level spacing 
distribution 

= ^ = /,d^-ixJ(x) ('-"'^ 

and can be derived from 

1 = hm — / dkS" 5{k - kn) (6.44) 

•^'J n=l 

in analogy to the above derivation. The invariant measure on the surface of 
section can now be read off easily as 

d^-^x u{x) = d^-'x J^f . (6.45) 

The integral (16.430 for a large quantum graph is in general too complex to 
be performed analytically. Barra and Gaspard |25j have used this integral to 
investigate the level-spacing distribution of small graphs. For large generic 
graphs approximations have to be used and it is not known how to proceed 
from this expression to understand the universal behaviour of generic large 
graphs that we will discuss in the remaining chapters. 
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Quantum chaos and 
universal spectral statistics 

More than twenty years ago it has been observed that spectral fluctuations in 
individual complex (chaotic) quantum systems are universal |136j . That is, 
all the spectral correlators defined in Chapter [H] (and also the level spacing 
distribution) of any chaotic quantum system are described by non-trivial 
system-independent functions. Once rescaled by the (system dependent) 
mean level spacing, they only depend on some general symmetry properties 
that follow the classification scheme of Wigner and Dyson [2501 \86\ \87\ [88] . 
This symmetry classification scheme is based on the behaviour of the sys- 
tem under a time-reversal operation. Systems either violate time-reversal 
invariance (symmetry class A in the notation of j255j ) or they are time- 
reversal invariant. In the latter case there are two symmetry classes which, 
in fermionic systems, are realised by time-reversal invariant dynamics with 
conserved spin (symmetry class ^I) or with broken spin rotational invari- 
ance (symmetry class All). For a general overview of symmetry classes and 
their algebraic properties, see Appendix |Al The spectral statistics in in- 
tegrable systems is different from the chaotic case. Spectral statistics in 
integrable systems with at least two degrees of freedom is generically Pos- 
sonian [l8]. The above statement of universal spectral statistics is based 
on an overwhelming basis of experimental and numerical evidence and has 
been promoted to a conjecture for chaotic (hyperbolic) Hamiltonian systems 
in the semiclassical regime by Bohigas, Giannoni and Schmit (see also 
[l9l[69]). Since, proving the Bohigas-Giannoni-Schmit conjecture and un- 
derstanding the physical basis of this type of universality has been one of 
the major challenges in the field of Quantum Chaos [1341 1233] 11331 [T93| [5T] . 
The dominant method that has been applied is semiclassical periodic orbit 
theory based on Gutzwiller's trace formula |133l I132j which expresses (the 
oscillatory part of) the density of states as a sum over periodic orbits of the 
corresponding classical Hamiltonian dynamics. In spite of significant recent 
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progress, universal spectral correlations in individual Hamiltonian systems 
are not fully understood. Moreover, some counterexamples of classically 
chaotic quantum systems have been identified which do not show universal 
spectral statistics |148t [56| 1253) . While in these systems the strong devia- 
tion from universality is fairly well understood it raises the question what 
are the precise conditions to find universality in the spectra of classically 
chaotic Hamiltonian systems. Today, large quantum graphs is the only class 
of individual quantum systems for which universality has been proven and 
precise sufficient conditions can be stated |127tll28] . We will give an outline 
of this proof later in Chapter [9l 

7.1 Universal correlators 

and random-matrix theory 

It has been known for some time that universality can be proven with an 
(additional) average over ensembles of systems. The first successful ap- 
proach was random-matrix theory [1761 1130^ I119j where a complex quantum 
system is described by a hermitian matrix with random entries. Random- 
matrix theory has become an important tool to predict spectral statistics, 
wave function statistics and transport statistics in complex quantum sys- 
tems (e.g. disordered mesoscopic systems) [246[ [S^ . Wigner and Dyson 
\250\ [86| [871 EH] proposed three Gaussian random-matrix ensembles, one 
for each symmetry class: the Gaussian Unitary Ensemble (GUE, symme- 
try class ^4), the Gaussian Orthogonal Ensemble (GOE, symmetry class Al) 
and the Gaussian Symplectic Ensemble (GSE, symmetry class All). The 
universal spectral correlators for these ensembles have been calculated an- 
alytically in the limit of infinite matrix dimension |176j . See appendix iBl 
for a definition of the ensembles and for their correlation functions. The 
success of random-matrix theory in predicting universal correlators is re- 
flected in the physical literature where 'universality' is frequently replaced 
by 'random-matrix behaviour'. 

If the spectrum of a quantum system is known (either by experimen- 
tal measurements or by a numerical calculation) a simple way to see, if its 
spectral statistics follows the universal predictions of random-matrix theory 
is to plot a histogram of its level-spacing distribution against the Wigner 
surmises (an energy interval with a few hundred eigenvalues is usually suf- 
ficient). The Wigner surmises ()B.9P are the random- matrix results for the 
level-spacing distributions of the Gaussian random-matrix ensembles of di- 
mension 2x2 (GUE and GOE) or 4 x 4 (GSE). They only deviate very 
slightly from the exact universal result for infinite matrix dimension |134j 
and the difference can only be resolved with very large data sets. In figure 
17.11 the level spacing distribution of a completely uncorrelated spectrum is 
plotted against the Wigner surmise of the GOE. Uncorrelated spectra are 
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realisations of a Poisson process. Poissonian eigenvalue statistics is generic 
for integrable systems (with at least two freedoms) [48^ I134j and the level 
spacing distribution is just an exponential decay 



Pis) 



(7.1) 



The strong deviation of the Wigner surmise from the Poisson distribution 
is a signature of the non-trivial correlations in complex quantum systems. 
The most obvious difference of the two distributions is that small values of 
the level spacing are favoured in a Poisson spectrum (level clustering) but 
are strongly suppressed in the Wigner surmise. This level repulsion in 
complex quantum systems can be understood using a perturbative approach 
|134j . Also large values of the level spacing occur much more frequently in 
a Poisson spectrum than in a universal complex quantum system. That is 
universal spectra are much more rigid than Poisson spectra. 



Pis) 




Figure 7.1: Level spacing distributions of a complex quantum system 
(Wigner surmise for the GOE, full line) and of an uncorrelated spectrum 
(Poisson distribution, dashed line). 



7.2 Analytical approaches to universality 

Apart from random-matrix theory there have been two other successful at- 
tempts to understand universality for ensembles of quantum systems: a par- 
ticle in a random potential (disordered system) [891 1179j and the so-called 
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Pechukas- Yukawa gas where the parametric dependence of eigenvalues and 
eigenvectors on a parameter A in a Hamiltonian of the form H = Hi + Ai?2 
is investigated [611 11991 12521 I134j . In spite of their success and signifi- 
cance these studies (and also random-matrix theory) cannot explain the 
deep quantum-to-classical correspondence between spectral correlations in 
a Hamiltonian systems and the dynamic properties of the underlying classi- 
cal flow. 

This correspondence is the main focus in quantum chaos where semiclas- 
sical periodic-orbit theory has been used [211 [501 [52l EH ESI 1228]. We will 
not go into the details of the semiclassical periodic orbit approach here but 
only summarise qualitatively some of the main results and ideas. For an in- 
troduction we refer to textbooks on Quantum Chaos [1331 11341 12331 1193] and 
to the cited literature. The periodic-orbit expansions for quantum graphs 
which we will present in the next chapter are in most parts analogous to 
semiclassical periodic orbit theory. While periodic-orbit theory on graphs is 
not always trivially generalised to Hamiltonian systems, its main advantage 
is that it allows for exact analytical treatments with controlled approxima- 
tions. 

The semiclassical theory distinguishes between a number of relevant time 
scales which have either classical or quantum mechanical origin. The classi- 
cal scale terg is determined by the inverse Lyapunov exponent. The shortest 
quantum time scale is the Heisenberg time tn = 27r/i/A. This is the time 
scale on which the discreteness of the spectrum is resolved. In the semiclas- 
sical limit in ^ tcrg- 

The spectral form factor has been the main focus of semiclassical periodic 
orbit theory of spectral statistics. The universal spectral form factors for the 
three Wigner-Dyson classes are shown in Fig. 17.21 The non-trivial spectral 
correlations in the form factor are apparent as deviations from x^°^^^"^{t) = 
1. 

Using Gutzwiller's trace formula |1321 1133] the form factor K{t) for a 
classically chaotic Hamiltonian system can be approximated semiclassically 
as a sum over pairs of periodic orbits with period tp = rtn 

p,p 

where Ap/p' are amplitudes corresponding to the periodic orbits and Tp p> are 
repetition numbers. The magnitude of A depends on the classical stability, 
and its phase on the classical action of the orbit p. 

The contribution of periodic-orbits in the range terg ^ tp <^ tn to the 
spectral correlators is quite well understood and it reproduces the universal 
form factor. In a seminal paper Berry |50] showed that a reduction of the 
double sum in (j7.2p to pairs which are either the same or time reversed (the 
diagonal approximation) can be performed on the basis of the sum rules 
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Figure 7.2: The spectral form factors for the three universality classes de- 
scribed by the GUE, GOE, and GSE. 

jl37j and leads to an agreement with the linear universal behaviour of the 
form factor for short times. Since, it has been a challenge to understand the 
next-to-leading order in the short time expansion 

-K'GUE(r) =r 

i^GOE(r) =2r-2r2 + 0(r3) ^^3^ 

i^GSE(r) =^ + ^ + 0(t3) 

of the form factor or to understand the long-time behaviour r > 1. Many 
interesting results on action correlations [2281 [2H [73] in terms of so-called 
'pseudo-orbit' expansions have been obtained [521E5]. Sieber and Richter 
made a crucial progress by identifying pairs of correlated periodic orbits 
whose contribution gives the next-to-leading order term in (17.30 |205[ 12251 
1226] (see also [9l [lOl [H]). This approach has been developed further in 
a number of recent articles [HOl [El EH [232]. The present status 
of semiclassical periodic orbit theory singles out in a systematic way the 
pairs of periodic orbits contribute to progessivly higher terms in the short- 
time expansion of the form factor [1851 1186] . However, due to the essential 
singularity of the form factor at the Heisenberg time r = 1 this result cannot 
be extended beyond r > 1. Another drawback of the present theory is that 
not all periodic orbits are considered. An estimate of the contribution of 
the omitted orbits to the form factor is not given. 
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There exists yet another attempt to explain universahty in individual 
chaotic quantum systems. It relies on a mapping to a supersymmetric field 
theory. It was hoped that these so-called 'ballistic cr-models' would allow a 
direct access to universality via their mean field approximations [U EJ [20} 
I187j . However, there are serious problems in regularising these theories, and 
so far they could only be established under the protection of an additional 
average over (very weak) disorder [183]. 

For individual quantum graphs periodic-orbit methods and methods 
from random-matrix theory and disordered systems like mappings to super- 
symmetric field theories can all be applied rigorously. The inherent disorder 
introduced into the quantum graph by a fixed choice of rationally indepen- 
dent bond lengths is sufficient to treat a quantum graph in a similar way as 
an ensemble of disordered systems or as an ensemble of unitary matrices. In 
the following chapters we will give a detailed account of the periodic-orbit 
theory and its supersymmetric counterpart to quantum graphs. 



Chapter 8 



Periodic orbit theory for 
spectral statistics 

In this chapter we describe the periodic-orbit theory for the spectral form 
factor for large graphs. We will not try to give the most general account 
but will restrict our discussion to star graphs. While the results can be 
generalised to other topologies, the main ideas can be explained in a clearer 
way for star graphs, where one avoids a few technical problems that arise 
in the general case. For the same reason we will only consider the discrete 
version of the form factor which corresponds to two-point correlations in 
the eigenphases of the quantum map. As shown in Section 16.51 eigenphase 
correlations are equivalent to spectral correlations for large graphs with bond 
lengths which are incommensurate and which are confined to a moderately 
narrow interval of lengths. 

The quantum dynamics on a star graph is not affected by adding a 
magnetic field (|2.18|) . The only way to violate time-reversal invariance is by 
choosing a central vertex scattering matrix in an appropriate way. Similarly 
one can add a spin degree of freedom which enables us to study the resulting 
spectral statistics and compare them with the predictions of random-matrix 
theory. 

One of the important features of random-matrix theory is that the statis- 
tics derived for ensembles of N x N matrices tend to well-defined limiting 
distributions in the limit N ^ oo, provided that the spectral parameters 
are properly scaled. In periodic-orbit theory the analogous limit is t// — > oo 
which is the semiclassical limit for Hamiltonian flows. For graphs the limit 
h ^ has no consequence. However, one can make the Heisenberg time tn 
arbitrarily large by increasing the number of bonds B. The disadvantage of 
graphs with respect to Hamiltonian flows is that there is no way to increase 
the number of bonds and at the same time preserve the classical dynam- 
ics on the graph. As a matter of fact there is no unique way to reach the 
limit B ^ oo. In the case of star graphs the limit — > oo is topologically 
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straight forward. However, one has to prescribe a sequence of central vertex 
scattering matrices. For this reason it is more appropriate to consider fixed 
large graphs and expand the spectral functions as powers scries in B~^. In 
the cases of Neumann or DFT vertex scattering matrices which we study 
here as examples one can perform the limit i? ^ oo in a consistent way. 



8.1 Spectral form factor for star graphs 

One of the main technical advantages of star graphs is the fact that one 
may reduce the dimension of the quantum evolution operator from 2B to 
B. This follows from the observation that a wave which scatters into a bond 
is totally reflected at the peripheral vertex. This idea can be implemented 
in the following way. The 2B x 2B graph scattering matrix for star graphs 
has the form 

where is the B x B star scattering matrix. Wc assume only dynam- 
ically connected star graphs for which S^, cannot be brought into block- 
diagonal form by a permutation of bond indices. It is equal to the central 
vertex scattering matrix up to a minus sign which we introduce for con- 
venience. We have chosen Dirichlct boundary conditions at the peripheral 
vertices. These appear as minus the identity matrix in the graph scattering 
matrix. The secular function of a star graph can be written as 

Ceik) = det (las - Usik)) = det {1b - Ue^ik)) (8.2) 

where 

UbM = n{k)S^ ■ (8.3) 
with Ti,{k)i,^y = 6i)h'e'^'^^'' is the B x B star quantum evolution map. The 
star bond propagation matrix T^{k) describes the phase accumulated during 
the propagation from the centre to the peripheral vertices and back. The 
minus signs picked up during the scattering at the peripheral vertices are 
cancelled by the minus sign of the central scattering matrix —S^,. Thus the 
star quantum evolution map UBi,{k) describes the process of scattering at 
the centre and subsequent propagation along the bonds until coming back 
to the centre. This way we have removed the direction index to from the 
description. 

The topoligical length n of a trajectory in this reduced picture is equal 
to the number of peripheral vertices (or undirected bonds) visited. This 
corresponds to a length 2n in the general approach based on directed bonds. 
Thus the Heisenberg period is now nn = B. The discrete form factor at 
time r = ^ is 

K,,, = h\ti S:\\ = Yl S{rpLp,rp,Lp,)-^A^rX;'* (8.4) 

p,p'eV{n) ^ P 
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where a periodic orbit p = bi,b2, ■ ■ ■ , is now specified by the sequence of 
undirected bonds b = 1, . . . , B and its amplitude is the product of the ampli- 
tudes of the star scattering matrix Ap = 3-^^^^^,^^ . . . S^hg^^jS'^hj,?)!- Another 
way of writing the same expression is 



, . ' r) 



pe{q} ^ 



(8.5) 



where {q} is the degenercy class of periodic orbits of length L = 2 Yl^=i QbLb- 
A sufficient condition for a star graph to be invariant under time-reversal 
(symmetry class Al) is that the star scattering matrix is symmetric 

= Si (8.6) 

Otherwise, time-reversal symmetry is in general broken (symmetry class 
A). More precisely, any star graph with a scattering matrix of the form 
5* = SqD where 5*0 is symmetric Sq = Sq and D = diag(e'^'') is time- 
reversal invariant. Under this condition every periodic orbit has the same 
amplitude as its time-reversed partner. Note also, that the secular equations 
for a graph with scattering matrix S^, = SqD and a graph with the scattering 
matrix S'^ = D^^^SqD^^"^ are completely equivalent. 

Time reversal invariant star graphs of symmetry class All can be con- 
structed by adding a spin freedom to the wave function [125J. This will be 
implemented in Section [8. 4[ 

We will use two examples of scattering matrices for star graphs to illus- 
trate our results: i. Neumann star graphs ^381 [391 El US HSl 1149] 

2 

'S'Neumann* 66' = ^^66' " "g ; (8-7) 

and a. DFT star graphs [ESI fT26] 



5'nFT*66' — ^e' B 



1 ■,2nbb' 

DFT* 66' = , 

\/B 

for which the the central vertex scattering matrix is taken as the discrete 
Fourier transform matrix (see Section 13.30 . Both scattering matrices are 
symmetric, thus both families of star graphs are time-reversal invariant. 
The main difference between the two families is revealed by comparing the 
probability of backscattering b ^ b to the probability to be scattered into 
any other bond b ^ b' ^ b. For DFT star graphs an incoming wave 
packet is scattered into any bond with the same probability -Pdft*6«-6 = 
-f*DFT*6^6' = 1/5. For Neumann star graphs with large B back scattering 

-f Neumann* 6^6 = (1 

— -|)^ is favoured, while scattering into any other bond 
has a small probability -PNeumann* 6^6' = In the limit S — > oo this dif- 

ference has a strong impact on spectral statistics. While DFT star graphs 
have the canonical universal spectral statistics of the GOE the Neumann 
star graphs belong to a different universality class [38 ] 141 ] ri49j . 
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8.2 The diagonal approximation 

We shall start the discussion of the periodic-orbit theory by considering the 
first leading term in the short time expansion of the form factor K{t). This 
limit is defined by the requirement that r = n/ B is much smaller than the 
Heisenberg time r ^ r// = 1. 

In (j8.5p we have expressed the form factor K^^n as a sum over degen- 
eracy classes {q} of periodic orbits which share the same length Lq. Each 
degeneracy class {q} gives a contribution 



'b 



\ ^ Ap 



pe{q} P 




E 4^ b (s-i*) 



P,P'6{q}:p7^p' ^ ^ 



If the phases of the amplitudes Ap were completely random one would expect 
that the purely diagonal term which is the first term on the right hand side 
of (|8.9|) dominates the form factor. However, the phases of the terms are not 
random and the extraction of the form factor as a power series in r amounts 
to unraveling the phase correlations in a systematic way. 

Let us now consider the time-reversal invariant case where S^, = S'J. In 
general each periodic orbit has then the same amplitude as its time-reversed 
partner. Therefore time-reverse pairs add coherently to the form factor 
which results in doubling the contribution of the pure diagonal approxima- 
tion. This is not exactly true because there always exist self-retracing 
orbits, for which by definition the time-reversed partners are identical. 

It is natural to devide the form factor into diagonal and off-diagonal 
contributions 

k,^ = kf:^ + kf^-'''^^ (8.10) 

where the diagonal part 

p&V{n) P p&V(n):p=p P 

consists of all pairs p,p' £ V{n) of primitive periodic orbits that are either 
equal p' = p or time-reversed to each other p' = p. The second term is a 
sum over self -retracing periodic orbits which have been counted twice in 
the first term. 

The off-diagonal contribution 

^°n~'"^ = ^ E SirpLp,rp,Lp,)^l;^ (8.12) 

p ,p' aV {n):p' ^p ,p' ^p ' 

to the form factor amounts to all remaining pairs of periodic orbits. It is 
responsible for the quantum interference effects which are not present in the 
diagonal part. 
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The diagonal approximation [50j consists of writing the form factor in 
terms of the diagonal part only. This is justified for short times n <^ B be- 
cause of the following argument. For n <^ B the majority of periodic orbits 
visit any bond at most once. In that case only the orbits which contribute 
to the diagonal approximation contribute coherently to the sum irrespec- 
tively of the chosen star scattering matrix Si, and one may expect that the 
off-diagonal terms are suppressed due to the remaining phase factors. 

At larger times the probability to visit a bond twice is non-negligible - 
this probability grows proportional to the time for n <C -B. When a bond is 
visited twice (or more often) some off-diagonal orbit pairs have correlated 
phases and add up coherently. Eventually these give higher order corrections 
in r = n/B. 

We will come back to a discussion of the off-diagonal part later in Section 
18.31 For the rest of this section we will discuss the diagonal approximation 
to the form factor in detail. Note that 

is just the classical weight of a periodic orbit ()4.24p for the classical evolution 
map 

Mi^w = \Si,wf . (8.14) 

Replacing \Ap\'^ — > Wp in the diagonal part ()8.1ip and comparing it to the 
classical periodic orbit expansion ()4.25p of the return probability Un = trTW" 
one obtains the identity 

p&V{n):rp>2 P p(^V{n):p=p P 

in terms of the classical weights. In the classical periodic orbit expansion 
of the return probability (14.250 repetition numbers Vp > 2 have a different 
factor in front of the weight Wp^ compared to the form factor. This differ- 
ence is neglected in the first term but corrected in the second term of ()8.15p . 
Below we will show that only a very small fraction of orbits p G 'P(n) has 
Tp > 2. The last summand in ()8.15p is the sum over self-retracing periodic 
orbits. 

The first term in ()8.15p is the graph analogue of a general semiclassical 
relation which expresses the short-time form factor as the product of 2t with 
the classical probability to return |84t . 

The total number of periodic orbits of period n grows exponentially 
oc B^/n = e^^^/n with the toplogical entropy At = \nB for a star graph. 
The number of periodic orbits p S V{n) which are repetitions of shorter peri- 
odic orbits, that is rp > 2, is only proportional to B'^/^^rp/n oc e^'^^/^'^rp/n. 
Though this number also grows exponentially it is only an exponentially 
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small fraction of all periodic orbits. For self-retracing orbits a similar ar- 
gumentation is straight forward because one half of the corresponding code 
is just the (time) reverse of the other half. Altogether the total number of 
periodic orbits that contribute in the second and third summand of ()8.15p 
only grows as e^^"/^/n with half the topological entropy. 

Let us now use the classical sum rule (I4.27P to estimate the diagonal 
approximation. The sum rule states ti„ ^ 1 for periods n ^ ?T.mix- On 
the same time scale the second and third part of (|8.15|) which contain the 
contributions from repetitions of shorter orbits and self-retracing orbits are 
exponentially small because the number of these periodic orbits cannot com- 
pensate the small weight of the repeated orbit, which is Wp^ oc (l/B)^ in the 
mean. The largest contributions to the second and third parts come from the 
finite number of short primitive orbits p with a period smaller than the mix- 
ing time rip < n™^. Their repetitions have a weig ht Wp" < e-"/""" which 
can also be neglected. In summary, a mixing classical evolution implies the 
quantum mechanical sum rule 

:^^diag 2 (8.16) 

n 

for the diagonal part of the form factor for time-reversal invariant quantum 
graphs. The deviations die out on the mixing time scale n™^. This time 
scale can be replaced by the (in general shorter) ergodic time scale if an 
additional time average over a small interval An is applied. 

Since the diagonal approximation is only valid for quantum mechanically 
short times n <^ B the sum rule (j8.16p can only be used effectively to predict 
the form factor if there is a time scale separation such that mixing sets in 
much faster than the Heisenberg time n^^^ <C nn = B. In that case there is 
a time regime, n'^^^/B <^ t <^ 1, where the spectral form factor 

Tl ~ 

{K{t = -))r « {K)n « 2t (8.17) 

follows the universal prediction from random-matrix theory up to corrections 
of order 0{t'^). 

In the limit of large graphs 5 ^ oo at a fixed time r = a more 
detailed account is necessary because the number of classical modes (which 
equals the number of bonds) grows. The spectral gap = l/rf"^^ defined 
in (I4.17P may also depend on B. If the spectral gap is bounded from below 
by a constant > (5 > equation (|8.17p becomes exact for B ^ oo for 
every fixed time r. If the spectral gap approaches zero more care has to be 
taken for those classical modes with eigenvalues i/^ — > 1 |235j . In that case 
the corrections 

B 

e=2 



8. Periodic orbit theory for spectral statistics 71 



to the classical sum rule ()4.27p have to be considered in more detail. One 
may estimate these corrections by 

{\un - l|)n <{B- 1)(1 - A,r - B e-^^^ (8.19) 

which shows that the corrections in the limit B ^ oo can only be neglected 
if AgB — > oo. Otherwise deviations from the universal behaviour remain 
in the diagonal approximation for the short-time form factor. Based on 
this observation, Tanner conjectured [235j that large quantum graphs have 
universal spectral statistics (also beyond the diagonal approximation) if the 
spectral gap either remains finite or decays as 

Ag ~ with < a < 1 (8.20) 

in the limit B ^ oo. We will discuss this condition again in connection with 
the supersymmetry approach to quantum graphs in Chapter O We shall 
conclude this general discussion with the following remarks. 

i. If the gap condition is violated, the contributions of repetitions and self- 
retracing orbits ()8.15p to Kn^^ should be examined because there is no a 
priori reason to neglect them. We shall give an example below. 

ii. The gap condition is consistent with the existence of the time regime 1 ^ 
T ^ nP^^/B. It has been shown that the gap conditions holds generically 
for large star graphs [1T | I258j . 

Let us illustrate the general discussion by two examples - the DFT star 
graph and the Neumann star graph. The DFT star graph is known to show 
universal spectral statistics |125j - the classical map is simply MDYi^bV = 
\/B for which all eigenvalues apart from vi = \ vanish exactly. Thus the gap 
condition is fulfilled maximally. The Neumann star graph is an example with 
spectral statistics which belongs to a universality class that is not described 
by the GOE [Ml ESl EE] , here MNeumann* fef,' = (4 + 5w{B^ - ^B))/B^ for 
which all eigenvalues (except for z/i = 1) have the same value = (B — 4)/B 
such that the gap is Ag = 4/ B which does not fulfil the gap condition. Thus 
we also have to consider repetitions of short orbits and self-retracing orbits. 
In the case of Neumann star graph it suffices to consider the repetitions of 
the shortest periodic orbits with primitive period Up = 1. These are at the 
same time self-retracing orbits. Altogether one arrives at |160j 

1 ^ 

^Neumann.n =^^^n + (1 " 2n) J] (-M Neumann. 66 )" 

6=1 

e"^" + 2t (1 - e-^") = 1 - 4t + 0{t^) 

which shows that the contributions from repetitions and self-retracing orbits 
dominate the short-time behaviour. 

We will now study the form factor in the case of a weakly broken time- 
reversal symmetry such that /S* — 7^ is small. The contributions from 



(8.21) 
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time reversed pairs to the diagonal part of the form factor aquire phases and 
they do not add up coherently any more. In that case ()8.15p is replaced by 

k^n-TiirMI + iTU'l) (8.22) 

where 

Tli^bb' = Slf^if^S^^w (8.23) 

and we have neglected the contributions from repetitions and self-retracing 
orbits. The first term in ()8.22p is the contribution of equal pairs of periodic 
orbits. For these the above discussion remains valid and we may replace 
tiM^ — > 1 for n ^ nerg if we assume that the gap condition holds for 
the classical map A^*. The second term in ()8.22p is the contribution from 
time-reversed pairs of periodic orbits. The matrix TZ^, has generally complex 
entries. This reflects the fact that broken time-reversal invariance leads 
to destructive interference of counter-propagating orbits. In the present 
context only the leading eigenvalue the matrix TZ^, has to be considered. This 
can be estimated using perturbation theory. As a result the contribution of 
the time-reversed pairs decays exponentially 

tr7^"cxe~"/""" (8.24) 

in n on a time scale set by 

lK"*«||tr5:(5.-5j)| . (8.25) 

The inverse time scale is an obvious measure for the violation of 

time-reversal invariance. For large graphs S — > oo and fixed time t = n/B 
the contribution of counter-propagating orbits vanishes if 1 /n^ oc B"' where 
< a' < 1. This leaves 

{K{Tf^% = {K^':^)„^T (8.26) 

which is the leading order of the universal result (|7.3p for broken time- 
reversal invariance. 

8.3 OfF-diagonal contributions 

The diagonal approximation only gives the leading order in the short-time 
expansion of the form factor. We will now discuss how universality is built up 
order by order when more and more off-diagonal pairs of periodic orbits are 
included in the expansion of the form factor. We will discuss the limit of large 
star graphs — > oo at a fixed time r = n/B <C 1. The leading off-diagonal 
correction to the diagonal approximation for general large graphs has been 
given by Berkolaiko, Schanz, and Whitney [l2] following the method by 
Sieber and Richter [225^ I226j for Hamiltonian systems. 
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For simplicity we will assume that the elements of the quantum map of 
the star graph are all of the same order \UBik)i,bb'\ ^ B^^^"^ as for example in 
the DFT star graph. As a consequence the gap condition (j8.20p holds with 
a = 0. With some extra effort the following discussion can be generalised 
to graphs where the gap condition holds with < a < 1. 

For time-reversal invariant graphs random-matrix theory predicts a lead- 
ing off-diagonal correction K{t) — K'^^^^{t) = — 2t^ -|- 0{t^) to the diago- 
nal approximation of the form factor. For broken time-reversal symmetry 
random-matrix theory predicts that there are no corrections to the diagonal 
part for short times r < 1. We will start with the time-reversal invariant 
case and come back to broken time-reversal at the end. 

A family of pairs of periodic orbits can conveniently be described in terms 
of diagrams which show the geometry of the two orbits. A diagrammatic 
language has first been introduced in [IJl We will follow closely \125\ 
I126j where a variant which is suitable for star graph dynamics has been 
developed. The art of finding the diagrams which contribute to a given 
order in r to the form factor is significantly simplified by our assumptions. 

To introduce diagrams let us remind that a pair of periodic orbits p and 
p' only contributes to the (discrete) form factor Kn if both orbits have the 
same length Lp = Lpi. Rationally independent bond lengths imply that the 
two orbits visit the same bonds with equal multiplicity in a permuted order. 
That is the code of one orbit is a permutation of the code of the other orbit. 
All periodic orbit pairs for which this permutation is equal (up to a cyclic 
permutation) share a similar geometry which can be expressed as a diagram. 

The diagrammatic language is most easily introduced by first considering 
the n-th trace 

peV{n) 

as a sum over periodic orbits of period n. For a fixed period, say n = 4, one 
may write this sum in terms of the following diagram 

B 

61,62,63,^4=1 

The diagram on the right hand side represents the sum over all periodic 
orbits of period n = 4. The sequence of bonds visited by the periodic orbit in 
the star graph is represented by a sequence of vertices in the diagram (these 
should not be confused with the vertices of the graph) which is traversed 
in the direction indicated by arrow. Each vertex in the diagram on the 
right hand side is translated on the left hand side as one summation index 
b = 1. . . . ,B which is summed over. Each (directed) line b ^ b' is translated 
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as a factor Si^y^. A diagram needs not be labelled with the summation 
indices hj as above. We will do so sometimes to facilitate the translation. In 
the diagrams for the traces s„ with higher period n it is useful to define a 
line element for matrix elements of the n-th power S"" by writing the power 
n explicitly next to the line. This leads to the diagram 



B 

E 

b=l 



{SI 



Ibb 




(8.29) 



which completes the introduction to diag rams for the trace Sn- 

Now we shall introduce diagrams that contribute to the form factor Kn- 
Instead of immediately giving a full set of rules how to draw and translate 
diagrams for the off-diagonal part we will first consider the diagonal part 
where the partner orbits are either the same or time-reversed. Neglecting the 
overcounting of repetitions of shorter orbits and the double counting of self- 
retracing orbits the following two diagrams give the dominant contribution 
to the diagonal approximation of the form factor (here explicitly for n = 4) 



^diag 




4 „ ^ 

/ ^ 'S',»r6ife45'^fe4b35'^ti3ti2'5'itb2fcl('^*6l'>4'^*b4fe3'^*fe3ti2'^*fefcl) ~^ 
bl,b2,b3,b4, = l 

B 

^ ^ Si,bib4 Si^h^h^ Si,b^b2 S-kb2bi{S-kb4bi Si,b3b4 S^^h^b^ Si,bib2 ) 
61,62,63,^4=1 
B 

^ -A^*6i64-A^Vc6463-A^*6362-^*626i • 



4 

B 



61,62,63,64=1 



(8.30) 



The first diagram represents the sum over all equal pairs and the second 
the sum over time-reversed pairs of periodic orbits. As before each vertex 
in the diagram is translated as a summation index bj (the explicit labels in 
the diagram are not necessary). A full line from vertex b to b' is translated 
as the corresponding matrix element and a dashed line as its complex 
conjugate. If a full and a dashed line are parallel (or antiparallel) one may 
translate both lines as one matrix element of the classical evolution map 
M-wb as in the last line (I8.30p . An overall prefactor is added where 1/B 
stems from the definition of the form factor and 4 is the number of cyclic 
permutations of one orbit with respect to the other. 
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In general the two diagrams for the diagonal approximation to the form 
factor Kn^^ has n vertices. The translation rules for lines and vertices are 
the same as in the example above and the prefactor is For drawing the 
diagrams it is useful to replace a long stretch of parallel (antiparralel) full 
and dashed lines by a single pair of parallel (antiparallel) lines and put a 
label next to the line which gives the length of the stretch. For the diagonal 
approximation this leads to the diagrams 



(8.31) 

6=1 

In general a stretch of parallel (antiparallel) lines from vertex b to b' is 
translated as a matrix element {M"')-kb'b of the n-th power of the classical 
map. Ergodicity sets in for n 3> n^'^^ and one may replace (A^")+6';, —>■ 
in the limit n,B oo. The corrections due to classical decaying modes 
vanish exactly in that limit and one obtains the form factor Kn^^ 2t as 
predicted by random-matrix theory. 

We will now introduce diagrams that contribute to the off-diagonal part 
of the form factor. The diagrams have to be ordered by the power of r to 
which they contribute in the short-time expansion of the form factor. A 
lot of care has to be taken when more and more diagrams are added to the 
form factor since certain pairs of periodic orbits occur in more than one 
diagram. Off-diagonal diagrams and the carefull considerations that have to 
be taken when they are added to the form factor are most easily introduced 
by considering the diagram 




(8.32) 



as a pedagogical example. 

We will first describe how this diagram is translated into a sum over 
pairs of periodic orbits. Both periodic orbits in the diagram have period 
n = ni + n2 + 4:. The full line is a periodic orbit with code 6i, . . . , 6„. The 
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dashed line is parallel to the full line along the long stretches of length ni 
and n2 on the left and right. In the code of the partner orbit two indices 
have interchanged bi 6ni+3 and, for this reason, we will call this diagram 
a transposition diagram. The contribution of this diagram to the discrete 
form factor is the sum over all pairs of periodic orbits with the given 
geometry. For ni ^ n2 the diagram is translated as the following sum 




bni+2 bn,+:>, /j,ti+4 

— 




4 E iM:'K+.b,{M':-K^,br.- 

Q* o* c** o* 

6„^+46i bib„^+2 * b2b„^+3 bn^+ab ■ 



(8.33) 

Apart from an additional factor n in the prefactor the translation follows 
the same rules as in the diagonal approximation. The additional factor n 
is explained by the obeservation that interchanging the indices 62 6ni+4 
in the code bi . . . 6.„ describes a different partner orbit than interchanging 
bi ^ &ni+3- The diagrams for both transpositions lead to the same diagram 
- only the labels at the vertices are changed. For ni 7^ n2 all transpositions 
bj ^ with j = 1, . . . , n lead to different partner orbits with the same 

contribution which results in the factor n. For ni = n2 (thus n = 2ni + 4 is 
even) only bj 5j+„^_)_2 with j = 1, . . . , n/2 lead to different partner orbits 
due to the symmetry of the diagram. In that case the additional factor is § 
instead of n. In general the prefactor of any diagram (including the diagonal 

2 

diagrams) is given by ^ for a diagram with s-fold symmetry (e.g. s = n 
for the diagonal diagrams). 

The total contribution of the transposition diagram is now computed 
using the unitarity of the graph scattering matrix 5'* and replacing M2\b' ~^ 

for j = 1, 2. For fixed lengths ni, n2 = n — ni — 4 of the two loops one 
gets 









n2 1 2 



-T^ , (8.34) 



where s = 2 for ni = n2 and else s = 1. 

This shows that a transposition diagram contributes to order in the 
form factor which is the leading correction to the diagonal part. We will 
now add all diagrams of this type to the two diagonal diagrams. Eventually 
it will turn out that their total contribution vanishes once pairs of periodic 
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orbits that occur in more than one diagram have been considered correctly. 
Later we will consider other diagrams that contribute to the same order 
in the form factor and give a non-vanishing contribution which is consistent 
with the predictions of random-matrix theory. 

Taking into account all different transposition diagrams which differ by 
the size of the loops one gets an overall contribution of which di- 

verges as n,B oo 0. This apparent divergence is resolved due to the 
fact that certain pairs of periodic orbits occur in different diagrams and are 
thus double (or multiply) counted. In particular, each of the transposition 
diagrams contains diagonal periodic orbit pairs which have to be subtracted 
before performing the limit of large graphs. In the explicit summation (18.330 
diagonal pairs of periodic orbits appear whenever the two transposed bond 
indices are the same (6i = bm+s)- The subset of pairs of periodic orbits 
which appears both, in the left diagonal diagram in (j8.3ip and in the trans- 
position diagram ()8.33p can be expressed by the diagram 



which has the same value as the transposition diagram ()8.34p . Each trans- 
position diagram is thus cancelled by subtracting the doubly counted orbit 
pairs. As a consequence the transposition diagrams do not contribute at 
all to the form factor in spite of the fact that each diagram is of order 
r^. A transposition of two indices to the right diagram of the diagonal ap- 
proximation in (j8.3ip (thus changing the direction of the dashed arcs in the 
transposition diagram) also leads to new diagrams which are cancelled in the 
same way. This observation can be generalised. Adding a diagram which 
differs from the transposition diagram by a further interchange of two bonds 

2 

inside one of the two long stretches of parallel lines gives a contribution 
which is a factor smaller than the transposition diagram. Subtracting all 
periodic orbits which have been counted twice - once in the transposition 
diagram and once in the new one - cancells the overall contribution of these 
orbits. 

In the transposition diagrams discussed above the two periodic orbits 
are parallel in two long stretches on the left and right of the diagram. We 
will call such long stretches loops because they describe a parallel propa- 
gation of the two trajectories in phase space. On these the motion is free 
of interference and described by the classical map. A long loop of length 

^Transpositions of an element with its neighbour (formaUy ni = —1) or next-neighbour 
(ni = 0) leads to diagrams which look a little bit different but have the same value. 





1 



.2 



= — T 
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rii in a diagram just gives a factor (A^"'^)^^, — > 1/B. Interference occurs 
only in the small region (composed of six bond indices) near the transposed 
indices. It seems natural to consider diagrams with long parallel stretches 
of classical motion with a few regions where quantum interference is active 
in order to find quantum corrections to the diagonal approximation. This 
leads directly to the so-called loop expansion which orders the diagrams 
according to the number of classical loops. The diagonal approximation is 
represented by a single loop. We will show that the universal correction 
to the diagonal approximation is obtained by all diagrams with two loops 
(including, of course, the correction diagrams for multiply counted pairs). 
We have already discussed two types of diagrams with two loops: the trans- 
position diagrams where the periodic orbits traverse the loops in the same 
direction, and the same type of diagram for loops which are traversed in op- 
posite (time-reversed) direction. None of these gave any contribution after 
the doubly counted orbits have been subtracted. There is a third type of 
diagram with two loops. For these one loop contains time-reversed trajec- 
tories and the other parallel trajectories. The diagram is easily drawn and 
calculated 

.^''^^^ 




(8.36) 



This eight-shaped diagram is the analog of the semiclassical Sieber-Richter 
\205\ 12251 [226j pairs for Hamiltonian systems. The left loop has a minimal 
length ni > 2 since the cases ni = 0, 1 have already been included as 
transpositions in the time-reversed diagonal diagram. For the same reason 
the right loop has a minimal length n — ni — 2 > 2. These minimal lengths 
of the loops are the graph analog of the minimal length for a classical loop 
in Hamiltonian flows. Altogether there are n — 5 new diagrams of the same 

2 

value -g^-. Let us now consider the doubly counted orbit pairs which have 
diagrams 




(8.37) 



with n2 = 2, 3, . . . , n — 2. These pairs have either already been included in 
diagonal or transposition diagrams or they appear in two new diagrams with 
different ni . As a consequence we have to subtract n — 3 correction diagrams 

2 

with value Altogether the form- factor in the two-loop expansion is 



Kn = 2T + r'^in - 5 - (n - 3)) ^ 2r - 2t'^ 



(8.38) 
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This is equivalent to the short-time expansion of the random-matrix predic- 
tion for the form factor in time-reversal invariant systems (symmetry class 
^I) in the same order. 

Breaking time-reversal invariance (symmetry class A) effects all anti- 
parallel loops. Instead of the classical evolution map A4 the dynamics along 
the anti-parallel loops is governed by the complex matrix TZ ()8.23p . A di- 
agram which contains an anti-parallel loop is suppressed exponentially and 
does not contribute in the limit B ^ oo. Thus only the first diagonal dia- 
gram remains while all other diagrams in the two-loop expansion are either 
cancelled by doubly counted pairs or vanish due to anti-parallel loops. The 
two-loop result Kn = ^ = r is consistent with the random-matrix predic- 
tion. 

It takes some effort to write down and calculate all three-loop diagrams 
for the next order [H]. The building blocks of the loop expansion are 
all known but it has been performed explicitly only for a few orders. For 
Hamiltonian systems the loop expansion in the semiclassical form factor 
works analogously and has recently been extended to all orders [T85l I186j . 
One should be aware that the loop expansion is valid only for r <C 1. Near 
the Heisenberg-time r = rj^ = 1 the probability that a bond is visited more 
than once is unity and long loops have negligible importance. Similar expan- 
sions which work near or beyond the Heisenberg time are not known. The 
loop expansion also cannot estimate the error due to all neglected diagrams. 
For a fixed n the expansion in diagrams is exact and finite, but in the limit 
i? — > oo, the number of orbit pairs grows exponentially while only a small 
fraction is included in the loop-expansion. 

8.4 Graphs with spin and their spectral statistics 

So far we have only discussed two of the three Wigner-Dyson symmetry 
classes on quantum graphs. The third symmetry class (^11) is relevant for 
fermionic systems with time-reversal invariant dynamics and broken spin- 
rotational symmetry. The time-reversed path of an electron in such a system 
travels along the same trajectory in opposite direction p i— > — p, and has, 
additionally, opposite spin direction s i— > —s. In quantum mechanics the 
time-reversal symmetry is described by an anti-unitary operator T which 
commutes with the Hamiltonian [H, T] = 0. The two Wigner-Dyson symme- 
try classes of time-reversal invariant systems are distinguished by = — 1 
in symmetry class ^11 in contrast to = 1 in symmetry class Al. 

The spectra of a quantum system in symmetry class All is doubly de- 
generate due to Kramers' degeneracy. If is an eigenstate of the system 
Hli') = Elv) the state Tlu) is an orthogonal eigenvector (vlTu) to the 
same energy HT\u) = TH\v) = ET\u). The different type of time-reversal 
symmetry has impact on spectral statistics [218J. The spectral fluctuations 
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which are displayed by complex quantum systems in symmetry class ^11 are 
described by the Gaussian Symplectic Ensemble (GSE). For more details we 
refer to the appendices 1X1 and iB] and to the literature |134j . 

For quantum graphs this symmetry class has first been studied by Bolte 
and Harrison who considered the Dirac-operator on the graph with appropri- 
ate boundary conditions at the vertices and discussed the spin contribution 
to spectral statistics. Quantum graphs with a spin degree of freedom have 
also been discussed in the context of localisation induced by spin-orbit cou- 
pling [37] . Here, we shall introduce spin degrees of freedom which are only 
effected by the scattering at the vertices. As in the previous section we shall 
discuss only star graphs [1251 1126j . The spin degree of freedom is added to 
a graph by considering wave functions ipbaixb) on the bond b which have a 
spin up a = 1/2 and a spin down a = —1/2 component. Spin rotational 
invariance is broken by allowing spin flips at the vertices. As in the previous 
we will limit ourselves to star graphs with spin and allow for spin flips only 
at the centre. 

The construction of the star quantum evolution operator in Section 18.11 
and spectral theory is completely equivalent to usual star graphs when one 
replaces the star graph with B bonds and a two-component wave function 
by a star graph with 2B bonds and a scalar wave function. 

Time-reversal symmetry of the star graph with spin requires that the 
unitary 2B x 2B star graph scattering matric 5* has the structure 

^^=(c Pi ^=-^^ ^=-^^ ^^-^^^ 

where the explicit structure is in the spin index such that the B x B matrix 
A (A^^) describes spin up (down) to spin up (down) scattering from one 
bond to another while B, and C describe scattering processes where the spin 
is changed. That is the scattering amplitude (6', a') {b, a) equals the one 
of the time-reversed process (6, —a) — > (6, —a) up to a minus sign which 
occurs when the spin component is flipped 

S*b(7,b' a' = { — ^ V Si,y ~a' ,b-a- (8.40) 

Since the number of spin flips along a periodic orbit p is even, its amplitude 
in the trace formula equals the amplitude of the time-reversed periodic- 
orbit Ap = Ap where the time-reverse p of an orbit p with the code p = 
(61, (Ti), (62, 172), . . . , {bn, (Tn) is now defined by a reverse order of the bond 
indices bj and flipped spins, p = {bn, — <7n), • • • , (^2) —^'2), (bi, — o"i). 

By convention every pair of degenerate eigenvalues in a spectrum that 
displays Kramers' degeneracy is only counted once in the density of states. 
This leads to an overall factor 1/2 in the trace formula and the mean level 
spacing between different eigenvalues is A = -5= as in a scalar graph with 
the same number of bonds. The equivalence of the spectral form factor with 
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the discrete time form factor of the eigenphases in large graphs leads to the 
relation 

71/ ~ 

{Kir = -))r = {Kn)n (8-41) 
for time-reversal invariant spin star graphs, where 

kr, = ^{\lxUBAkr\'')k- (8.42) 

We will discuss the form factor in the diagrammatic language introduced 
in the previous chapter and assume that the classical evolution operator 
M.i,hcT,b'(j' = \S-kh(j,h'a'\'^ is strongly mixing such that the decaying classical 
modes do not contribute to any of the considered diagrams in the limit of 
large graphs. 

The loop expansion follows the same rules as in the previous section. 
However one needs to add a sum over the 2n spin indices along the two 
orbits. We start with the two diagrams of the diagonal approximation. The 
first diagram can be written in the form 



n 




(8.43) 



where the matrix that appears on the right hand side describes the propa- 
gation of two spins from one bond to another. It contains the classical map 
as a 2B X 2B submatrix which describes the propagation of two parallel 
spins that remain parallel. Anti-parallel spins are propagated by the matrix 

^6(t,6'(t' = -S** 6 a,b' a' K j, (8.44) 

while anti-parallel spins and parallel spins are coupled by 

Jha.Va' = Si,ha,Bi a'S*^ha,h' -cj' Cha ,h' a' = S h a ,B' a' S'l h - a fi' cj' ■ (8.45) 

Being defined in terms of a unitary matrix /C, £ and J all have eigenvalues 
inside the unit circle. Their contribution to the trace is thus exponentially 
suppressed. Only the ergodic mode of the classical map survives in the ra-th 
trace when n is larger than all decay times. In the limit B ^ oo the first 
diagonal diagram gives a contribution r/4 to the form factor. The second 
diagonal diagram gives the same contribution - with the only difference that 
the ergodic modes corresponds to antiparallel spins along the loop. Minus 
signs occur in addition to the classical propagator whenever both spins flip 
but since the number of spin flips is even along a periodic orbit they all 
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cancel. We thus have 



{K{t 



^\diag\ 




+ 




T 

2 



^.46) 



in accordance with the random-matrix result (|7.3p . 

Loops in other diagrams can be treated in a similar way. Each loop 
just gives a factor 1/2B and only transports parallel spins for a loop that 
is traversed in the same direction by both orbits and anti-parallel spins for 
loops that are traversed in opposite direction. For time-reversed loops an 
additional factor —1 appears if the number of spin flips along the loop is 
odd. The corrections to the diagonal approximation are found again in 
the diagrams with two loops. While the transposition diagrams cancel with 
the doubly counted orbits in just the same way as for scalar graphs one 
again gets a non-vanishing contribution from the eight-shaped diagrams of 
the Sieber-Richter type 




E 

^1 1^2 ,''3 i''4,0"l ,0'2 jCa ,0-4 



n2(-l) 



16^2 



^-kh2a2,hl<^l '-'*fe40-4,fe30"3'-'Tlr 63—0-3,610-1 6404,62—02 



n 



8B^ 



iA7) 



The main difference to the scalar case is the different overall sign (the other 

factors are just due to Kramers' degeneracy). The different sign also appears 

2 

in the diagram for the doubly counted orbits which has the same value —^r- 
With n — 5 eight-shaped diagrams (not counting the transposition diagrams) 
and n — 3 diagrams accounting for double counting the two loop correction 

2 2 

is = ^ in accordance with the universal result (|7.3p . 



8.5 Andreev graphs 

and non-standard symmetry classes 

Recently the threefold symmetry classification of quantum systems intro- 
duced by Wigner and Dyson has been extended to a ten-fold classifica- 
tion |247l 12481 \TE[ I255j . This has been necessary for certain physical 
systems, such as a Dirac particle in a random Gauge field, quasiparticles 
in in a hybrid superconducting-normalconducting structure or quasiparti- 
cles in a disordered superconductor. The new feature in these systems is 
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a electron-hole or particle-antiparticle symmetry. The one-particle excita- 
tions are described by a Hamiltonian (for instance, the Dirac operator or the 
Bogoliubov-de Gennes operator) that have a positive spectrum correspond- 
ing to (quasi-) particle excitations and a negative spectrum corresponding to 
antiparticle excitations (while the excitation energies of the full many-body 
field theory remain positive). That is if \v) is an eigenstate with energy 
the particle-antiparticle symmetry yields another eigenstate \v') with energy 
—E. There are two types of such particle-antiparticle symmetry operators. 
They are either unitary or antiunitary. A unitary symmetry operator V 
that anti-commutes with the Hamiltonian \P,H]j^ = is a chiral sym- 
metry operator and an anti-unitary operator C which anti-commutes with 
the Hamiltonian [C, = is a charge- conjugation symmetry. Both 
types of symmetry have impact on spectral statistics near = and lead to 
new quantum interference effects which are not present in the Wigner-Dyson 
classes. Sufficiently far away from the symmetry point E = Q these systems 
are usually described well by the standard Wigner-Dyson classes. Combin- 
ing charge-conjugation and chiral symmetries with time-reversal symmetries 
one can show that there are seven non-standard symmetry classes beyond 
the Wigner-Dyson classes |255) . A full account of these is given in appendix 
El In each non-standard symmetry class random-matrix ensembles can be 
constructed which predict universal statistical properties of their spectra in 
complex quantum systems. 

In this section we will concentrate on two non-standard symmetry classes, 
which have been called C and CI (see appendix |A|), and on chaotic systems 
with spectral statistics which is dominated by the universal predictions of 
Gaussian random-matrix ensembles (see appendix IB]) . The symmetry classes 
C and CI can be realised by quasiparticles in a hybrid superconducting- 
normalconducting structure|14[ll5j. The wave function obeys the Bogoliubov- 
de Gennes equation 

- — 2^ — - - a, (a;) 




.48) 

where ^e(a^) and ^h{x^) are the (quasi-) electron and hole components of 

the wave function, ^— — — — is the kinetic energy of the electron or hole 
in the presence of a magnetic field A, Ep is the Fermi energy, and As(a;) is 
the pair potential of the superconducting condensate which couples electron 
and hole components. The pair potential vanishes in the normal conducting 
region. There is an anti-unitary charge conjugation symmetry C 



which transforms an eigenstate with energy E io a different eigenstate with 
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energy —E. The Bogoliubov-de Gennes equation is time-reversal invariant 
with = 1 (symmetry class CI) if the magnetic field vanishes and the pair 
potential is real (for broken time-reversal invariance the symmetry class is 
C). 

In a normalconducting region electrons and holes are not coupled. In 
a hybrid normalconducting-superconducting system electrons and holes are 
coupled by Andreev scattering (see figure ET]) at a normalconducting-super- 
conducting interface [TH]. An electron (or hole) with energy \E\ <C |As| <C 
Ep in the normalconducting region (A^ = 0) hits the interface and is 
retroflected as a hole (electron) of the same energy in the direction from 
which the incident quasiparticle came (as opposed to normal specular re- 
flection at a potential wall or normalconducting-isolating interface). The 
physical process can be described as an electron which enters the supercon- 
ductor and recombines with an electron-hole pair such that the two electrons 
build a new Cooper pair and the hole is ejected out of the superconductor. 
Up to corrections of order Ej^fEp the momentum is conserved in this pro- 
cess, not only in magnitude but also in direction (again in contrast to normal 
reflection at a potential wall). Note, that the velocity of a hole is in oppo- 
site direction to its momentum due to the overall minus sign in the kinetic 
energy (or 'negative mass') of the hole. 




Figure 8.1: Andreev reflection at a normalconducting-superconducting in- 
terface and normal reflection at a potential wall. 

If a non-magnetic chaotic system is coupled to a superconductor Andreev 
reflection will destroy chaos near the Fermi level [E = 0) since almost every 
trajectory will hit the superconducting interface and the incident particle 
travels back along the same trajectory. Exponential divergence of nearby 
trajectories can only survive on a short time scale between two Andreev 
reflections. Still signatures of chaos in the normal system remain which are 
a topic of interesting present research. For details we refer to the literature 
|29t 11771 [T6| I237t I129j . Here we will only be interested in the regime where 
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the combined electron-hole dynamics is chaoticH ~ this can be achieved by a 
magnetic field which bends both electrons and holes in the same direction 
such that hole trajectories do not travel back along the electron trajecto- 
ries. In general this breaks time-reversal invariance (symmetry class C) but 
certain reflection and point symmetries can restore (a non-conventional) 
time-reversal invariance (symmetry class Cljl. 

The charge conjugation symmetry in a chaotic system effects all spectral 
correlation functions near E = and the deviations from the Wigner-Dyson 
classes decrease on the scale of the mean level spacing A. That is for E ^ A 
chaotic systems in the symmetry classes C and CI are described by the GUE 
and GOE (symmetry classes A and Al). The effects can thus not be seen by 
any type of spectral average in an individual system. Instead the spectral 
average has to be replaced by an average over some system parameter like 
the Fermi energy, the magnetic field or the geometry of the system. The 
deviations from Wigner-Dyson behaviour turn out to be universal and are 
themselves described by Gaussian random matrix ensembles which we will 
call the C-GE and CI-GE. 

In contrast to the Wigner-Dyson symmetry classes there are universal 
interference effects also in the mean density of states for the non-standard 
symmetry classes (averaged over some system parameter). These will be in 
the focus of the following discussion. 



C-GE 



CI-GE 



(<*(«» 1 





Figure 8.2: The mean density of states in the random matrix ensembles C- 
GE and CI-GE (full lines). The dashed lines give the flat density of states 
of a system without charge-conjugation symmetry and the same mean level 
spacing. 



^In the literature on Andreev billiards it has become a convention to call a system 
chaotic or integrable according to the dynamics of the normal system where electrons and 
holes are not coupled. We will always refer to the combined electron-hole dynamics in 
presence of the superconductor. 

^Another option to restore chaos is to introduce many point scatterers (disorder). Both 
symmetry classes can also be realised in a completely different context like two coupled 
spins where chaos in the semiclassical limit (large spins) is not prevented by Andreev 
reflections. 
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In figure [8T2I the mean density of states for the Gaussian random-matrix 
ensembles C-GE and CI-GE [141 [TS] is shown (for explicit formulae see 
appendix [B|) as a function of the energy s = E/AE in units of the mean 
level spacing. Quantum interference leads to a dip in the mean density of 
states at s = where {d{s)) = 0, that is there are no states on the Fermi 
level. For s ^ 1 (or E ^ A) the mean density of states approaches the 
value {d{s)) 1 (or {d{E)) 1/A which defines the mean level spacing 
here). With 



{d{s)) = 1 + 



drF(r) cos 2ttts 



.50) 



we introduce the Fourier transform of {d{s)) — 1 which we will call the 
first- order form factor (or just form factor in this chapter) because its 
periodic-orbit treatment is to some extent analogous to the (second-order) 
form factor in the Wigner-Dyson classes. We will show that the short-time 
expansion 



Firf = - 1 

F{rf' = -l + l + 0{t') 



^.51) 



of this form factor for a quantum graph in the corresponding symmetry class 
is given by a loop expansion analogously to the loop expansion of the second- 
order form factor in the Wigner-Dyson case. In figure 18.31 the universal 
form factors are shown together with a numerical average over Andreev star 
graphs of the corresponding symmetry classes. The form factor is a strongly 
oscillating quantity. To compare it to the universal result an additional time 
average over a small interval At <C 1 is needed to tame the oscillations. 



-I).2- 
-0.4- 
-11.6- 

-n.8- 




Figure 8.3: Universal first-order form factor as predicted by the random 
matrix ensembles C-GE and CI-GE (full lines). The dashed lines are an 
average over 10000 realizations of Andreev star graphs with i? = 50 bonds. 

Andreev graphs jl25l I126j are quantum graphs which are coupled to 
a superconductor such that Andreev reflection occurs. Such graphs can be 
constructed in the symmetry classes C and CI by assigning a two-component 
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wave function '^i,f,-{xh) to each bond where a = e,h corresponds to the 
electron and hole components on the bond b with additional restrictions 
on the boundary conditions at the vertices. We will not give the most 
general definition, instead we will look at star graphs. Each bond of the star 
graph is coupled to a superconductor at the peripheral vertex where Andreev 
reflection couples the electron and hole components. Along the bonds and at 
the central vertex the graph is assumed to be normal conducting (electrons 
are not coupled to holes). 

The quantum evolution map of an Andreev star graph with B bonds in 
symmetry class C (broken time-reversal invariance) has the form 

where the unitary B x B matrix S*^ is the electron-electron scattering matrix 
of the central vertex. Charge conjugation symmetry then requires that the 
hole-hole scattering matrix is the complex conjugate 5"* . The diagonal B x 
B matrices Teh{k) {The{k)) contain the phase factors acquired for a hole 
propagating from the centre to a peripheral vertex where it is scattered 
back as an electron which propagates back to the centre. Altogether 

Tehbb'ik) = -i5,,,e'^''e'2'=^'' n,,,,{k) = -\5,ye-''^e'^'^^\ (8.53) 

where e''^^^'' is the phase acquired duringthe propagating along the bond b 

J4| 



from the centre to the periphery and baclo and — /e (— /e~ *>) is the phase 
acquired during the Andreev reflection for a hole (an electron) hitting the 
superconductor - e"'^*" is the phase of the superconductors pair potential. 
An average over system parameters is easily obtained by taking the phases Ob 
as independent random variables which are equidistributed on < 0;, < 27r. 

For a time-reversal invariant Andreev star graph (symmetry class CI) 
there are additional restrictions: the central electron-electron scattering ma- 
trix (and thus also the hole-hole scattering matrix) has to be symmetric 
Si, = 5j and the pair potential has to be real, thus e'^*" = ±1. For time- 
reversal invariant Andreev graphs the different signs of the pair potentials 
will be taken as independent random variables in a system average. Note, 
that it is sufficient to break either the symmetry of the central scattering 
matrix or the reality of the pair potentials to break time-reversal invariance. 



*At first siglit this does not seem to be in accordance with the Bogoliubov-de Gennes 
equation. Indeed an electron (a hole) propagating along a bond with energy E would 
aquire a phase e'^^^*" {e~ikhLb) according to the Bogoliubov-de Gennes equation where 
electron (hole) momentum is fee = y/Ep + E (fej, = y/Ep — E) - the sign in the phase 
of the hole is different since it propagates in opposite direction to its momentum. Since 
we are interested in the limit Ep 3> E one may expand the momenta fee h = y/ep ± E — 

,3/2 , 



"vAEf ± — j= + 0{ 3/2 )- Adding the hole and electron phases (fee — kh)Lt the leading 
part cancels. Keeping fe — - — f= fixed in the limit Ep 00 only 2kLb remains. 
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Indeed, in the numerical calculations performed in figure [8^ we have taken 
S'jf = Sdft which is symmetric and one finds very good agreement of the 
form factor with the universal result. 

We will assume that the S"^ is a full matrix and the magnitude of all ma- 
trix elements is of order 1 / ^/B. Since we do not perform any spectral average 
there is no fundmental difference between rationally dependent and incom- 
mensurable bond lengths here. For simplicity we choose all bond lengths 
equal Lfc = L. The mean level spacing is A = ii/2BL. Setting k = As in 
the oscillating part of the density of states dP^^{k)dk = dP^'^{s)ds we arrive 
(after a few minor manipulations) at 

oo 

^°''(^) = B cos(2^s|)tr {S.DS:D*r (8.54) 

n=l 

where Dhty = Sbb'^^^*' ■ One can immediately read off the time averaged form 
factor 

{F{t = = {F^)n (8.55) 

where 

F„ = 2(tr iS.DS:D*r)9 (8.56) 

is a discrete variant. The trace appearing in Fn is equivalent to a sum over 
periodic orbits of period 2n which are scattered at the centre of the star 
alternately as electrons or holes. 

Averaging over the phase factors e'^*" (for broken time-reversal invari- 
ance) only those periodic orbits survive in F^ which visit each bond an even 
number of times - one half as electrons being Andreev reflected back to the 
centre as holes the other half as holes being Andreev reflected as electrons. 
The sum over the remaining orbits is still very complex and cannot be per- 
formed in a closed form. It resembles the sum over pairs of periodic orbits 
in the (second-order) spectral form factor where only those pairs survive 
the average where both orbits visit the same bonds. Since there is only 
one orbit there is no diagonal and off-diagonal terms in Fn- However there 
are self-dual periodic orbits which are invariant under electron- hole inter- 
change |123j . That is, a self-dual periodic orbit visits the same sequence of 
bonds twice with electrons and holes interchanged at the second traversal. 
The self-dual approximation takes into account only the coherent contribu- 
tion of self-dual orbits. Since self-dual periodic orbits only exist for odd n 
one has 

^self-dual ^fo for Odd n, 

I — 2tr7W" for even n, 

where Aibb' = is the corresponding classical map on the graph which 

does not distinguish between electrons and holes. We can use the classical 
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sum rule (|4.27p and replace trA^" — > 1 and an average over n yields the 
leading order of the universal term 

^^(^)sclf-dual^^ = (F^'=i^~'^"^')„ = -1. (8.58) 

The overall minus sign is just a consequence of the 2n Andreev reflections 
each of which gives a factor — i. While the diagonal approximation to the 
second-order form factor Kn is linear in r the self-dual approximation starts 
with a constant. The reason for this difference is that the contribution of a 
pair of periodic orbits in a spectral two-point correlator gets an extra factor 
n for all cyclic permutations of the second orbit with respect to the first 
one. In contrast we here deal with a single orbit where such a factor cannot 
arise. 

If time-reversal is broken the contribution of all remaining periodic orbits 
have to vanish for universal spectral statistics. In contrast time-reversal 
invariant Andreev graphs in symmetry class CI display corrections to the 
self-dual approximation - the first order being linear in r. The average in a 
time-reversal invariant Andreev graph is just over e'^*" = ±1. The orbits that 
survive this average outnumber the orbits that contributed at broken time- 
reversal invariance. Every bond is still visited an even number of times - but 
there is no restriction that half of them should come from the centre as an 
electron and the other half as a hole. There are additional orbits already on 
the level of the self-dual approximation. However their contribution cancels 
as we will show now. If an orbit first visits a sequence of n bonds and then 
traverses the same bonds backwards with electrons and holes interchanged it 
is self-dual with respect to charge-conjugation combined with time-reversal 
invariance. This is most easily seen by drawing diagrams 

n n 



self —dual 




+ 




^.59) 



where the first diagram contains the self-dual orbits that we have already 
discussed for broken time-reversal symmetry and the second diagram con- 
tains the additional orbits which only survive in the time-reversal invariant 
case. The first diagram gives a contribution —1 and the second diagram 
gives (— )"2n. The additional factor n stems from the reduced symmetry of 
the diagram with respect to the first self-dual diagram. Averaging over a 
small time interval the contribution of the second diagram cancels due to 
its alternating sign. 
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Higher order corrections in both symmetry classes can be obtained by 
drawing diagrams with more than one loop |126j . One should be aware that 
both types of loops parallel and anti-parallel are present in the time-reversal 
invariant case but only the parallel loops survive the ensemble average in 
the case of broken time-reversal invariance. The linear correction to the 
self-dual approximation is due to diagrams with two loops. We will not 
present a full account of these diagrams but refer to the literature. It has 
been shown that all two-loop diagrams in symmetry class C give a vanishing 
contribution. In contrast there are two-loop diagrams with a non-vanishing 
contribution in symmetry class CI - these have one parallel and one anti- 
parallel loop. Subtracting doubly counted orbits they give the value r/2 
which is the linear order of the universal result. 

We have shown that periodic-orbit theory in Andreev star graphs can 
account for universal spectral statistics in non-standard symmetry classes. 
The remaining five non-standard symmetry classes can equally be treated 
with star graphs but their treatment does not contain any new ideas and we 
refer to the literature [126j . 



Chapter 9 



The Supersymmetry 
approach to Quantum 
Graphs 



In this chapter we wih come back to universal spectral statistics in the stan- 
dard Wigner-Dyson symmetry classes and present a proof that the statistics 
in large well-connected time-reversal invariant quantum graphs follows the 
predictions of random-matrix theory [1271 1128j . This proof involves a field- 
theoretic description of spectral correlation functions in a specially adapted 
version of the supersymmetric non-linear cr-model. While the theory can be 
presented for general simple graphs we restrict again to star graphs which 
allows for some technical simplifications. Similar methods have been used 
efficiently in disordered systems j89l 11791 HSOt 11181 1246| where an average 
over a disorder ensemble leads to supersymmetric cj-models and universal 
spectral statistics can be proven (for an ensemble of systems). The phases 
gifcLj, acquired during the propagation through a bond are a source of dis- 
order in an individual quantum graph with incommensurate bond lengths. 
This type of disorder eventually allows for an exact mapping to a variant of 
the nonlinear cr-model as an additional powerful tool to analyse their spec- 
tral statistics. We will present the theory and the proof of universality for 
the two-point correlator of the eigenphases of the quantum evolution map 
of large star graphs. Generalisations to general graphs can be found in the 
literature [T271[T2H] . 
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9.1 An exact supersymmetric model for spectral 
correlators 



The starting point of the supersymmetry approach to the two-point corre- 
lator of eigenphases of the quantum evolution map is the expression 



Ms) = A 



87r2 dj+dj- 



j±=o 



where 



with 



Cij+,j-;s) 



CB{e+,B;k)CB{0-,B;ky 



(9.1) 



(9.2) 



(9.3) 



^(j+, j„; s) is a generating function expressed as a quotient of secular func- 
tions. A similar expression can be used for the spectral two-point correlator 
R2{s) which we will not discuss here (for moderate bond length fluctua- 
tions we have shown in Chapter [8] that the two correlators are essentially 
equivalent for large graphs). 

The generating function expressed as a quotient of determinants is an 
ideal starting point for supersymmetry approaches since it can easily be writ- 
ten as a Gaussian integral over commuting and anti-commuting variables. 
Powerful methods are available to perform averages on Gaussian integrals 
which renders such an expression a desirable object. We will not give an in- 
troduction to the supersymmetry method and the notions of superdetermi- 
nant sdetA of a supermatrix A or of its supertrace str ^ = tr — tr Af 
which is available in many textbooks and reviews [89l [TM [1801 EsQ, but 
we will give the main steps leading to a convenient Gaussian expression for 
the generating function. 

Defining the supervectors 



ZN 

Xi 
\XN/ 



and Tp = {zl ... z"^ xi ■■■ Xn) , 



(9.4) 



^One should be aware that there are different conventions for the definition of a su- 
pertrace which might differ by an overall sign and consequently for the superdeterminant 
which may be defined via sdet^ = e"*'''"^. There are also different conventions for the 
integration over anti-commuting numbers which differ by an overall constant. 
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where Zi are complex commuting (bosonic) variables while Xi Xi ci-re 
independent anti-commuting (fermionic) numbers, the quotient of deter- 
minants of an X matrix and a (positive definite) N x N matrix 
can be expressed as a Gaussian integral 



det Ai 
det Ai 

Here, 



(sdetA)-^ = J d{il;,^)e-^^^. (9.5) 



is a block-matrix in superspace (boson-fermion space) and the measure is 
given by 

N 

d{i>, tP) = TT-^ H dRe{zi)dlm{zi)dxidxi (9.7) 

i=l 

with / dxi Xi = 1 and / dxi = 0. 

Before applying this to (j9.2p . it will be convenient to double the matrix 
dimensions using 



CB{0;k) =det(l - e-'^^/B.(^) 



(9.S) 



At this point the doubling of dimension seems arbitrary - it leads to simplifi- 
cations at a later stage. Now we write the generating function as a Gaussian 
superintegral 

C{j+,j^;s) = I d(V;,V)e-''^e-^[^'^l (9.9) 

where 



J ( 1 T,{k/2y\ . 



(9.10) 



Here, V' = {^a,s,x,fe} is a 8-B-dimensional supervector where, a = it distin- 
guishes between the retarded and the advanced sector of the theory (com- 
ponents coupling to ^ or respectively). The index s = F,B refers to 
complex commuting and anti-commuting components (determinants in the 
denominator and numerator, respectively), and x = 1,2 to the internal 
structure of the matrix kernel appearing in ()9.10p . The 2x2 matrices 
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are diagonal matrices in superspace containing the appropriate phase factors 
in the boson-boson and fermion-fermion sector. 

To account for the (optional) time-reversal invariance of the scattering 
matrix, we introduce the doublets 

^ (9.12) 

where a^^ = ( o -"i ) is the Pauli matrix in superspace. Notice that the lower 
components of ^ emanate from the upper component by the time reversal 
operation (transposition). For later reference, we note that the new fields 
are related to each other through 

The matrix r is defined by 

r =EBaf' - i£;FCT2^ (9.14) 

where o"*'" are Pauli matrices in the newly introduced 'time-reversal' space 
and -Bb/f the projectors on the bosonic and fermionic sectors, respec- 
tively. All we will need to know to proceed is that r obeys the conditions 

= and = a^^ . (9.15) 

The appearance of the matrix r in (j9.12p suggests to introduce the gener- 
alised matrix transposition 

A^ = tA^t-\ (9.16) 

Using ()9.15|) and (A^)'^ = a^^Aa^^ for a supermatrix [89], one finds that 
the generalised transposition is an involution, 

{A^Y = A- (9-17) 
For later reference we also note that 

= ^'^a^^A^^l = (9.18) 
With all these definitions, the action (j9.1Up now takes the form 



(9.19) 
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where the matrix structure is again in the auxihary index x and we have 
introduced the matrix 

5=(t It) P.20) 

where the matrix structure is in time-reversal space. 

We are now in a position to subject the generating functional to the 
spectral average. We do this by invoking the method discussed in Section 
16.31 whereby the average over k is replaced exactly by phase space averaging 
on the -B-torus. The only fc-dependence is in T^[k/2) which is replaced by 
replace r*({06}) = diag(e"'^''). The averaging is written explicitly as 

{aU,j^;s)U = e--''¥ J d(^,V^)e-^on / (9-21) 

b=l 

Here, 

So =^'+^i^'+_i + + 
is the phase-independent part of the action and 

So far, we have not achieved much other than representing the spectral de- 
terminants by a complicated Gaussian integral, averaged over phase degrees 
of freedom. The most important step in our analysis will now be to subject 
the generating function to an integral transform known as the colour-flavour 
transformation [2541 1256^ I257j . The colour- flavour transformation amounts 
to a replacement of the phase-integral by an integral over a new degree of 
freedom, Z. Much better than the original degrees of freedom, the Z-field 
will be suited to describe the long time behaviour of the system, which is 
equivalent to the low energy sector s ^ 1 of the field theory. 

In a variant adopted to the present context (a single 'colour' and F 
'flavours') the colour-flavour transformation assumes the form 

I d^^,^ei^.++.!:e-V ^ I ^^2,Z) sdet (l - ZZ) ^vlzv-+-^_z.^^ (9 

where 'q± and v± are arbitrary 2F dimensional supervectors and Z, Z are 
2F-dimensional supermatrices. The boson-boson and fermion-fermion block 
of these supermatrices are related by Zbb = -^bB' ■^ff = — ■^fF' "while 
the entries of the fermion-boson and boson-fermion blocks are independent 
anti-commuting integration variables. The integration d{Z,Z) runs over 
all independent matrix elements of Z and Z such that all eigenvalues of 
^BB-^BB ^^^^ than unity and the measure is normalised such that 

j d{Z, Z) sdet(l - ZZ) = 1. (9.25) 



(9.22) 
(9.23) 
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We apply the colour-flavour transformation B times - once for each 
phase (pb- As a result, we obtain a B-fold integral over supermatrices Zb- 
There is one flavour corresponding to the time-reversal index t = 1,2. We 
combine all matrices Zb (Zb) into a single block-diagonal 4i?-dimensional 
supermatrix Z (Z) such that 

Zbts,b't's' = Sb,b'Zbts,t's'- (9.26) 
The averaged generating function now has the form 

(e(j+,j-;^))=e"^ / d(v5,^) / diz,z) 

J J (9.27) 

sdet(l-ZZ)e^^(*'*'^'^) 

where 

Si^,^,Z,Z)=^iQ^ 

^ /z^cSt Z^ \ ^ 

z ._5r^' 



(9.28) 



and we used 2^iZ^'i = ^iZ^'i + ^iZ'^^i, 2^'2-^^2 = '^2Z^2 + ^2Z''^2- 
Here, the indices 1, 2 refer to the auxiliary index x, and the matrix structure 
is in advanced-retarded space. Integrating the Gaussian fields ^' and ^' we 
arrive at the (exact) representation 

{aU,Us)) = I d(Z,Z)e-^(^'^) (9.29) 

where the action is given by 

S{Z, Z) = - str log (1 - ZZ) + ^str log (1 - Z^ Z) 

^ (9.30) 
+ -str log (1 - Sz+Z^z*_S^Z), 



where the prefactor e ' b has cancelled. 



What have we gained with expression (|9.30p apart from an exact refor- 
mulation of the two-point correlator in terms of a complicated supersym- 
metric field theory? This question seems quite urgent since instead of an 
integral over B phase factors e"^*" we now have to deal with an integral over 
B pairs Zb, Zb of supermatrices - each of size 4x4. The main difference 
of expression (j9.3U|) with respect to the definition (|9.2p is the direct cou- 
pling of the retarded and advanced sectors of the theory. In the defining 
expression the retarded sector contributes via (products of) periodic orbits 
of type tr (T^((/))5'>,)" while the advanced sector contributes further factors 
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involving instead. Such periodic orbits gather quasi-random phases. In 
contrast, after the colour- flavour transformation we obtain a theory which 
may be expressed as a sum over periodic orbits str ^Z'^SZS^^ of a quite 
different type where forward and backward scattering alternate. The phase 
of such an orbit can be expected to be less random. We will show in the 
next section how a saddle-point analysis can be used to extract the univer- 
sal contribution to spectral statistics together with sufficient conditions for 
a dominance of universal spectral correlation over small deviations. It is 
not clear if the supersymmetric method may also be useful for graphs which 
deviate strongly from universality like the Neumann star graph. 



9.2 The mean-field approximation 
and universality 

The expressions (I9.29P and (I9.30p for the generating function of the two-point 
correlation function of the eigenphases of the quantum map of a quantum 
graph are an exact identity. While the integral over the modes Z, Z cannot 
be performed analytically in a closed form this expression is an ideal starting 
point for saddle-point analysis for large graphs B ^ oo. We will only be 
interested in the correlator for s <C -B which allows us to expand the sources 
z± defined in (I9.1ip as 

^± = 1 - if {<rrJ± T I) • (9.31) 

Higher orders will vanish in the limit i? — > oo. The resulting action can be 
written as a sum 

S{Z, Z) = Sq{Z, Z) + Ss{Z, Z)s + S+{Z, Z)j+ + Z)j_ (9.32) 

where 

Sq{Z, Z)=- str log [l - ZZ^ + 

+ ^str log (1 - Z^Z) + ^str log (l - SZ^S^Zj 



Ss{Z, Z) = — i— str 



B l-SZ^S^Z (9-33) 



+ ^ ' ^ B 1-SZ^S^Z 
~ .vr SZ-<yfS'^~Z 

S-{Z, Z) =- l-rrStr -rr. 

B 1 _ SZ-^S^Z 



The saddle-point manifold (which is also called the zero-mode for rea- 
sons to become clear below) is equivalent to a mean-field description in 
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which all system-dependent features (encoded in the matrix S) drop out. 



We will show that a reduction of the action (19.330 to the saddle-point man- 
ifold yields an exact expression for universal spectral two-point correlators. 
Later, in Section f9. 31 we will analyse the validity of this approximation and 
give sufficient conditions under which deviations from the mean-field can 
be neglected in the limit B ^ oo. We will assume time-reversal invariance 
Si, = Sj (and thus S = S^, l*'" = 5"^) throughout - the effect of breaking 
time-reversal invariance will be discussed in Section [9.41 

Let us start with the saddle-point equations for the action (I9.33p . Note 
that the source terms j± are the only terms which break the supersymmetry 
- in other words ^(j+ = 0,j_ = 0; s) = 1 which can be seen directly from 
the defining expression (j9.2p . The two-point correlator as a derivative of the 
generating function with respect to the sources is a measure of the response 
of the supersymmetric action at j± = to breaking this symmetry. A saddle- 
point analysis is justified because the largest response can be expected where 
the supersymmetric action is small. This will be shown explicitly when we 
analyse the validity of this approximation. 

The first saddle point equation at j± = and s = takes the form 

^ = = str - str (9.34) 

SZ 1-ZZ 1-Z^Z ^ ' 

and has the solution 

Z = Z^. (9.35) 
The second saddle-point equation reads 

(55o Z SZS^ 

—i^ = = str — — - str ^ 9.36) 

where we have used explicitly the solution Z = Z'^ oi the first saddle-point 
equation and the time-reversal invariance of the system via S'^ = S. It is 
solved for field configurations Z that commute with the scattering matrix 
[5, Z] = 0. This implies equidistribution of the field Z over the bond index 
h. The saddle-point manifold is thus given by the mean-field configurations 

Zobts,t's' = yts,t's' and Zobts,t's' = Yts,t's', (9.37) 

where Y = is a 4 x 4 supermatrix. The commuting parts of Y obey 
= Ibb = ^BB ^"^^ ^FF ~ = ~^FF while the non-commuting 
entries of Y are all independent integration variables. The fermion-fermion 
part is integrated over IR^ ~ while boson-boson part is restricted to the 
compact region where all eigenvalues of ^b^bb are less than unity. 
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Reducing the action (j9.33p to the zero-mode the first contribution van- 
ishes exactly S'o(.^O) -^o) = while the remaining term become^ 

5*""^(y, Y) = -ivrs str ^ + mj+ str ^ - i7rj_ str ■ ^ 



1-YY 1-YY 1-YY 

(9.38) 

Restricting the integration to the saddle-point manifold, we obtain 

(e(j+,j-;«)) -e^°^(j+,j-;«)) - / d(y,y)e-^'^°"(^'^), (9.39) 



where the denotation ^^lOE ^j^dicates that the matrix integral over Y obtains 
but an exact representation of the GOE correlation function. The right hand 
side can be represented in more widely recognisable form. Let us define the 
8x8 supermatrix 



y 1 / \y 1 , , , 

9.40 

1 -h2yy/(i -yy) -2y/(i-yy) \ 
2y/(i-yy) -1 - 2yy/(i - yy)y ' 

where = ( J ) . It is then a straightforward matter to show that the 
action SiY , Y) takes the form of Efetov's action [89j for the GOE correlation 
function 

e^OE(i+,j„;.) = j dQe'«(«) (9.41) 
where the measure is given by 6Q = d(y, y), 

5(g) = |str(Q-S,)e (9.42) 

and ^ = ~ (^"'^'^0 j cr^*'-^ ) ■ "^^'^ ^ discussion of the integral (I9.41|) . and 
the ways random matrix predictions are obtained by integration over Q, we 
refer to the textbook [891. 



9.3 Validity of the mean-field approximation and 
sufficient conditions for universality 

So far, we have shown that the reduction of the exact supersymmetric field 
integral for the generating function ^ to an integral over the saddle-point 



^Note, that for s 7^ or j± 7^ the second saddle-point equation is changed which leaves 
only two saddle points in the reduced action (|9.38p . These can serve as starting points 
for perturbative treatments of the correlator. The loop expansion in the diagrammatic 
periodic-orbit treatment of the form-factor in Chapter [5] is essentially equivalent to an 
expansion around the saddle-point Z = Z = Q. Keeping the full saddle-point manifold of 
So we are able to go beyond such perturbative expansions. 
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manifold (that is over mean-field configurations) results in the GOE spectral 
correlations. However, we have not yet shown under which conditions this 
reduction is actually legitimate. Let us turn to this subject now. 

In a full saddle-point analysis one writes the fields as a sum Z = Zq -|- 
bZ where Zq parameterises the saddle-point manifold and bZ describes the 
fluctuations around the saddle-point manifold. The reduction of the exact 
field integral to the integral over the saddle-point manifold is obtained by an 
expansion of the action S{Z^ Z) in the fiuctuations 6Z to second order and a 
subsequent Gaussian integral over these modes. In general this will lead to 
deviations from the mean-field result presented above. We will investigate 
under which conditions these deviations are small and vanish in the limit 
i? ^ oo of large graphs. 

As we will show later it is sufficient for this purpose to consider the 
expansion of the full action (j9.33p to second order in the fields Z, 



S^^\Z,Z) =sf\z,Z) - ^str (SZ^S^Z) + 



B 



str [afSZ^S^Zj - ^str (sZ^afS^z) , 



(9.43) 



where 



S^^\Z,Z) =str (^ZZ - ^Z^Z - ^SZ^S^Z^ 



1 

=-str 



{Z - Z^){Z - Z^) + Z^Z - SZ'^S^Z 



(9.44) 



Physically, the quadratic action describes the joint propagation of a retarded 
and an advanced Feynman amplitude along the same path in configuration 
space. It thus carries information similar to that obtained from the diagonal 
approximation in the periodic orbit approach. More precisely it is equiv- 
alent to neglecting deviations from Kn^^ = 2ntr A4"/i? which come from 
repetitions and self-retracing orbits. The second order expansion is justified 
if the fluctuations of the fields Z are massively damped (in the sense that 
the matrix elements of Z effectively contributing to the integral are much 
smaller than unity.) Under these conditions, the integration over matrix el- 
ements of Z may be extended to inflnity and we obtain a genuine Gaussian 
integral. In fact, one is forced to extend the integration to inflnity in order 
to preserve ^{j±; s) = 1 at every level of approximation. 

(2) 

The eigenvalues m£ of the quadratic form appearing in Sq determine 
the damping - or the mass, in a field theoretical jargon ~ inhibiting 
fiuctuations of an eigenmode. As indicated by its name, the zero-mode Zq 
carries zero mass. With the observation that the quadratic form involves 
the classical map M-i,bb' = |5'v,;,(,'P of the graph via 

stvSZ^S^Z = str S,bb'Zl,SU,,Zk = str ^ Z,M,,k'Zb' (9.45) 

bb' bb' 
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the action is easily diagonalised by an orthogonal transformation = 
ObiZ'g (we will drop the primes of in the sequel) and can then be 
written as a sum over separate contributions 



B 

:(2) 



S^'Hz,Z) = Y,sf{Ze,Zi) (9.46) 



(2) 

where S), is the contribution 



Sf{Ze, Zi) =-str {Z, - ZJ){Ze - ZJ) + (1 - ui)ZJZi 



vri 

VI— six 



sZlZ, - [UafZl-j_Zlaf)Z, 



(9.47) 



of the eigenmode of the classical map M with eigenvalue ui. Let us now 
focus on the integral over one mode Z^, Z^. We can split the action of a 
single mode into contributions from a pair fields Zi, Zi that satisfies the 
first saddle-point equation Z = Z'^ and contributions violating it. One can 
show that only the configurations Z = Z'^ contribute non-trivially while 
the integral over configurations violating this condition just give a factor 
unitjU. We are left with half the number of integration variables and an 
action = ^-^strZZ = ^siv ZZ for vanishing sources j± = s = 0. In 
this expression we see that the eigenvalues vi of the classical map are related 
to the masses = 1 — of the supersymmetric field theory. Performing 
the remaining integral we will see that the mass mi sets the scale for the 
response to the sources i± which break the supersymmetry of the model. 
That is, a large mass mi corresponds to a strongly damped contribution to 
the two-point correlator. The remaining Gaussian integral over one mode 
gives 



h= j d(Z^,Z^)e- 



ri I ,- 7r(m.<"l)(g+js) i2ri I ■ 7r(mf^l)(s-js) l2 (9.48) 
_ V^^ miB J "I" < miB J 

ri I ■ 7r(m^-l)(s+iA) -|2M i ■ 7r(mf-l)(s-jA) 12 
L miB J L ' rriiB ' 

where j'a = J+ — j- and js = j+ + j- . The generating function in this 



^ For the integral over anti-commuting numbers this can easily seen from the action 
(|9.47|) where the shifts ^^.bf/fb bf — -^Jbf/fb decouple the anticommuting parts of 

Ze from the anti-commuting part of Z. For the commuting entries one can write down the 
integral over all real and imaginary parts of the entries bb/ff and ^^bb/ff ~ shifting 
the entries and changing the contours of integration leads to a similar decoupling such 

where 



that the action at s = j± = is of the form Se — |str ^ cvon^f 
and Ztadd = 



even / 
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approximation is just the product 

B 

^^'\U,j.;s) = l[h (9.49) 
e=i 

over the contributions from each mode. Differentiating with respect to the 
sources j± we finahy obtain the quadratic approximation to the correlation 
function, 



(9.50) 



E 



(m2S2 + 7r2(m£- 1)2^2)2 



The contribution of the zero mode (mi = 0) is given by —:;^t^ and coincides 
with the diagonal approximation to the GOE correlation function. Later on 
we shall see that in the case of broken time reversal invariance, one half of 
the matrix elements of Zq become massive implying that the contribution 
of the zero mode reduces to the GUE expression — 2^^^^ ■ 

In the limit B — > oo, the s-dependence of the contribution of massive 
modes to the correlation function is negligible for our purpose, i.e. individual 
modes contribute maximally as ~ (m£ - lf/2m}B'^ ~ (miB)-'^. Only 
modes of mass ~ i?^", where q is a non-vanishing positive exponent, 
can survive the limit B ^ oo. The contribution of an individual mode 
is negligible if the exponent < a < 1. There are at most 0{B) nearly 
massless modes, and we are led to require that i?2o-i j^^^g^ vanish in the 
limit of large graphs B ^ oo, or that < a < 1/2. 

This is a slightly stronger condition than the one discussed within the 
context of the periodic-orbit analysis of Section [821 where we have seen that 
for fixed r the form factor is universal if < a < 1. The two results are 
fully consistent as the condition < q < 1/2 implies the stronger statement 
that deviations to the two-point correlator vanish uniformly in s. 

In the intermediate region 1/2 < a < 1 - permissible by Tanner's cri- 
terion - non-universal corrections vanish only if the number r of classical 
modes with a small mass remains constant (or does not grow too fast) such 
that S2A2/r oo. If, however, the number of low energy modes is ex- 
tensive, r B, the stricter condition < a < 1/2 has to be imposed. An 
example of a graph with 0{B) almost massless modes is the Neumann star 
graph, for which all the classical modes (apart from the zero-mode mi = 0) 
have a mass m ~ 1/B. The Neumann star thus strongly violates the condi- 
tion for universal GOE statistics - as a single mode with mass m ~ 1/B is 
already sufficient to give non- vanishing corrections. 

Above we have shown that in the limit B ^ oo only the zero mode 
effectively contributes to the correlation function (provided, of course, the 
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master condition Ag ~ is met.) While the zero mode integral must be 
performed rigorously, all other modes are strongly overdamped and may be 
treated in a quadratic approximation. This is the a posteriori justification 
for the quadratic approximation on which our analysis of the mass spectrum 
was based. 



9.4 Breaking time- reversal invariance 

The analysis above applies to time reversal invariant graphs. In this section 
we discuss what happens if time reversal invariance gets gradually broken. 
We assume full universality, i.e. -BA^ S> 1 such that only the zero-mode 
contributes to R2{s). Our aim is to derive a condition for the crossover be- 
tween GOE-statistics in the time-reversal invariant case and GUE-statistics 
for fully broken time-reversal invariance. 

The substructure of the Z-fields with Z = Z'^ in time reversal space is 
given by 

where .^d/C6 Zj^/cb 2x2 supermatrices subject to the constraint 

■^D/CBB = -^c/DBB ^^'^ ^c/UFF = ~-^c/DFF' while the non-commuting 
entries of these matrices arc independent integration variables. The sub- 
scripts D(C) allude to the fact that in disordered fermion systems, the modes 
Zd (Zc) generate the so-called diffusion (cooperon) excitations. Physically, 
the former (latter) describe the interference of two states as they propagate 
along the same path (the same path yet in opposite direction) in config- 
uration space. Cooperon modes are susceptible to time reversal invariant 
breaking perturbations. 

Substituting this representation into the quadratic action, we obtain 

= str (Zd(1 - M^)Zo + Zc(l - 7^*)Zc) (9.52) 

for the action of the zero mode at j± = s = 0. Here, Mi^w = {\iS*b,b'\'^} 
is the classical evolution map and TZ-^w = S-^b,b'Sli^,^. For broken time- 
reversal invariance A^* 7^ TZ^, and the symmetry of the action in time reversal 
invariance space gets lost. 

Noting that B = j,, \S^i,b'\'^, we conclude that the Cooperon zero mode 
Zc,b = Yc acquires a mass term ~ Bmc str(Yc^c)5 where the coefficient 



S*bb'{S^bb' ~ ^*b'b) 



bb' 



(9.53) 
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measures the degree of the breaking of the symmetry S = S'^ . We have 
encountered the Cooperon mass before in Section 18.21 as the inverse time 
riT = 1 /mc on which a pair of time reversed orbits decays in the diagonal 
approximation. The Cooperon mode may be neglected once Bmc ^ oo as 
B ^oo. 

9.5 On higher order correlation functions and the 
gap condition 

The supersymmetry method described above can be generalised in a straight 
forward manner to higher order correlation functions with a generating func- 
tion that contains an appropriate number of spectral determinants in its 
enumerator and denominator. A reduction to a saddle-point manifold gives 
the universal correlations known from random-matrix theory. It all works 
the same way as for the two-point correlator and the condition that this re- 
duction becomes exact in the limit of large graphs also remains unchanged. 
For a finite graph this prooves universal spectral correlators upto small de- 
viations in correlators of order n where n <^ B. 
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Appendix A 



The symmetry classes of 
quantum systems 

Based on earlier ideas of Wigner |250j , Dyson introduced a three- fold classifi- 
cation of quantum systems according to their behaviour under time-reversal, 
spin and rotational invariance [86^ [871 I88j . This symmetry classification 
turned out to be very useful, for instance in semiclassical, disordered and 
random-matrix approaches to complex quantum systems. In this appendix 
we give a short summary of this three-fold symmetry classification and a 
recent extension to a ten-fold classification which is related with the ten 
classes of Riemannian symmetric spaces. We only present a description of 
the ten symmetry classes and refer to the literature |255j for a proof that 
this classification is complete. Table ETTl summarises this appendix. 

In the literature the symmetry classes are often denoted by the Gaussian 
random-matrix ensemble of hermitian matrices respecting the correspond- 
ing symmetry restrictions. Such a notation is misleading and often quite 
confusing. Here, we stick to the notation introduced by Zirnbauer [255j . 

According to a theorem by Wigner a symmetry of a physical system is 
represented in quantum mechanics either by a unitary operator or by an 
anti-unitary operator that commutes (or anti-commutes) with the Hamil- 
ton operator H. Unitary symmetry operators which commute with H are 
associated with constants of the motion. In the eigenbasis of such a uni- 
tary symmetry operator the Hamilton operator is represented by a block 
diagonal matrix - each block connected to one eigenvalue of the symmetry 
operator which serves as a quantum number of the state. In the following 
we will always assume that any unitary symmetry has been used to reduce 
the Hilbert space. 

Anti-unitary symmetry operators which commute with H do not lead 
to a constant of motion, yet they effect the form of the representation of 
the Hamilton operator in a suitable basis. More importantly, they have an 
effect on statistical properties of the spectrum and of the wavefunctions in 
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complex quantum systems. The same is true for (unitary and anti-unitary) 
symmetry operators that anti-commute with H. 

Though every symmetry of a physical system is represented by a uni- 
tary or anti-unitary operator in quantum mechanics the reverse is not true. 
An arbitrary unitary or anti-unitary operator U in general does not qual- 
ify as a symmetry of the system. E.g. for any Hamilton operator H one 
may uniquely define a anti-unitary operator A by its action AJ2n^n\n) = 
^^a*|n) in the eigenbasis H\n) = En\n) of the Hamilton operator. Triv- 
ially, A commutes with the Hamilton operator [A, H] = without being 
connected to any symmetry of the system. In contrast a symmetry operator 
can be defined representation independent which includes that it can be de- 
fined independent of the value one formally assigns to h. This immediately 
rules out the trivial operator A since the eigenbasis will generally depend 
on the value of h. 

A.l Time-reversal invariance 

Consider a physical system represented by the Hamilton operator H. The 
system is time-reversal invariant if an anti-unitary symmetry opera- 
tor T, the (generalised) time-reversal operator exists which commutes 
with the Hamilton operator 

[■H,T] = 0. (Al) 

It can be shown [134^ I125j that such an operator obeys either 

= 1 or = -1. (A.2) 

and T reverses the direction of time in the Schrodinger equation. 

It follows that there are three symmetry classes, the Wigner-Dyson 
classes connected to time-reversal invariance: 

Class A: Systems with broken time-reversal invariance. 
Class ^I: Time-reversal invariant systems and = 1. 
Class ylll: Time-reversal invariant systems and = — 1. 



A particle in a potential described by the Hamilton operator H = ^ + 
V{x) is time-reversal invariant. The corresponding time-reversal operator 
Tai is defined to be the anti-unitary operator obeying Tai x = x and 
TaipT^j^ = —p. In coordinate representation 7^41 is represented as complex 
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conjugation of the wave function {x\Tai\'4!) = {x\iIj)* . Obviously, = 1 
and such systems belong to class Al. 

For particles with spin s = 1/2 and broken spin rotation invariance one 
may define the time-reversal operator TaiI; which acts as Tai on x and p 
and additionally obeys TausT^^ = —s. For this operator T^jj = — 1 and 
the system is in class ^11 if \H, Taii\ = 0. 

In both cases time-reversal is broken by the presence of a magnetic field. 
Thus a particle in a magnetic field generally belongs to class A. Note how- 
ever, that we made a special choice of the time-reversal operator in our 
examples. 

A. 2 Kramers' degeneracy 

The energy spectrum of any physical system which is time-reversal invariant 
with = — 1 is (at least) two-fold degenerate. If |n) (normalised such that 
(n|n) = 1) solves the stationary Schrodinger equation H\n) = En\n) so does 
T\n). Since {Tn\Tn) = {n\n)* = 1 the state T\n) does not vanish. At the 
same time |n) andT|n) are orthogonal (n|Tn) = {Tn\'T'^n)* = — (n|T?7-) = 
such that En is two-fold degenerate. 

A. 3 Chiral symmetries 

and charge conjugation symmetries 

We will now consider certain unitary symmetry operators V and anti-unitary 
operators C which anti-commute with the Hamilton operator, 

[H,V]+ = or [H,C]+ = . (A.3) 

One may have to shift the Hamilton operator by a constant to reveal such a 
kind of symmetry which is tacitly assumed in the following. Equation (jA.3p 
leads to an energy spectrum which is symmetric. If E is in the spectrum, so 
is —E. Since the Hamilton operator of a proper quantum mechanical system 
should be bounded from below such symmetries either occur in systems with 
finite-dimensional Hilbert space (e.g. a finite number of coupled spins) or 
in the framework of one-particle equations of a field theory. Examples are 
the Dirac and Bogoliubov-de-Gennes equations which are formally equiva- 
lent to the Schrodinger equation, that is they can be written in the form 
i^^l^'(t)) = H\^{t)). The negative part of spectrum of these operators 
is related to positive energy excitations of the corresponding field theory. 
Since we will not consider the complete field theory here we will call the 
one-particle operator H the Hamilton operator of the system. 

A chiral symmetry operator V is unitary and obeys V'^ = \. A phys- 
ical system with a chiral symmetry can have an additional time-reversal in- 
variance. In that case the time-reversal operator T and the chiral symmetry 
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operator V are required to commute \T ,V] = 0. Thus there are three chiral 
symmetry classes: 



Class ^III: 



Class BDl: 



Class CII: 



Systems with a chiral symmetry and broken 
time-reversal invariance. 

Time-reversal invariant systems with a chiral 
symmetry and = 1. 

Time-reversal invariant systems with a chiral 
symmetry and = — 1. 



The chiral symmetry classes have first been considered in Quantum 
Chromodynamics as Dirac particles in a random gauge potential [2471 1248j . 
They also occur for Bogoliubov-de-Gennes quasiparticle excitations in the 
certain types of disordered superconductors [IT]. 

Finally a charge conjugation symmetry is described by an anti- 
unitary operator C that anti- commutes with the system Hamilton operator 
and which obeys either = 1 or = — 1. Accounting for the possibility 
of additional time-reversal invariance there are six symmetry classes which 
have a charge conjugation symmetry all of which can be found in certain 
types of disordered superconductors. Two of them have been presented 
as the time-reversal invariant chiral symmetry classes. Indeed, if P is a 
chiral symmetry operator and T with [V,T] = a time-reversal operator 
C = VT is a charge conjugation operator with [C,T] = and = T^. 
The remaining four symmetry classes have not been given any satisfactory 
and consistent name. We will call them charge conjugation symmetry 
classes though the symmetry operator C may not be related to any physi- 
cal charge conjugation (see below). In any case the distinction between the 
three chiral symmetry classes and the four charge conjugation classes has 
historical origins and should not be taken as a mathematical or physical 
relevant distinction. The four charge conjugation symmetry classes are: 
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Systems with broken time-reversal invariance 
and a charge conjugation symmetry obeying 

C2 = -l. 

Time-reversal invariant systems with a charge 
conjugation symmetry obeying both, = 1 
and = -1. 

Systems with broken time-reversal invariance 
and a charge conjugation symmetry obeying 
= 1. This symmetry class is denoted as class 
D if the Hilbert space has even dimension. 

Time-reversal invariant systems with a charge 
conjugation symmetry obeying both, = — 1 
and = 1. 

They have first been discussed in connection to hybrid normalconducting- 
superconducting structures and disordered superconductors [T^ [TB]. 

Note, that the notation for these symmetries stems from the systems for 
which they have first been discussed in detail where they were connected to 
the charge conjugation of a quasiparticle or to the chirality of a Dirac parti- 
cle. In general these symmetries need not have this physical interpretation: 
any of the seven symmetry classes can be realised by a Hamilton opera- 
tor for two coupled spins Si^2 with spin quantum numbers si^2 where the 
(generalised) chiral symmetry operator or (generalised) charge conjugation 
operator are not related to the chirality or electric charge of a particle. 

A. 4 The symmetry classes for quantum graphs 

The symmetries of the Hamilton operator H of a physical system can be 
translated to symmetries for the quantum evolution map [126J Usik) = 
T{k)S{k). It will be more convenient to define an equivalent quantum evo- 
lution map by 

UBi/2{k) = T{k/2)S{k)T{k/2) (A.4) 

and discuss this symmetrised variant. Here, the right factor T{k/2) propa- 
gates the wave function from the centre of a directed bond to the next vertex 
where it is scattered by S{k) to the next directed bond. Finally, the left fac- 
tor T{k/2) propagates the wave function along the next directed bond to 
its center. Some symmetry classes are more naturally discussed for graphs 
with a multi-component wavefunction (spin or particle-hole components). 
A time-reversal operator T acts on the evolution map as 

UBi/2{k) ^ T[UBi/2{k)\ = fUBi/2{kTf^ (A.5) 



Class C: 



Class CI: 



Class BD: 



Class Dili: 
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where T is a unitary matrix with the additional property T*T = ±1. Time- 
reversal invariance of a graph means that the evolution map is equal to its 
time-reversed under the action of T 

UB{k) = T[UB{k)^ = fUB{kff-\ (A.6) 

For symmetry class AI one requires T*T = 1 whereas symmetry class ^411 
requires T*T = —1. The natural choice for the time-reversal symmetry 
operator in symmetry class AI is given by 

r = af" (A.7) 

which acts like the corresponding Pauli-matrix on the direction indices uj 
while bond indices b remain uneffected. Equivalently, the evolution map of 
a time-reversal invariant graph in class AI has the property 

^B(^)l/2a/3 =^b(^)i/2/3o • (^-8) 

Time-reversal invariant graphs with = — 1 (class AIT) can be realised 
as graphs with spin s = 1/2 and 

f = af^ cj^P'" . (A.9) 

Thus time-reversal involves a change in direction and in spin. 

The chiral and charge conjugation symmetry classes on graphs can all 
be realised by either putting the Bogoliubov-de-Gennes or Dirac equations 
on the graph. Thus they always involve electron and hole components of a 
wavefunction (one needs an additional spin a = ±1/2 in some classes). We 
will not give a complete construction of the symmetries for all seven classes, 
for the classes C and CI an example for a star graph is given in Section 18.51 

A chiral symmetry operator V acts on the evolution map as 

Ubi/2{E) ^V[Ubi/2{E)] = VUbii2{-E)V-^ (A.IO) 

where "P is a unitary matrix with V"^ = \ and E is the energy. A quantum 
graph has a chiral symmetry if 

UbiI2{E) =V[Ubi/2{E)^] = VUb{-E)^V~^ (A.ll) 

where taking 1/(b{—E)^ reflects the anti-commuting character of the Hamil- 
ton operator. It is easy to see from equation (jA.lip that the corresponding 
energy spectrum is symmetric. If IAbi/2{E) has an eigenvalue unity and 
there is a chiral symmetry Ubi/2{~E) also has an eigenvalue unity. 

Finally, a charge conjugation symmetry operator C acts on the evolution 
map as 

Ubi/2{E) ^C[Ubi/2{E)\ = CUBi/2{-EYe-^ (A.12) 
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where C is a unitary matrix with = ±1. A quantum graph with a charge 
conjugation symmetry obeys 

Ubi/2{E)=C[Ubi/2{E)^] =CUBi/2{-EfC-\ (A.13) 

Note, that wc have not used the momentum k but the energy E as 
the argument of the evolution map. For the chiral and charge conjugation 
symmetry classes one should be aware that the momentum is in general 
not a constant of motion and thus the energy has to be used. E.g. in 
the Bogoliubov-de-Gennes equation an electron and a hole moving in the 
positive x-direction at energy E are described by the wavefunctions 

^electron _ ^^^JI+E x ^^^^^%x) ~ (^^ g"'^^^ (A. 14) 

where n is the Fermi energy. li fx ^ E one may expand k± = ^/JT^^ ^ 
lb 2^ = fc_F =b 5k. In the limit ^ — >■ oo one may thus replace E ^ 5k 
and the symmetry relations of the evolution map hold for 5k. 
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symmetry class 


T 


V 


c 


symmetric space 


A 











U{N) 


Al 


+1 








U{N)/0(N) 


All 


-1 








U{2N)/Sp{N) 


Alll 





+1 





U{p + q)/U{p)xU{q) 


BDl 


+1 


+1 


+1 


SO{p + q)/SO{p) X SO{q) 


Cll 


-1 


+1 


-1 


Sp{p + q) / Sp{p) X Sp{q) 


C 








-1 


Sp{N) 


CI 


+1 


-1 


-1 


Sp{N)/U{N) 


BD [D) 








+1 


SO{N) 


Dill 


-1 


-1 


+1 


SO{2N)/U{N) 



Table A.l: The ten symmetry classes of quantum systems. If a symmetry 
class obeys time-reversal symmetry or a spectral mirror symmetry the entry 
±1 in the corresponding column indicates if the symmetry operator squares 
to ±1. The entry indicates that the corresponding symmetry is broken. 
The last column gives the corresponding Riemannian symmetric space (of 
compact type). 



Appendix B 



Some relevant results of 
random-matrix theory 

Random matrices were first applied to physical systems by Wigner and 
Dyson |25m [5^1 [571 [55] as a model for complex nuclei. Since, random- 
matrix theory has had an enormous impact on various areas in physics and 
is today one of the main tools to describe statistical properties of complex 
quantum systems. There are many books and reviews on random-matrix 
theory [T30l 11761 1^5 1 [62] and its applications in physics. In this appendix we 
summarise some of relevant results in the context of this review. We will 
neither discuss the various methods to calculate properties of random matri- 
ces nor give any proof of the presented results, and we will only summarise 
results on the Gaussian ensembles of Hermitian matrices. 

B.l Universality and universality classes 

The success of random-matrix theory relies on the fact that many statistical 
properties of the spectrum and of the wavefunctions for a large class of 
complex quantum systems are universal (system independent) yet not trivial. 
When properly scaled they only depend on some general properties of the 
system such as time-reversal invariance. In appendix [A] we have described 
the ten symmetry classes of quantum systems which was completely general. 
Here, we relate these symmetries to the statistical properties of complex 
quantum systems described by a random Hamilton operator. 

A universality class is a subset of a symmetry class which shares the 
same statistical properties - or at least some universal correlation functions 
upto small system dependent deviations. In general there are many univer- 
sality classes within one symmetry class. Each of them can be described 
(and defined) by some ensemble of random matrices in the limit of large 
matrices. Usually there will be a lot of different ensembles that share the 
same universal statistical properties [l"35| . 
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We will focus on the so called ergodic universality classes that can 
be described by Gaussian ensembles of Hermitian matrices in each of the 
ten symmetry classes which describe strongly chaotic systems or disordered 
systems in the delocalised regime. For each of the three Wigner-Dyson 
symmetry classes {A, Al, and ^11) and for two of the charge conjugation 
symmetry classes (C and CI) there is a unique ergodic universality class. 
In contrast for the other symmetry classes there is an additional parameter. 
For instance, the symmetry class BD can be realised in an even or in an 
odd dimensional Hilbert space and this will effect the spectral statistics. 
In the odd-dimensional case the symmetry fixes one eigenvalue at E = 
which changes the statistical properties of the spectrum near E = while in 
the even-dimensional case (also denoted as class D) the charge conjugation 
symmetry does not fix a vanishing eigenvalue. Similarly, a chiral symmetry 
operator may fix some integer number u = 0,1,2 .. . of eigenvalues E = 0. 
The number u is known as the topological quantum number and is important 
for applications of random-matrix theory to quantum chromodynamics. We 
will only present the results for those ergodic universality classes where no 
eigenvalue £^ = is fixed by the symmetry operators V or C (thus 1^ = for 
the chiral classes). 

B.2 The Gaussian ensembles 
of random-matrix theory 

A Gaussian random-matrix ensemble |176j consists N x N Hermitian 
matrices H = H'^ with a Gaussian distribution 

P{H)6H = e-^''^''(M . (B.l) 

This form is the same for all ergodic universality classes. However the sym- 
metries of the Hamilton operator restrict the number of independent matrix 
elements. The flat measure AH = 7\A~^di/iidRe(-ffi2)dIm(ifi2) . . . runs over 
all independent real and imaginary parts of the Hermitian matrix H with 
the additional symmetry requirements of the symmetry class. The real pa- 
rameter A > rescales the spectrum of the system and is a normalisation 
constant. 

If there are no further symmetry restrictions on H = H'^ the ensemble 
is called the Gaussian unitary ensemble (GUE). It describes ergodic 
systems with broken time-reversal invariance in class A. Time-reversal in- 
variance with = 1 (class AX) restricts the matrices H = H'^ to be real 
(and hence symmetric) H = H* = . The corresponding ensemble is 
called the Gaussian orthogonal ensemble (GOE). 

Time-reversal invariant systems in symmetry class All are realised only 
in even dimensional Hilbert spaces due to Kramers' degeneracy. The random 
Hamiltonian is a 2N x 2N matrix restricted by the symmetry condition 
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(^-ij^ V ) ~ (-Iat ^o' ) corresponding ensemble is called 

the Gaussian symplectic ensemble (GSE). Note, that we have made a 
special choice for the time-reversal operators T for the GOE and GSE. 

The Gaussian ensembles for the chiral and charge conjugation symme- 
try classes can easily be constructed from the three Wigner-Dyson ensembles 
GUE, GOE, and GSE by reducing the number of independent matrix ele- 
ments further as required. In the chiral case the ensembles have been called 
chiral Gaussian unitary ensemble (chGUE) for class ^III, the chiral 
Gaussian orthogonal ensemble (chGOE) for class BDl, and the chiral 
Gaussian symplectic ensemble (chGSE) for class CII. The Gaussian 
ensembles for the charge conjugation symmetry classes do not have any 
established name. We will refer to them as C-GE, CI-GE, D-GE, and DIII- 
GE. 



B.3 Spectral statistics for Gaussian random ma- 
trices 

The mean density of states for a random-matrix ensemble is defined as 

{d(s)) = -{tr6(s- H)) = - [ dHP(H)tr6(s-H) (B.2) 
9 9 J 

where g = 2m symmetry classes with Kramers' degeneracy (every eigenvalue 
is counted only once) while g = 1 in other ensembles. For large matrices in 
the Wigner-Dyson ensembles the mean density of states is given by Wigner's 
semicircle [2501 UMl IIZS] law 

(d(^))GUE,GOE,GSE = /l"^ = l + (B.3) 

where we have scaled the spectrum such that {d) = 1 in the centre of the 
semi-circle by an appropriate choice of the scaling parameter A in the defi- 
nition (jB.ip . Here is the size of the matrix (N x N for GUE and GOE, 
2N X 2N for GSE). The result (jR3|) is exact in the limit s,N ^ oo where 
s/N is kept constant. The relevant limit for our purposes is ^ cxd keep- 
ing s constant. In that limit the density of states is flat and the mean level 
spacing. As = {d)~^ = 1. 

In the presence of chiral or charge conjugation symmetries, there are 
deviations from Wigner's semicircle law near s = due to the symmetric 
spectrum. On the scale of the mean level spacing the mean densities of 
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states for these ensembles are given by [HI 11131 12491 11881 1143j 



^2 1 



(c?(s))chGUE =— ^ (^0 (^S) + Jl (^S)) 

(d(s))chGOE =(d('5))chGUE + 2 "^oCtTs) y-~ J '^iMCl 

('i(s))chGSE =(<i(2s))chGUE - T7-'o(27rs) / d^Jo(6 

^ JO 



sin 27rs 

(«!(^)>C-GE =1 - 



(B.4) 



TT 

(d('S)>Cl-GE =(d(s))chGUE - 2'^o(7rs)Ji(vr|s|) 

, , ^ sin27rs 

('i(s))D-GE =1 + 

('i(s))Diii-GE =((i(2s))CT-GE + |ji(2vr|s|) . 

For all seven ensembles the deviation from the flat density of states is pro- 
nounced at s = and decays for s ^ 1. The deviations near s = are 
universal interference effects that can be seen in complex quantum systems 
of the corresponding symmetry class. For quantum graphs this is discussed 
in Section [831 

The spectral two-point correlation function for an ensemble of ma- 
trices is defined as 

s s 
R2{s; So) ={d[so + -)d{sQ --))-! 

I s s (B-5) 

=-(tr 5{so + --H)tT 6(80 ---H))-l. 
9 2 2 

For the Wigner-Dyson ensembles they do not depend on the central energy 
So (in the limit — > cxd; s, sq = const, as before) and are given by [17611134] 



-R2(s)gue =^{s) 



sin^ TTS 



(7r|s| cos TTS - sin7r|s|) (2 Si(7r|s|) - tt) ,^ „x 
-K2(s)goe =-K2(s)gue H 2tt^ 

/ \ \ 27r|s| COS 27rs — sin27r|s| , 
-K2(s)gse =-K2(2s)gueH r^-^ Si(27r s ) . 

The Fourier transform 

/oo 
ds/?2(s)e-2"'^" (B.7) 
-00 

is known as the spectral form factor. For the Wigner-Dyson ensembles 
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(B.8) 



r| for |r| < 1, 

1 for |r| > 1, 

r| (2 -ln(2|T| + 1)) for |r| < 1, 

2 - |r|ln|^ for |r| > 1, 

^(2-ln|l- |t||) for|r|<2, 
1 for |r| > 2. 

Since we do not discuss second-order correlation functions for systems 
with chiral or charge conjugation symmetry we will not state further results. 
Note however, that all spectral correlation function go over to the correlation 
functions of the Wigner-Dyson ensembles for energies s 3> 1 

Another frequently used statistics that depends on correlation-functions 
of any order is the level spacing distribution. Ordering the spectrum such 
that Sn < Sn+i the level spacings the differences Sj = Sj+i — Sj between two 
subsequent eigenvalues. Their distribution for the Wigner-Dyson ensembles 
can be approximated very well by the Wigner surmises 

^^(s)gue = — 5-e 
P(s)GSE=-^e 9.- 



which are exact for 2 x 2 (GUE, GOE) or 4 x 4 (GSE) matrices. 
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